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Vorwort 


Dieses Buch entstand aus einer Sammlung von Programmen, die 
geschrieben wurden, um den vorhandenen Computer zu nützlicheren 
Dingen zu überreden als Strom-, Gas- und Wasserabrechnungen bzw. 
Gewichtsverfolgungen der neugeborenen Tochter. Was lag also näher als 
zu versuchen, den Rechner für Probleme einzusetzen, mit denen man sich 
ansonsten lang und breit (und mit einer gehörigen Portion Frust) «zu 
Fuß» auseinandersetzen mußte; so entstanden Programme zur Aus¬ 
gleichsrechnung, zur Fehleranalyse und zur Funktionsbetrachtung. Da 
das Schreiben dieser Programme gerade für den CPC 464 sehr viel Spaß 
machte und dabei immer wieder neue (gute!) Seiten dieses Computers 
festgestcllt werden konnten, wurde die Programmsammlung ergänzt um 
Problemkreise, die laut Lehrplan bzw. aus eigener Erfahrung in der 
Oberstufe der Gymnasien, in den Fachoberschulen und in den ersten 
Semestern eines Ingenieurstudiums anliegen. 

Demzufolge ist dieses Buch gedacht für alle Schüler ab 11. Klasse, für 
alle Studenten, die sich im Grundstudium mit Mathematik herumschlagen 
müssen, aber auch für alle «fertigen» Techniker und Ingenieure. Es soll 
aber keine reine Programmsammlung darstellen; vielmehr werden die 
mathematischen Grundlagen zu jedem Programm möglichst eindeutig 
erläutert, damit man auch weiß wieso, warum, weshalb. Hierzu gehört die 
Durchrechnung von Beispielen ohne Benutzung des Rechners. Vorausge¬ 
setzt werden in etwa Mathematikkenntnisse der 10. Klasse sowie - in 
einigen Kapiteln - Grundkenntnisse der Differentialrechnung. Des 
weiteren werden vielfach Literaturangaben gemacht, die mit einem G für 
mehr grundlagenorientierte oder mit einem W für weiterführende 
Literatur versehen sind. 

Der Schwerpunkt des Buches liegt aber ganz eindeutig bei den Pro¬ 
grammen selbst; diese können selbstverständlich auch ohne vorheriges 
Durcharbeiten der Grundlagenabschnitte benutzt werden. Außerdem 
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wurden viele REM-Anweisungen eingebaut und auf gute Dialogtechnik 
geachtet (bei dem Speicherplatz des C PC 464 ist das ja kein Problem); die 
Bedienung der Programme ist also recht einfach, und eventuelle 
Änderungen sollten kein Problem sein. Diesem Zweck dienen auch die 
beigestellten Programmbeschreibungen. 

Möge dieses Buch dazu beitragen, das Grauen von der Mathematik ein 
wenig abzubauen. Und wie sagte schon C. F. Gauss: 

Der Mangel an mathematischer Bildung gibt sich durch nichts so 
auffallend zu erkennen wie durch maßlose Schärfe im Zahlenrechnen. 

Daher nicht erst auf die 16-bit-Computergeneration warten (die sind 
sowieso viel zu teuer), sondern ran mit den 8 bits von Schneider! 


Appen 


Harald Baumgart 



1 

Funktionen 


1.1 Darstellung beliebiger Funktionen 

Das Ziel einer jeden Kurvendiskussion ist es, einen prinzipiellen 
Überblick über den Verlauf einer Funktion zu erhalten. Üblicherweise 
wird so vorgegangen, daß besondere Punkte wie Nullstellen, Extrem¬ 
werte, Pole und Asymptoten berechnet werden und als Abschluß der 
Graph der Funktion gezeichnet wird. 

Mit Flilfe eines Computers kann man diese Kurvendiskussion auch um¬ 
gekehrt durchführen. Mathematiklehrer sind darüber zwar nicht sonder¬ 
lich erbaut, aber es erleichtert die Arbeit ungemein, wenn man sich als 
ersten Schritt den Funktionsgraph erzeugt. Das hat zum einen den 
Vorteil, daß man einfach weiß, wie die Kurve verläuft, zum anderen kann 
man für Nullstellen usw. zumindest Schätzwerte entnehmen und so die 
nachfolgenden Rechnungen überprüfen. Für einen Computer spielt die 
Anzahl der Rechenschritte nur eine untergeordnete Rolle. Demzufolge 
geht das Programm so vor, daß die Funktionswerte in sehr kleinen 
Intervallen berechnet und auf dem Bildschirm dargestellt werden. Hierbei 
kann der CPC464 seine hervorragenden grafischen Fähigkeiten voll 
ausspielen. 

In den Zeilen 
Zeilen 10 

bis 270 ist das eigentliche Hauptprogramm abgelegt. Hierbei dienen 

Zeile 20 

und 30 dazu, den Rechner zu veranlassen, bei einem beliebigen 
Fehler, der bei der Rechnung auftritt, nicht mit einer 
Fehlermeldung auszusteigen, sondern den Fehler zu igno¬ 
rieren und weiterzumachen. Fehler bei Berechnungen 
können des öfteren auftreten, da Funktionen durchaus 
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Zeile 


Zeilen 

bis 


Zeilen 

bis 

Zeilen 

bis 

Zeile 


Zeile 

bis 


nicht im gesamten Wertebereich definiert sein müssen. 
Erinnert sei hier an Pole oder Lücken von gebrochen 
rationalen Funktionen, aber auch an die Wurzelfunktion, 
die innerhalb der Menge der reellen Zahlen nicht definiert 
ist, wenn der Wert unter der Wurzel in den negativen 
Bereich übergeht. In 

60 ist die Funktion abgelegt. Diese Zeile wird als Unterpro¬ 
gramm aufgerufen, wenn ein Funktionswert berechnet 
werden soll. 

Da der CPC 464 Schwierigkeiten beim Potenzieren nega¬ 
tiver Zahlen hat, sollte man z.B. x 3 nicht als x A 3, sondern 
als x*x*x eingeben. Auf dieses Problem wird später noch 
eingegangen. 

Wenn eine Funktion auf dem Bildschirm dargestellt werden 
soll, sind bestimmte Grunddaten für das Achsenkreuz wie 
Darstellungsbereich bezüglich x- und r-Achse, Beschrif¬ 
tungsabstand usw. nötig. In den 

80 

110 werden diese Daten per DATA-Anweisung eingegeben. 
Daneben besteht auch die Möglichkeit, diese Daten im 
Dialog einzugeben; die Abfrage, welche Art der Datenein¬ 
gabe gewählt werden soll, liegt in den 

140 

180 Sollen die Daten per DATA-Anweisung eingelesen werden, 
so erfolgt zunächst der Sprung ins Unterprogramm in den 

280 

310 Bevor das eigentliche Lesen durchgeführt wird, muß 
zunächst noch der Zeiger für die READ-Anweisung durch 

290 auf Zeile 100 gesetzt werden, damit auch bei mehrmaligem 
Durchlauf die richtigen Daten verwendet werden. 

Sollen die Daten dagegen im Dialog eingegeben werden, 
springt das Programm ins Unterprogramm 

320 

530 Aus diesen Zeilen kann man auch entnehmen, welche 
Bedeutung die einzelnen Größen haben, da in der READ/ 
DATA-Anweisung bzw. im Dialog natürlich die gleichen 
Variablennamen verwendet werden. An dieser Stelle sei 
außerdem noch darauf hingewiesen, daß an den großen 
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Zeilen 

bis 


Zeile 

bis 


Zeile 

bis 


Zeilen 

bis 


Zeilen 

bis 


Zeile 

bis 


Zeile 

bis 


Strichen der x- und j-Achse, also dort, wo beide Achsen 
grob geteilt werden (dxg und dyg), auch eine Beschriftung 
durch einen Zahlen wert erfolgt; die Skalierung wird hierbei 
automatisch durchgeführt. In den 
540 

770 besteht zunächst die Möglichkeit, sich für den eingegebenen 
Wertebereich der x-Achse den maximalen und den mini¬ 
malen >>-Wert berechnen zu lassen. Wenn der Rechner 
hiermit fertig ist, werden in 
700 

710 beide Werte sowie die entsprechenden Werte für die j-Achse 
ausgegeben, die man ja vorher entweder über READ/DATA 
oder im Dialog dem CPC 464 mitgeteilt hat. In 
720 

760 besteht jetzt die Möglichkeit, die Werte für die j-Achse zu 
ändern; hat man sie im Vorwege zu groß oder zu klein 
gewählt, kann man sie hier noch anpassen, damit die Kurve 
auf dem Bildschirm optimal dargestellt wird. In den 
780 

800 erfolgt die Festlegung des Koordinatenursprungs auf dem 
Bildschirm. Hier wird versucht, diesen so gut wie möglich 
auszunutzen. In den 
870 

920 wird zunächst das Koordinatenkreuz gezeichnet und nach¬ 
folgend die Achsen eingeteilt und beschriftet. Im Unter¬ 
programm 

1490 

1610 werden schließlich die Funktionswerte berechnet und nach 
Anpassung an das Bildschirmformat graphisch ausgegeben. 
Mit Hilfe des letzten Unterprogramms in 

1620 

1800 besteht noch die Möglichkeit, Text in das Diagramm 
einzugeben; dies kann z. B. irgendeine Überschrift, die 
Funktionsbezeichnung, das Datum o.ä. sein. Durch Cur¬ 
sorbewegung besteht noch die Möglichkeit, Text in das 
Diagramm einzugeben; dies kann z.B. irgendeine Über¬ 
schrift, die Funktionsbezeichnung, das Datum o.ä. sein. 
Durch Cursorbewegung wird der Text an eine beliebige 
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Stelle des Bildschirms gebracht und dort ausgedruckt. Der 
Vorgang kann beliebig oft wiederholt werden; auf diese 
Weise kann man auch z.B. Achsenbezeichnungen löschen 
oder ändern. 

Den Abschluß des Programms bildet 
Zeile 270 Durch sie verschwindet der Cursor. Beendet wird das 
Programm durch zweimaliges Drücken der ESC-Taste. 
Ohne auf jede einzelne Zeile hinzuweisen, sollen an dieser 
Stelle noch einige Standard-Programmierweisen bespro¬ 
chen werden, wie sie auch in allen anderen Programmen 
dieses Buches Vorkommen. 

Die erste sei die Eingabe einer Codezahl, aufgrund der das 
Programm Entscheidungen trifft und z.B. in bestimmte 
Unterprogramme verzweigt. Im Programm "Funktion 1.0" 
kommt das in den 

Zeilen 150 

bis 180 vor; sehr viel ausgeprägter ist es in den Menüprogrammen 
ab Kapitel 4. Die Codezahl wird mit einer Kurzbeschrei¬ 
bung ihrer Funktion ausgegeben, diese Aufgabe erfüllen 

Zeile 150 
und 160 

Wegen der INPUT-Anweisung in 

Zeile 170 wartet das Programm dann auf die Eingabe der Codezahl; 
dies erfolgt hier über die Variable a. 

Gleich nach Eingabe der Codezahl wird geprüft, ob sich a 
innerhalb des vorgegebenen Wertebereiches befindet, da 
ansonsten eine Fehleingabe vorliegt. In diesem Fall wird der 
gesamte Text erneut ausgegeben, so daß eine neue Eingabe 
erfolgen muß. Wird statt einer Zahl z. B. ein Buchstabe oder 
ein sonstiges Zeichen eingetippt, erkennt der Rechner diesen 
Fehler von sich aus und verlangt die Eingabe einer Zahl. 
Wenn diese Prüfung durchgeführt worden ist, kann 
Zeile 180 in die entsprechenden Unterprogramme verzweigen. 

Die zweite besprochene Standardprogrammierweise sei eine 
"Ja/Nein-Abfrage"; diese kommt z.B. in den 

Zeilen 1750 
bis 1780 vor. 

Auch hier wird zunächst ein Text ausgegeben, aufgrund 
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dessen die Ja/Nein-Entscheidung ausgeführt wird. Das 
Programm erwartet die Eingabe von a%, also einer String¬ 
variablen. In diesem Buch wird generell die Ja-Anwort 
durch Betätigung einer der beiden ENTER-Tasten gegeben. 
Dies hat den Grund, daß zumindest die eine ENTER-Taste 
beim CPC 464 erfreulich groß ausgefallen ist, so daß eine 
Betätigung besonders komfortabel ist. Im allgemeinen 
kommt die Ja-Antwort wesentlich häufiger vor als die Nein- 
Antwort; diese gibt man durch Betätigen der n-Taste. 
Falls die ENTER-Taste gedrückt wurde, springt 
Zeile 1760 zur entsprechenden Programmstelle; das gleiche macht 
Zeile 1770 für die n-Taste. Wichtig ist hier, daß beide möglichen 
Antworten getrennt abgefragt werden. Wurde nämlich 
versehentlich ein anderes Zeichen eingegeben, gelangt das 
Programm wegen der vorhergehenden Abfragen nach 
Zeile 1780 und die Eingabe wird nicht akzeptiert, sondern der 
Abfragetext wird erneut ausgegeben. Fehleingaben sind 
somit ausgeschlossen. 

Um die Arbeitsweise des Programms zu zeigen, ist in Zeile 60 die 
Funktion 

V = 3jc 3 - 10x 2 + 3 

abgelegt. Wählt man die Darstellung über die DATA-Anweisung, so wird 
die Funktion im Wertebereich — 2 bis +4 wiedergegeben. Die y-Achse 
geht von - 15 bis + 5. 

Länge des Programms; etwa 6,7 K. 
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10 REM Programm Funktion 1.0 
20 ON ERROR GOTO 30 
30 RESUME NEXT 
40 GOTO 120 

50 REM Zelle 60 = reserviert fuer die Funktlonsglelc 
hung 

60 y = 3*x*x*x-10*x*x+3 
70 RETURN 

80 REM 100=DATA fuer die Grunddaten Koordinatenkreuz 

90 REM Ablegeweise: DATA xmin,xmax,ymln,ymax,dxg,dxk 

,dyg,dyk,darstx,darsty 

100 DATA -2,4,-15,5,6,2,4,5, 1,3 

110 REM Moegllche Werte fuer darstx bzw. darsty: 0 = 
keine Darstellung » 1 = IDezimalste1le »2=2 Dezi 

malstellen > 3 = keine Dezlm 
alstelle 
120 MODE 2 

130 PRINT'Programm Funktion 1.0 zur Darstellung eine 
r beliebigen Funktion ":PRINT"ln einem Koordlnatenkr 
e uz " 

140 LOCATE 1,10 :PRINT“Llegen die Grunddaten fuer da 
s Achsenkreuz <xmin,xmax etc.) dem Programm ln e 

inerDATA-Anweisung vor oder 
sollen sie im Dialog eingegeben werden ?" 

150 PRINT: PRINT TAB<20)J"DATA";TAB<50);"<1)" 

160 PRINT:PRINT TAB<20);"Dialog“JTAB(50)»"<2) " 

170 INPUT a : IF a<l OR a>2 THEN 140 
180 ON a GOSUB 280,320 

190 GOSUB 540 : REM Berechnung und Ausgabe von ymln 
und ymax 
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200 GOSUB 780 :REM Berechnung der Bildschirmpunkte f 
uer den Koordinatenursprung 

210 GOSUB 870 : REM Zeichnung des Koordinatenkreuzes 

220 GOSUB 930 : REM Einteilung ln dxg,dyg 

230 GOSUB 1030 : REM Einteilung ln dxk,dyk 

240 GOSUB 1130 : REM Beschriftung der Achsen 

250 GOSUB 1490 : REM Berechnung und Ausgabe der Funk 

t lonswerte 

260 GOSUB 1620 : REM Texteintragung ins Diagramm 
270 GOTO 270 

280 REM UP Lesen der Grunddaten Koordinatensystem pe 
r DATA 

290 RESTORE 100 

300 READ xmln,xmax,ymln,ymax,dxg,dxk,dyg,dyk,darstx, 

darsty 

310 RETURN 

320 REM UP Eingabe der Grunddaten Koordinatensystem 
lm dlalog 

330 CLS : PRINT’MIn-Wert der X-Achse : ’5:INPUT xmin 

340 PRINT:PRINT"Max-Wert der X-Achse : ";:INPUT xmax 

350 PRINT:PRINT"Min-Wert der Y-Achse : ";:INPUT ymin 

360 PRINT:PRINT"Max-Wert der Y-Achse : ";:INPUT ymax 

370 PRINT:PRINT:PRINT:PRINT"Wie oft soll die X-Achse 
grob (mit Beschriftung) geteilt werden : *;:IN 
PUT dxg 

380 PRINT:PRINT’Wle oft soll ein Feld noch fein 
(ohne Beschriftung) unterteilt werden : “:INPUT dxk 

390 PRINT:PRINT"Wie oft soll die Y-Achse grob ( 
mit Beschriftung) geteilt werden : ";:INPUT dyg 
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400 PRINT:PRINT"Wie oft soll ein Feld noch fein 
(ohne Beschriftung) unterteilt werden : INPUT dyk 

410 CLS:PRINT’Wle sollen die Zahlen an der x-Achse d 
argestellt werden ?*: 

420 PRINT:PRINT:PRINT TAB<20);"-gar keine Darstellun 
g";TAB<60);"(0)" 

430 PRINT:PRINT TAB<20);"-m1t einer Dezimalstelle"^ 
AB(60){* <1)" 

440 PRINT:PRINT TAB<20);"-mit zwei Dezimalstellen";T 
AB(60)*"< 2) “ 

450 PRINT:PRINT TAB<20) i "-ohne Dez 1 malste1 len" ; TAB(6 
0);"(3)" 

460 INPUT darstx : IF darstx<0 OR darstx>3 THEN 410 
470 CLS:PRINT'Wie sollen die Zahlen an der y-Achse d 
argestellt werden ?" 

480 PRINT:PRINT:PRINT TAB<20);"-gar keine Darstellun 
g“;TAB(60)»"(0)" 

490 PRINT:PRINT TAB(20);"-mit einer Dezlmalstelle";T 
AB(60>;"(1>" 

500 PRINT:PRINT TAB<20);“-m11 zwei Dezimalstellen"^ 
AB<60);"(2)" 

510 PRINT:PRINT TAB(20);"-ohne Dezimalste1 len";TAB<6 
0);"(3)" 

520 INPUT darsty : IF darsty<0 OR darsty>3 THEN 470 
530 RETURN 

540 REM UP Berechnung und Ausgabe von ymin und ymax 
550 CLS : PRINT'Sollen ymln und ymax berechnet werde 
n ( /n) :INPUT a$ 

560 IF a*=“" THEN 590 
570 IF a$="n" THEN 770 
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580 GOTO 540 
590 a=- 1 

600 PRINT:PRINT:PRINT"Bltte einen Augenblick Geduld. 
Ich arbe1te . “ 

610 PRINT:PRINT"Wenn Ich bei 500 bin, bin Ich fertig 
620 ymlnb=0 : ymaxb=0 

630 FOR x=xmin TO xmax STEP (xmax-xm1n)/500 
640 a=a+l 

650 LOCATE 1,23 : PRINT a 
660 GOSUB 60 

670 IF y<yminb THEN yminb=y 
680 IF y>ymaxb THEN ymaxb=y 
690 NEXT x 

700 CLS : PRINT“Blslang Im Programm verwendete Werte 
fuer die Einteilung der y-Achse : ":PRINT:PRINT"yml 
n = ";ymln,"ymax = "jymax 

710 PRINT:PRINT:PRINT"Soeben aus der Funkt 1onsgleIch 
ung berechnete Werte PRINT:PRINT"ym1n = "Jyminb," 

ymax = "»ymaxb 

720 PRINT:PRINT:PRINT"Sol len die Werte der Y-Achsene 
lntellung veraendert werden < /n> ";:INPUT a* 

730 IF a$=“" THEN 760 
740 IF a* = "n"THEN 770 
750 GOTO 700 

760 PRINT:PRINT"ymln = "J:INPUT ymln :PRINT:PRINT"ym 
ax = “;:INPUT ymax 
770 RETURN 

780 REM UP Berechnung der Bildschirmpunkte fuer den 
Koordinatenursprung 

790 REM Erläuterungen fuer die Umrechnung auf Bildsc 
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hirm s. UP Funkt 1onsdarste11ung 

800 xb= 70+<500/(xmax-xmln))*<0-xmin) 

810 IF xmln>=0 THEN xb = 70 

820 IF xmax<=0 THEN xb = 570 

830 yb = 40+<320/<ymax-ymln))*(0-ymln> 

840 IF ymln >= 0 THEN yb = 40 
850 IF ymax <= 0 THEN yb = 360 
860 RETURN 

870 REM UP Zeichnen des Koordinatenkreuzes 

880 REM Erläuterungen fuer die Umrechnung auf Bildsc 

hlrm s. UP Funktionsdarstellung 

890 CLS 

900 ORIGIN 60,yb : DRAWR 520,0 
910 ORIGIN xb,30 : DRAWR 0,340 
920 RETURN 

930 REM UP Einteilung In dxg,dyg 

940 FOR k=0 TO dxg 

950 ORIGIN 70+k*<500/dxg>,yb-10 

960 DRAWR 0,20 

970 NEXT k 

980 FOR k=0 TO dyg 

990 ORIGIN xb-10,40+k*<320/dyg) 

1000 DRAWR 20,0 
1010 NEXT k 
1020 RETURN 

1030 REM UP Einteilung ln dxk,dyk 

1040 FOR k=0 TO dxg*dxk 

1050 ORIGIN 70+k*(500/<dxg*dxk)),yb-5 

1060 DRAWR 0,10 

1070 NEXT k 

1080 FOR k=0 TO dyg*dyk 
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1090 ORIGIN xb-5,40+k*< 320/<dyg*dyk>) 

1100 DRAWR 10,0 
1110 NEXT k 
1120 RETURN 

1130 REM UP Beschriftung der ACHSEN 
1140 TAG 

1150 REM Beschriftung y-Achse 
1160 IF xb = 570 THEN 1240 
1170 FOR k=0 TO dyg 

1180 ORIGIN xb-70,46+k*(320/dyg),xb-70,xb-15,400,0 
1190 f=ymln+k*((ymax-ymln)/dyg) 

1200 IF f=0 AND xb<>70 THEN 1220 
1210 ON darsty GOSUB 1460,1470,1480 
1220 NEXT k 
1230 GOTO 1290 
1240 FOR k=0 TO dyg 

1250 ORIGIN 585,46+k*<320/dyg),585,640,400,0 
1260 f=ymin+k*<(ymax-ymln)/dyg) 

1270 ON darsty GOSUB 1460,1470,1480 
1280 NEXT k 

1290 REM Beschriftung x-Achse 
1300 IF yb = 360 THEN 1380 
1310 FOR k=0 TO dxg 

1320 ORIGIN 27+k*<500/dxg),yb-15,27+k*(500/dxg),87+k 

*(500/dxg),yb-15,yb-40 

1330 f=xmln+k*<<xmax-xmln)/dxg) 

1340 IF f = 0 AND yb<>360 THEN 1360 
1350 ON darstx GOSUB 1460,1470,1480 
1360 NEXT k 
1370 GOTO 1430 
1380 FOR k=0 TO dxg 
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1390 ORIGIN 27+k*<500/dxg),390,27+k*<500/dxg),87+k*< 

500/dxg),395,375 

1400 f=xmln+k*<(xmax-xmln)/dxg) 

1410 ON darstx GOSUB 1460,1470,1480 
1420 NEXT k 
1430 TAGOFF 

1440 ORIGIN 0,0,0,640,400,0 
1450 RETURN 

1460 PRINT USING "#####.#“> f >:RETURN 

1470 PRINT USING "####.##“;f5:RETURN 

1480 PRINT USING "#######"»f5:RETURN 

1490 REM UP Berechnung und Ausgabe der Funktionswert 

e 

1500 REM Funktionsdarstellung lm "wlndow" 70,570,360 
,40 

1510 REM Achsenkreuze werden jeweils um 10 Punkte la 
enger ausgezogen 

1520 REM delta x = 500, delta y = 320 
1530 ORIGIN 70,40 

1540 a=0 : y=0 : IF ymin>=0 THEN y=ymin 
1550 FOR x=xmln TO xmax STEP (xmax-xmln)/500 
1560 a = a+1 

1570 GOSUB 60 : REM Berechnung des Funktionswertes 
1580 b = 320/<ymax-ymln)*<y-ymln) 

1590 PLOT a,b 
1600 NEXT x 
1610 RETURN 

1620 REM UP Texteintragung Ins Diagramm 

1630 IF yb=40 THEN WINDOW#l,1,80,1,1 : GOTO 1650 

1640 WIND0W#1,1,80,25,25 

1650 PRINT#1,“Sol 1 in das Diagramm Text eingetragen 
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werden < /n> ■5:INPUT#1,a* 

1660 IF a$ = " " THEN 1690 
1670 IF a*="n" THEN 1790 
1680 GOTO 1650 
1690 CLS #1 

1700 PRINT#1,"Text"»:INPUT#1,a* 

1710 PRINT#1,"Wohin soll der Text? (Cursor bewegen, 
dann ENTER-Taste druecken)" 

1720 LOCATE 1,1 

1730 INPUTJb* : IF b* = "" THEN PRINT at : GOTO 1750 
1740 GOTO 1710 

1750 LOCATE 1,1 : PRINT" ": CLS#1 : PRINT#1,"Noch el 
n Text < /n> ";:INPUT#1,a# 

1760 IF a# = " “ THEN 1690 
1770 IF a*="n* THEN 1790 
1780 GOTO 1750 
1790 CLS #1 
1800 RETURN 

1.2 Berechnung einer Wertetabelle 

Neben der Darstellung des Funktionsgraphen ist auch eine Wertetabelle 
ein nützliches Hilfsmittel, um einen Überblick über den prinzipiellen 
Verlauf einer Funktion zu bekommen. Außerdem können für besonders 
interessierende x-Werte die dazugehörigen Funktionswerte entnommen 
werden, so daß inter- und extrapoliert werden kann und die genauen 
Werte für Extremwerte, Wendestellen usw. festgehalten werden können. 

Das Programm Wertetab 1.0 benutzt wie fast alle anderen 
Programme in diesem Buch den Mode 2, d.h., es werden 80 
Zeichen pro Zeile auf dem Bildschirm dargestellt. 

Die Funktionsgleichung ist in 

Zeile 70 abgelegt; genau wie im Programm Funktion 1.0 wird sie als 
Unterprogramm in das übrige Programm eingebunden. In 
den 
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Zeilen 90 

bis 140 wird der Bildschirm in fünf Teilbereiche eingeteilt, wobei die 
ersten drei für die Ausgabe der Wertetabelle benutzt werden, 
der vierte für die Eingabe von x-Werten und der fünfte für die 
Ausgabe von Text zwecks Bedienerführung. 

Durch die Benutzung dieser WINDOWS bietet das Pro¬ 
gramm die Möglichkeit, bis zu drei formatierte Wertetabellen 
auf dem Bildschirm auszugeben, wobei bis zu 20 Werte in eine 
Tabelle eingetragen werden können. Daneben können aber 
auch die berechneten Funktionswerte ohne Formatierung 
auf dem Bildschirm auszugeben. Die Abfrage hierfür findet in 
den 

Zeilen 150 

bis 200 statt. 

Für die formatierte Ausgabe muß man über die 

Zeilen 210 

bis 290 angeben, wie die x-Werte in der Tabelle dargestellt werden 
sollen; die entsprechende Codezahl wird als darstx dem 
Programm eingegeben. In den 

Zeilen 300 

bis 380 erfolgt das Entsprechende für die v-Werte; die Codezahl ist 
hier darsty. In den 

Zeilen 420 

bis 480 wird der Kopf der Wertetabelle gezeichnet; über die 
FOR/NEXT-Schleife in den 

Zeilen 410 

und 590 bzw. über die Berechnung der x-Koordinaten in den 

Zeilen 420 

bis 450 wird der gleiche Programmabschnitt für alle drei Werte¬ 
tabellen benutzt, obwohl sie an verschiedenen Plätzen auf 
dem Bildschirm erscheinen. In 

Zeile 500 wird schließlich der x-Wert eingegeben, für den ein Funk¬ 
tionswert berechnet werden soll; nach Überprüfung der 
Endbedingung (s.u.) in 

Zeile 520 geschieht dies in 

Zeile 540 Die Ausgabe erfolgt in den 

Zeilen 550 
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bis 570 wobei entsprechend der Codezahlen darstx und darsty die 
Formatierung gewählt wird, die in den Unterprogrammen in 
den 

Zeilen 610 

bis 720 abgelegt ist. Für x wird hierbei wegen der Rechtsbündigkeit 
PRINT USING benutzt, während für y ROUND Verwen¬ 
dung findet und die Zahlen somit linksbündig ausgegeben 
werden. 

Sollen die Funktionswerte unformatiert ausgegeben werden, 
müssen die Programmschritte ab 
Zeile 700 benutzt werden. Hierbei wird in 
Zeile 770 ein gesondertes Ausgabefenster definiert; in 
Zeile 790 erfolgt die Eingabe des x-Wertes und in 
Zeile 800 die Berechnung des Funktionswertes. 

Erfolgt bei der Berechnung des Funktionswertes ein Fehler 
(z.B. Division durch 0), so wird über 

Zeile 20 

und 30 der Funktionswert auf 999999999 gesetzt und normal 
weitergemacht. 

Programmlänge: etwa 2.8 K. 

Anhand des Programms Wertetab 1.0 sei beispielhaft das 
Beenden der Eingabe einer Anzahl von Meßwerten erläutert. 
Für diese Aufgabe gibt es diverse Lösungsmöglichkeiten. Am 
einfachsten ist es, wenn die Anzahl der einzugebenden 
Meßwerte bekannt ist; in diesem Fall baut man entweder 
einen Zähler auf oder benutzt gleich eine FOR/NEXT- 
Schleife. Wenn die Anzahl der Meßwerte bei mehreren 
Programmabläufen nicht konstant ist, muß man vor der 
eigentlichen Meßwerteingabe die für den jeweiligen Durch¬ 
lauf gültige Zahl als Variable eingeben und diese Variable im 
Zähler benutzen. 

Eine weitere Möglichkeit zum Beenden der Eingabe besteht 
darin, daß man im Programm einen Abbruchwert definiert. 
Dies kann eine beliebige Zahl sein; allerdings sollte sie so 
gewählt werden, daß sie mit Sicherheit nicht als Meßwert 
vorkommt. Wenn man jetzt einen Meßwert eingegeben hat, 
prüft das Programm zunächst, ob der Eingabewert gleich 
dem Abbruchwert ist. Ist das der Fall, hört das Programm 
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mit der Meßwerteingabe auf und führt mit der Weiterver¬ 
arbeitung o.ä. fort. 

Das soeben beschriebene Verfahren wird z.B. im Programm 
A usgleich (siehe Kapitel 10) verwendet. Der Abbruchwert ist 
dort 333.33 (Zeile 230). Hier im Programm Wertetab 1.0 wird 
eine Variante des Verfahrens angewendet. Als Abbruchwert 
wird keine Zahl, sondern der Buchstabe n verwendet. Das hat 
den Vorteil, daß der Abbruch wert garantiert nicht als 
Eingabewert vorkommt; man muß allerdings in Kauf 
nehmen, daß die Meßwerteingabe etwas komplizierter wird. 
Da der Abbruchwert ein String ist, muß auch der jeweilige 
Meßwert zunächst als Stringvariable eingegeben werden, da 
sonst kein Vergleich stattfinden kann. Dieser wird in 
Zeile 520 durchgeführt; wurde n eingegeben, beendet das Programm 
die Eingabe von Werten. 

Ist die Abbruchbedingung nicht gegeben, so wird der String 
x$ zunächst in die Zahl x umgewandelt; dies passiert in 
Zeile 530 unter Benutzung der VAL-Funktion. Anschließend macht 
das Programm dann mit Berechnung des Funktionswertes, 
Ausgabe usw. weiter. 


10 REM Programm Wertetab 1.0 
20 ON ERROR GOTO 30 
30 y=999999999:RESUME NEXT 
40 MODE 2 
50 GOTO 90 

60 REM Zelle 60 = Funktionsglelchung 
70 y=SQR(x) 

80 RETURN 

90 REM Bildschirm wird ln 5 Bereiche eingeteilt 
100 WINDOW#1, 1,27,1,22 
110 WINDOW #2, 28,54,1,22 
120 WINDOW #3,55,80,1,22 
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130 WINDOW #4,1,80,23,24 : REM fuer die Eingabe der 
werte etc. 

140 WINDOW#5,1,80,25,25 : REM Fuer Ausgabe von Text 
150 CL3 

160 PRINT"Programm Wertetabe1le“ 

170 PRINT:PRINT:PRINT:PRINT“Sollen die Werte ln eine 
r formatierte Tabelle ausgegeben werden ( 

1) ":PRINT"oder sollen die Y 
-Werte jeweils einzeln ausgegeben werden 
(2)“:INPUT a 

180 IF a<1 OR a>2 THEN 150 
190 IF a=1 THEN 210 
200 IF a=2 THEN 760 
210 CLS 

220 PRINT"Wle soll der X-Wert ln der Tabelle darsges 
teilt werden ?“:PRINT:PRINT 

230 PRINT TAB(20>;“-mit einer Dezimalstelle " » TAB 
(60)» * <1) “ 

240 PRINT TAB(20);"-mit zwei Dezimalstellen “;TAB (6 
0 ) 5 "( 2 ) “ 

250 PRINT TAB(20)J“-mlt drei Dezimalstellen ";TAB (6 
0);“<3)* 

260 PRINT TAB<20);"-mit vier Dezimalstellen ";TAB <6 
o>;"<4>" 

270 PRINT TAB< 20);* - 1 m E-Format"5 TAB <60);*(5)" 

280 PRINT TAB(20);"-so wie er eingegeben wird“;TAB ( 
60>;"<6>" 

290 INPUT darstx 
300 CLS 

310 PRINT’Wle soll der y-Wert in der Tabelle darsges 
teilt werden ?":PRINT:PRINT 
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320 PRINT TAB(20),"-mit einer Dezimalstelle "5TAB 
(60);"<1)" 

330 PRINT TAB(20);"-mlt zwei Dezimalstellen ";TAB (6 
0 );"( 2 )" 

340 PRINT TAB(20)5“-m1t drei Dezimalstellen “STAB (6 
o> ;“< 3>- 

350 PRINT TAB(20)»"-m1t vier Dezimalstellen ";TAB (6 
0);*<4>" 

360 PRINT TAB(20)» "-1m E-Format";TAB (60);"(5)" 

370 PRINT TAB(20);"-so wie er berechnet wird*;TAB (6 
0);"<6>" 

380 INPUT darsty 
390 CLS 

400 PRINT#5,"Beenden der Eingabe durch Eingabe von 
n ." 

410 FOR k=l TO 3 

420 ORIGIN k*216-55,375 : DRAWR -130,0 
430 ORIGIN k*216-120,400 : DRAWR 0,-32 
440 ORIGIN k*216-121,400 : DRAWR 0,-32 
450 ORIGIN k*216 — 1 19,400 : DRAWR 0,-32 
460 LOCATE #k, 8,1 : PRINT#k,"X" 

470 LOCATE #k, 16,1 : PRINT#k,"Y" 

480 PRINT#k 

490 FOR n=3 TO 22 

500 INPUT#4,"X-Wert "Jxt 

510 CLS #4 

520 IF x*="n" THEN GOTO 730 
530 x=VAL<x*) 

540 GOSUB 70 : REM Funktlonsglelchung 

550 ON darstx GOSUB 610,620,630,640,650,660 : REM Au 

sdruck des x-Wertes 
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560 ON darsty GOSUB 670,680,690,700,710,720 : REM Au 
sdruck des y-Wertes 

570 DRAWR 0,-16 : DRAWR -1,0 : DRAWR 0,16 : DRAWR -1 

,0 : DRAWR 0,-16 : DRAWR 2,0 

580 NEXT n 

590 NEXT k 

600 GOTO 600 

610 PRINT#k,USING ■########.#“;x i :RETURN 
620 PRINT#k,USING "#######.##";x; :RETURN 
630 PRINT#k,USING -######.###";x i :RETURN 
640 PRINT#k,USING "#####.####";x i :RETURN 
650 PRINT#k,USING *#.##### A ~ A *;x i :RETURN 
660 PRINT#k,x;:RETURN 

670 PRINT#k,TAB<14);ROUND <y,l> : RETURN 
680 PRINT#k,TABC14)5ROUND <y,2> : RETURN 
690 PRINT#k,TAB<14);ROUND <y,3> : RETURN 
700 PRINT#k,TAB(14);ROUND (y,4) : RETURN 
710 PRINT#k,TAB<14>;USING”#.##### AÄAA ‘;y : RETURN 

720 PRINT#k,TAB(14);y : RETURN 
730 CLS#4 
740 CLS#5 
750 GOTO 750 

760 REM Teilprogramm Einzelausgabe 
770 CLS 

780 WIND0W#1,1,80,10, 15 

790 PRINT#5,"Beenden der Eingabe durch Eingabe von n 

M 

800 INPUT#4,"X-Wert : ";x$ 

810 IF x* = "n" THEN 870 
820 x=VAL(x$ > 
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830 GOSUB 70 : REM Funktlonsglelchgung 
840 PRINT#1#"X : ’txi TAB<20);"Y : ";y 
850 CLS#4 
860 GOTO 800 
870 END 

1.3 Berechnung von Funktionswerten nach dem 
Horner-Schema 

Grundlagen 

Will man den Funktionswert von ganzen rationalen Funktionen 
bestimmen, ist es oft vorteilhaft, ein Schema anzuwenden, das nach dem 
englischen Mathematiker William George FlORNER (1786 bis 1837) 
benannt worden ist. Das Horner-Schema hat seine überragende Bedeu¬ 
tung allerdings verloren, seit diese Aufgabe mit Hilfe der Computer stark 
vereinfacht worden ist; muß man jedoch «zu Fuß» rechnen, ist es oft 
unmöglich, mit einigermaßen vertretbaren Rechenaufwand zu einem 
Ergebnis zu kommen, wenn man das Schema nicht anwendet. Da die 
Polynomdivision außerdem durch das Horner-Schema erleichtert wird 
(siehe Abschnitt 1.4), wird es an dieser Stelle besprochen. 

Bei der Anwendung geht man folgendermaßen vor: Die Koeffizienten 
der ganzen rationalen Funktion 

y = a„x n + + ... + a, x + a 0 (1-3.1) 

werden nebeneinander geschrieben, wobei fehlende Glieder durch 0 zu 
berücksichtigen sind. In einer zweiten Zeile steht an erster Stelle eine 0. 
Die weiteren Spalten sowie die Glieder einer dritten Zeile werden wie folgt 
berechnet: 

□ die in der gleichen Spalte stehenden Zahlen der ersten und zweiten 
Zeile werden addiert und die Summe als c„, c„ _ j, ..., c x , c 0 in Zeile 
drei eingetragen. 

□ diese Summe wird jeweils mit a, also mit dem Wert, für den der 
Funktionswert berechnet werden soll, multipliziert; das Ergebnis 
schreibt man um eine Zeile versetzt in die zweite Zeile. 
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Diese beiden Regeln werden abwechselnd angewendet, bis folgendes 
Schema entstanden ist: 


a n a n -1 «„-2 a i a o 

0 ac„ ac„_ | ... a-c 2 ac t 

C n C n - 1 C n~ 2 ••• C l C 0 


(1.3.2) 


Nach den oben genannten Regeln ist c n = a n \ c n _ t = a n _ , + a • c n usw. 
Der Wert c 0 schließlich stellt den gesuchten Funktionswert an der Stelle a 
dar. 


Beispiel 

Der Funktionswert der Funktion 

y = x 5 - 2.V 4 + 4x 2 + x-5 (1.3.3) 

soll an der Stelle x = 2 bestimmt werden. Als Horner-Schema ergibt sich: 

1 -2 0 4 1 - 5 

0 2 0 0 8 18 

7 0 ~ 0 4 9 l7 

Der Funktionswert an der Stelle x = 2 ist also y = 13. 


Literaturhinweise 

G: [5], S. 318 bis 320; G: [17], S.43 bis 46; G: [24], S. 57 f.; G,W: [22], 
S. 210f.; G,W: [21], S. 85 bis 92. 

Programmbeschreibung 

Das Programm ist relativ kurz; es belegt nur 0,7 K im 
Speicher des CPC 464. In 

Zeile 20 wird der Grad des Polynoms eingegeben (entspricht dem 
Wert n). Danach erfragt das Programm in den 

Zeilen 40 

bis 120 im Dialog die einzelnen Faktoren a n bis a 0 , wobei man 
fehlende Glieder durch Eingabe von 0 berücksichtigen muß. 
In 
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Zeile 130 wird schließlich der x-Wert eingegeben, für den die Berech¬ 
nung durchgeführt werden soll (entspricht a). Das Horner- 
Schema läuft in den 


Zeilen 140 
bis 180 ab. 

Zeile 190 gibt den Funktionswert aus, und die 
Zeilen 200 

bis 230 fragen ab, ob eine weitere Berechnung durchgeführt werden 
soll. 


10 REM Programm Horner 

20 CLS : INPUT'Grad des Polynoms : “Jgrad 
30 DIM xges(grad+l,3) 

40 PRINT:PRINT:PRINT"Geben Sie die Polynomfaktoren e 
ln 

50 FOR k=grad T0 0 STEP -1 

60 LOCATE 1,8 : PRINT"a*;kJ" : ";:INPUT xges(k,l> 

70 LOCATE 6,8 : PRINT" 

80 NEXT k 

90 PRINT:INPUT "Eingabe l.o. < /n> ";a$ 

100 IF a$="" THEN 130 
110 IF a*=“n" THEN 40 
120 GOTO 90 

130 CLS : INPUT "x-Wert fuer die Funktionsberechnung 
: ";x 

140 xgesfgrad,2)=0 

150 FOR k=grad TO 0 STEP -1 

160 xges(k,3)= xges<k,l> + xges(k,2) : IF k=0 THEN 1 
80 

170 xges(k-1,2)=xges(k,3)*x 
180 NEXT k 

190 CLS : PRINT"X = “JxJ TAB<20);"Y = " 


>xges<0,3> 
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200 PRINT:PRINT:INPUT"Noch eine Berechnung ( /n> ";a 
$ 

210 IF a$ ="" THEN 130 
220 IF a*="n" THEN 240 
230 GOTO 200 
240 END 


1.4 Polynomdivision 

Jede ganze rationale Funktion f(x) läßt sich als Produkt von sogenannten 
Linearfaktoren (x — a ,) schreiben, wobei a, jeweils eine Konstante 
darstellt. Ist a t eine Nullstelle der Ausgangsfunktion, so entsteht nach dem 
Abspalten des Linearfaktors (.v — «,), d.h. nach Division von f(x) durch 
(x — tij), eine neue ganze rationale Funktion ./„(.*), deren Grad um 1 
niedriger ist als der von/(jc). Spaltet man von f a (x) erneut eine Nullstelle 
ab, so ist der Grad der neuen Funktion wieder um 1 niedriger usw. Dieses 
kann man so lange wiederholen, bis die ursprüngliche Funktion ganz als 
Produkt von Linearfaktoren geschrieben wird. So hat z. B. die Funktion 

J\x) = x i -2x 2 - 15 x (1.4.1) 

drei Nullstellen; sie liegen bei 

- v oi = 5; x 02 = 0; ,v 03 = — 3 
Somit kann man f(x) auch wie folgt schreiben; 

f(x) = (x - 5) * (x - 0) * (.* + 3) 

Durch Ausmultiplizieren kann man aus dieser Form wieder Gleichung 
(1.4.1) hersteilen. 

Die praktische Nutzanwendung der besprochenen Regel liegt z. B. in 
der Bestimmung von Nullstellen. Sind eine oder mehrere Nullstellen 
bekannt, so erhält man nach Abspalten der entsprechenden Linearfak¬ 
toren eine Funktion niedrigeren Grades, von der die restlichen Nullstellen 
leichter zu bestimmen sind. Sind z.B. von einer Funktion 4. Grades 2 
Nullstellen bekannt, so erhält man eine Funktion 2. Grades, auf die man 
die bekannte Lösungsformel (siehe auch Abschnitt 1.7) anwenden kann. 
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Das Abspalten von Linearfaktoren kann man sich mit Hilfe des in 
Abschnitt 1.3 beschriebenen Horner-Schemas stark erleichtern. Man 
setzt lediglich als Wert a die bekannte Nullstelle ein und erhält die 
Koeffizienten der neuen Funktion direkt in der dritten Zeile (siehe hierzu 
auch [24], S. 58 f). 

Hat die Ausgangsfunktion den Grad n, so ist der Faktor in der ersten 
Spalte der dritten Zeile der Koeffizient von x" ~ 1 , der in der zweiten Spalte 
der Koeffizient von x” ~ 2 usw. In der letzten Spalte muß 0 stehen, da 
ansonsten a keine Nullstelle darstellt. 

Bezugnehmend auf das als Gleichung (1.3.2) aufgestellte allgemeine 
Schema erhält man für die Funktion 

y = a n x n + a„_ i x n ~ 1 + ... + a { x + a 0 (1.4.2) 

demnach nach Abspalten des Linearfaktors (x — ä) die neue Funktion 

/«(*): 

y =/<,(*)*(* - a ) 

= ( c n x n ~ 1 + c„ _ , x” ~ 2 + ... + c 2 x + c,) *(x - a), 
wenn a eine Nullstelle der Funktion darstellt. 

Beispiel 

Ausgangsfunktion: f(x) = x 3 — 2x 2 — 15x 
Nullstelle: x 01 = 5 

Horner-Schema: 1 —2 —15 0 

0 5 15 0 

1 3 0 0 

Also ergibt sich nach Abspalten des Linearfaktors (x — 5) die neue Funk¬ 
tionsgleichung 

f a (x) = x 2 + 3x 

Somit kann man die Ausgangsfunktion auch wie folgt schreiben: 
f(x) = (x 2 + 3x)*(x - 5) 

Die Gleichung f a {x) kann aber leicht gelöst werden; man erhält die 
Nullstellen 

x oi = 0 
X 02 = 


3 
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Literaturhinweise 

G: [5], S. 318 bis 320; G: [17], S.43 bis 46; G: [24], S.57f.; G,W: [22], 

S. 210 f.; G,W: [21], S.85 bis 92 

Programmbeschreibung 

Die 

Zeilen 10 

bis 190 entsprechen den ersten Zeilen des in Abschnitt 1.3 beschrie¬ 
benen Programms Horner. Lediglich die nachgeschaltete 
Ausgabe ist anders. In den 

Zeilen 200 

bis 410 erfolgt die Ausgabe der Funktion und des Ergebnispolynoms 

mit dem eventuellem Rest (wenn hier ein von 0 verschiedener 
Wert auftaucht, war a keine Nullstelle). Über die 

Zeilen 420 

bis 440 kann vom Ergebnispolynom ein weiterer Linearfaktor 
abgespalten werden. In diesem Fall wird in den 

Zeilen 480 

bis 530 das Ergebnispolynom in den Speicherplatz für das Ausgangs¬ 
polynom umgespeichert. 

Schließlich kann über die 

Zeilen 450 

bis 470 noch eine neue Rechnung durchgeführt werden. 


Programmlänge: etwa 1,9 K. 



34 


Funktionen 


10 REM Programm poldlv 

20 REM Dieses Programm spaltet von einem Polynom unt 
er Anwendung des Horner-Schemas einen Linearfaktor a 
b und berechnet das Restpoly 

nom. Dies kann beliebig oft wiederholt werden. 

30 CLS : INPUT"Grad des Polynoms : "igrad 
40 DIM xges<grad+1,3) 

50 PRINT:PRINT:PRINT"Geben Sie die Polynomfaktoren e 
ln . " 

60 FOR k=grad TO 0 STEP -1 

70 LOCATE 1,8 : PRINT"a";k5" : ";:INPUT xges(k,l> 

80 LOCATE 6,8 : PRINT" 

90 NEXT k 

100 PRINT:INPUT "Eingabe i.o. ( /n> ";a* 

110 IF a$ = " " THEN 140 
120 IF a$=*n" THEN 50 
130 GOTO 100 

140 CLS : INPUT “Nullstelle : "Jx 

150 xges(grad,2)=0 

160 FOR k=grad TO 0 STEP -1 

170 xgesCk,3)= xges(k,l) + xges(k,2) : IF k=0 THEN 1 
90 

180 xges(k-l,2)=xges(k,3)*x 
190 NEXT k 

200 REM Ausgabe von Ursprungspolynom und Ergebnispol 
ynom 

210 CLS : PRINT"Ausgangspolynom : " 

220 FOR 1=grad TO 0 STEP -1 

230 IF xges<1,1) = 0 THEN 300 

240 IF xges( 1,1) > 0 THEN PRINT" + 

250 PRINT xges( i, 1 ) ; 
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260 IF 1=0 THEN 300 

270 IF i=1 THEN 290 

280 PRINT'x*"; 1»: GOTO 300 

290 PRINT'x“5 

300 NEXT 1 

310 PRINT:PRINT:PRINT'Berechnetes Polynom (Nullstell 
e = " ; x;") “ 

320 FOR 1 = grad TO 1 STEP -1 
330 IF xges(1,3) = 0 THEN 400 
340 IF xges (1,3) > 1 THEN PRINT' + '; 

350 PRINT xges<1,3)5 

360 IF 1-1=0 THEN 400 

370 IF 1-1=1 THEN 390 

380 PRINT'x*";1-1;" GOTO 400 

390 PRINT'x 

400 NEXT 1 

410 PRINT' Rest : ';xges<0,3);* / < x -';x;')“ 

420 PRINT:PRINT:PRINT:PRINT'So 11 das berechnete Poly 
nom noch durch eine andere Nullstelle dividiert werd 
en < /n> ?":PRINT:PRINT'Der 
Rest kann dabei nicht beruecksIchtlgt werden.':INPU 
T a* 

430 IF a< = ” THEN 480 

440 IF a*=*n' THEN 450 ELSE 420 

450 CLS : PRINT'Soll eine weitere Rechnung mit einem 
anderen Polynom durchgefuehrt werden ( /n)":INPUT a 

* 

460 IF a$ = " THEN RUN 10 

470 IF a*="n* THEN 540 ELSE 450 

480 REM Umspeicherung des Polynoms 



36 


Funktionen 


490 xges(grad,l) = 0 

500 FOR i=grad-1 TO 0 STEP -1 

510 xges(l,l> = xges(l+l,3) 

520 NEXT 1 

530 GOTO 140 

540 END 


1.5 Berechnung der Nullstellen einer Funktion mit 
Hilfe der Regula Falsi 

Grundlagen 

Nullstellenbestimmungen nach der Regula Falsi gehören zu den 
Näherungsverfahren, die dann angewendet werden, wenn eine Gleichung 
nicht oder nur mit erheblichem Aufwand geschlossen lösbar ist. Wie der 
Name sagt, erhält man mit einem solchen Verfahren nur einen 
Näherungswert für die tatsächliche Nullstelle. Normalerweise wird 
iterativ vorgegangen, so daß man durch Erhöhung der Rechenschrittzahl 
sehr hohe Genauigkeiten erreichen kann. Bei Verwendung eines 
Computers, bei dem ja die Anzahl der Rechenschritte praktisch keine 
Rolle spielt, erhält man oft den genauen Wert der Nullstelle. 

Bei Anwendung der Regula Falsi müssen zwei Näherungswerte Xj und 
x 2 bekannt sein; für Xi muß der Funktionswert f{x x ) kleiner 0 sein, für x 2 
muß f(x 2 ) größer 0 sein. Bei einer steigenden Kurve liegt x x also links von 
der Nullstelle und x 2 rechts; bei einer fallenden Kurve ist es umgekehrt. 

Die Werte von x, und x 2 kann man z. B. aus einer Wertetabelle (s. Ab- 
schn. 1.2) oder aus der Zeichnung des Graphen der Funktion (s. Abschn. 
1.1) entnehmen. Wenn die beiden Näherungswerte x t und x 2 bzw. die 
dazugehörigen Funktionspunkte P { (x, ,>■,) und P 2 (x 2 ,y 2 ) bekannt sind, 
werden P l und P 2 durch eine Gerade verbunden. Diese Gerade stellt eine 
Sehne der Funktion dar; daher wird die Regula Falsi auch oft als Sehnen- 
Verfahren bezeichnet. 
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Bild 1.1 Konvergierende Regula Falsi 


Die Gleichung der Sehne lautet: 


y-y i _ y 2 -yi 

x - X, x, - X, 


(1.5.1) 


x = beliebiger x-Wert 
y = dazugehöriger _y-Wert 

Gleichung (1.5.1) sagt nichts anderes aus, als daß die Steigung der Sehne 
konstant ist. 

Wenn x = x s ist (siehe Bild 1.1), gilt für die Sehne: y = 0. Somit erhält 
man für den Schnittpunkt der Sehne mit der x-Achse aus Gleichung 
(1.5.1): 


x s = x \ — 


y 2 - y t 


(1.5.2) 
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Wenn man jetzt x s = „vif setzt, mit diesem Wert aus der zu untersuchen¬ 
den Funktionsgleichung den dazugehörigen Funktionswert/(xf) berech¬ 
net, diesen als j, und xf als x, in Gleichung (1.5.2) einsetzt und somit eine 
neues x s berechnet, stellt man fest, daß x s näher an die Nullstelle 
heranwandert. Daher wiederholt man das Verfahren so lange, bis die 
erforderliche Genauigkeit erreicht ist. 

Wenn man P, und P 2 umgekehrt wählt, daß also/(x,) größer 0 und 
/(x 2 ) kleiner 0 ist, konvergiert das Verfahren nicht, wie in Bild 1.2 
dargestellt ist. Die erste Sehne ist die gleiche wie in Bild 1.1; setzt man 
dann jedoch s s = x, und verfährt wie oben angegeben, liegt der 
Schnittpunkt der neuen Sehne mit der x-Achse rechts von der Nullstelle; 
man ist also «übers Ziel hinausgeschossen». Wählt man P, und P 2 jedoch 
wie ursprünglich gesagt, kann man das Verfahren auf jede beliebige 
Funktion anwenden, egal, ob rechts- oder linksgekrümmt, ob steigend 
oder fallend. 



Bild 1.2 Nicht konvergierende Regula Falsi 
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Beispiel 

Für die Funktion 

y = x 2 — 9 

soll die Nullstelle gesucht werden, die zwischen 2,5 und 3,5 liegt, und zwar 
mit einer maximalen Abweichung von 0,1 bezogen auf den y-Wert. 


x = 2,5 

/(-v) = 

-2,75 

x =3,5 

f(x) = 

3,25 

Daraus folgt die Festlegung: 



II 

l_/i 

y i = 

-2,75 

»n 

cn 

II 

rs 

** 

y 2 = 

3,25 


Aus Gleichung (1.5.2) erhält man für den Schnittpunkt der Sehne mit der 
x-Achse 


- Y s 


= 2,5 + 2,75 


3, 5 - 2,5 
3,25 + 2,75 


= 2,958 


= x* 

Hierfür erhält man den Funktionswert 
/(x*) = 2,958 2 - 9 
= - 0,25 


Erneut eingesetzt in Gleichung (1.5.2) erhält man 


- v s 


= 2,958 + 0,25 


3,5 - 2,958 
3,25 + 0.25 


= 2.997 


Der Funktionswert hierfür ist — 0,018, so daß die geforderte Genauigkeit 
bereits in zwei Schritten erreicht ist. 


Literaturhinweise 


G: [4], S. 314; G,W: [5], S.684 bis 688. 
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Programmbeschreibung 


Zeile 50 
Zeilen 70 
bis 110 

Zeile 120 


Zeilen 130 
und 140 

Zeile 150 
Zeile 160 


Zeilen 200 
und 210 


Auch das Programm Regula Falsi ist sehr kurz; Es hat eine 
Länge von 0,9 K. In 

ist die zu untersuchende Funktionsgleichung abgelegt. In den 

erfolgt die Eingabe der Schätzwerte sowie der Genauigkeits¬ 
grenze. Die Größe n in 

stellt einen Zähler dar, mit der die Anzahl der Durchläufe 
festgestellt wird, die der Rechner benötigt, um die geforderte 
Genauigkeit zu erreichen. In den 

werden für beide Schätzwerte die entsprechenden Funktions¬ 
werte berechnet; die Regula-Falsi-Formel läuft in 
ab. In 

erfolgt die Prüfung, ob die gezielte Genauigkeit ausreichend 
ist; falls nicht, wird eine neue Rechnung durchgeführt, 
nachdem der Zähler n um 1 erhöht wurde. Reicht die 
Genauigkeit dagegen, wird die Nullstelle mit der Anzahl der 
benötigten Durchläufe sowie dem Funktionswert ausge¬ 
geben. In den 

wird schließlich abgefragt, ob noch eine weitere Nullstelle 
bestimmt werden soll. 
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10 REM Programm Regula Falsl 
20 MODE 2 
30 GOTO 70 

40 REM Zelle 50 = Funktlonsglelchung 
50 y = x*x - 9 
60 RETURN 
70 CLS 

80 INPUT ‘Schaetzwert 1 fuer xO <f(xO> < 0 !!!!) 
x 1 

90 PRINTrINPUT ‘Schaetzwert 2 fuer xO <f<xO> > 0 !! 
! ! ) "Jx2 

100 PRINT:PRINT:PRINT“Mlt welcher Genauigkeit soll d 
le Nullstelle Berechnet werden ? " 

110 PRINT:INPUT "Geben Sie die maximale Abweichung a 
n “Jgrenze 

120 n=l : REM Durchlaufzaehler 

130 x=xl : GOSUB 50 : yl=y 

140 x=x2 : GOSUB 50 : y2=y 

150 xs = xl - yl*<<x2-xl>/<y2-yl)> 

160 x = xs : GOSUB 50 : IF ABS<y> < ABS(grenze) THEN 1 
90 

170 xl=xs : yl=y 
180 n=n+1 : GOTO 150 

190 CLS: PRINT'x = ‘;xs;TAB<20);‘y = "Jy 

200 PRINT:PRINT”Zur Berechnung wurden "»n;“Durchlaeu 

fe benoetlgt" 

210 PRINT:PRINT"Noch eine Nul1ste1lenberechnung derg 
lelchen Funktion < /n) “>:INPUT a$ 

220 IF a* = ” THEN 70 
230 IF a* = "n* THEN END 
240 GOTO 210 
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1.6 Berechnung der Nullstellen einer beliebigen 
Funktion 

In diesem Abschnitt wird ein Programm beschrieben, das eine ganz 
typische Rechneranwendung darstellt, da eine Vielzahl von Rechnungen 
ausgeführt werden, die von Hand mit vertretbarem Zeitaufwand gar nicht 
durchgeführt werden könnten. 

Das Programm berechnet die Nullstelle einer gegebenen Funktion, 
indem die Kurve einfach mit einer sehr kleinen Schrittweite abgetastet 
wird. Wechselt das Vorzeichen des jeweiligen Funktionswertes, wird die 
Schrittweite halbiert und die Kurve in umgekehrter Richtung abgetastet. 

Dieses Verfahren wird so lange angewendet, bis die erforderliche 
Genauigkeit erreicht ist. Von der gefundenen Nullstelle geht es dann mit 
der ursprünglichen Schrittweite weiter. 

Gegenüber dem in Abschnitt 1.5 beschriebenen Regula-Falsi-Ver- 
fahren hat diese Methode vor allen Dingen den Vorteil, daß keine 
Schätzwerte für die Nullstelle benötigt werden. 

Die Funktionsgleichung ist in 
Zeile 70 abgelegt. Die 
Zeilen 90 

bis 140 dienen der Eingabe der Intervallgrenzen, der Schrittweite und 
der erforderlichen Genauigkeit. 

Das Abtasten der Funktion erfolgt in den 

Zeilen 150 

bis 230 wobei in 

Zeile 210 der beschriebene Vorzeichenwechsel stattfindet. In 
Zeile 240 erfolgt die Ausgabe der gefundenen Nullstelle, und über die 
Zeilen 250 

bis 260 wird die Suche fortgesetzt, nachdem eine Nullstelle gefunden 
worden ist. 

Länge des Programms: etwa 1 K. 

10 REM Programm NullRechner 

20 REM Die Nullstelle wird berechnet. Indem die Kurv 
e mit vorgegebener Schritt- weite abgetastet wird, b 
el Vorzelchenwechse1 die Sch 
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rittweite halbiert wird etc. 

30 ON ERROR GOTO 40 
40 RESUME NEXT 
50 GOTO 90 

60 REM Zelle 60 = Funktlonsglelchung 
70 y=-2*x*x*x*x*x + 10*x*x*x - 8*x 
80 RETURN 

90 CLS : PRINT“Geben Sie den Startwert der Nullstell 
ensuche ein.“ 

100 INPUT'Es wird ln positiver Richtung gesucht!“;xs 
110 PRINT : PRINT"Bls zu welchem x-Wert soll die Nul 
lstellensuche fortgesetzt werden “;:INPUT xe 
120 PRINT:INPUT'Mlt welcher Schrittweite soll vorgeg 
angen werden ";swO 

130 PRINT:PRINT"Mlt welcher Genauigkeit soll die Nul 
lstelle berechnet werden ?" 

140 INPUT"Geben Sie die maximale Abweichung ein. ";g 
renze 

150 x=xs : GOSUB 70 : yl=y 
160 k=2 : sw = swO 
170 x=x+sw 

180 IF x>= xe THEN 270 
190 GOSUB 70 

200 y2=y : IF ABS(y) < ABS(grenze) THEN 240 

210 IF SGN(y1)=SGN(y2) THEN 220 ELSE 230 

220 y1 = y2 : GOTO 170 

230 sw=-sw/k : yl=y2 : GOTO 170 

240 PRINT:PRINT'xO = “5x,"y = ";y 

250 x=x+swO : GOSUB 70 

260 sw = swO : yl=y : GOTO 170 

270 END 
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1.7 Lösen quadratischer Gleichungen 

Grundlagen 

Die allgemeine Form einer quadratischen Gleichung lautet 

a 2 x 2 + a x x + a 0 = 0 (1.7.1) 

Mittels Division durch a 2 führt man diese Gleichung in die sogenannte 
Normalform über: 


x 2 + ^-x + ^ = 0 (1.7.2) 

a 2 «2 

Durch Termersetzung erhält man 


x 2 + px + q = 0 bzw. 

x 2 + px = -q (1.7.3) 

Zu dieser Gleichung addiert man auf beiden Seiten die quadratische 
Ergänzung (p/2) 2 . 


x 2 + px + 




Somit kann man aber nach dem binomischen Lehrsatz die linke Seite 
folgendermaßen umformen: 



Jetzt wird auf beiden Seiten die Wurzel gezogen. 



Hieraus erhält man die beiden bekannten Lösungsformeln: 
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Zur Lösung einer quadratischen Gleichung stellt man also zunächst die 
Normalform her und setzt dann die Faktoren p und q in Gleichung (1.7.4) 
und (1.7.5) ein. 

Beispiel 

4x 2 + 16x — 48 = 0 
Normalform: x 2 + 4x — 12 = 0 

Hieraus folgen direkt die Lösungen: 

x, = -2 + Jl 2 + 12 = 2 

x 2 = -2- Jl 2 + 12 = -6 

Literaturhinweise 

G: [15], S. 145 bis 162; G: [5], S. 251 ff.; G: [6], S. 69 bis 72 
Allgemeine Literatur zu Kapitel 1 

G: [15], S. 107 ff.; G: [5], S.224 bis 225; G: [12], S. 176ff.; G,W: [22], 
S. 179 ff.; G,W: [6], S. 83 ff.; G,W: [24], S. 27ff. 

Program mbesch reibung 

Zeile 20 

bis 110 Eingabe der Faktoren a 2 , a , und a 0 
Zeilen 120 

und 130 Termersetzung 

Zeile 140 Für den Fall, daß der Wert unter der Wurzel kleiner 0 ist, 
erfolgt keine Berechnung, sondern der Sprung nach 

Zeile 250 
Zeile 150 

und 160 Berechnung der Lösungen 
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Zeile 170 Ausgabe 
Zeile 180 

bis 200 Abfrage, ob eine weitere Rechnung durchgeführt werden soll. 
Programmlänge: etwa 0,8 K. 


10 REM Programm QuadratgleIch 

20 CLS : PRINT'Geben Sie die Faktoren der quadratlsc 
hen Gleichung ein .* 

50 PRINT:PRINT*y = a2*x A 2 + al*x + aO“ 

60 LOCATE 1,6 : INPUT “a2 : “;a2 
70 INPUT “al : “ial 
80 INPUT’aO : "JaO 

90 PRINT:INPUT “Eingaben l.o. < /n> ";a* 

100 IF a* = ““ THEN 120 

110 IF a* = “n“ THEN 10 ELSE 90 

120 p=al/a2 

130 q=a0/a2 

140 IF <p*p/4 - q) < 0 THEN 250 
150 xOl = -p/2 + SQR<p*p/4 - q> 

160 x02 = -p/2 - SQR(p*p/4 - q> 

170 PRINT:PRINT:PRINT"xOl = “;x01 : PRINT “x02 = “;x 
02 

180 PRINT:PRINT"Noch eine Gleichung < /n> "5:INPUT a 
$ 


190 IF a*=*“ THEN 10 

200 IF a*="n“ THEN END ELSE 180 

250 PRINT:PRINT"Dle Gleichungen hat lmaglnaere Loesu 
ngen. Loesung leider nicht moeglich.“:GOTO 180 
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Differentialrechnung 


2.1 Bestimmung der Ableitung 
von rationalen Funktionen 

Grundlagen 

Eine ganze rationale Funktion hat die Form 

y = a { x"' + a 2 x" 2 + ... + a m .x nm (2.1.1) 

Die 1. Ableitung dieser Funktion berechnet sich wie folgt (s. auch [6], 
S. 198 ff): 

/ = a t n l x”‘~ l + a 2 n 2 x n2 ~ l + ... + a m n m x*"~ l (2.1.2) 

Die 2. Ableitung erhält man, indem man die gleiche Regel auf die 
Funktionsgleichung der 1. Ableitung anwendet. Der Grad von x 
verringert sich beim Bilden der Ableitung also um 1. Flieraus folgt, daß die 
Ableitung einer Geraden gleich einer Konstanten ist: 

y =a { x' 

y' = öj • 1 • x° = a, 

Die Ableitung einer Konstanten ist dagegen gleich 0, da die Ableitung 
einer Funktion ja gleichzusetzen ist mit der Steigung der Funktion; eine 
Konstante ist aber eine Parallele zur .v-Achse und hat somit keine 
Steigung. Das gleiche Ergebnis erhält man bei «sturer» Anwendung der 
Ableitungsregel: 

y = a 0 = a 0 ■ x° 
y' = 0 • a 0 ■ x~ I = 0 
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Auch auf gebrochene Exponenten von x läßt sich obige Regel an wenden: 
y = a, • y 1/2 = a, • l/x 


y' = l / 2 ‘ a \ ' x 


1 / 2 _ 


fli 



Beispiel 

y = 5y 4 — 2y 2 + 3 

/ = 20 .y 3 - 4y 
= 60 y 2 - 4 
/" = 120y 


Program mbesch reibung 


Die Funktion wird per Dialog in den 

Zeilen 30 

bis 140 eingeben. Auch bei Konstanten muß hierbei der Exponent 
von y eingegeben werden (er ist hier gleich 0). 

Die Berechnung der Ableitungen erfolgt in den 

Zeilen 180 

bis 250 Die restlichen Zeilen dienen lediglich der Ausgabe. Die 
Funktion selbst wird in den 


Zeilen 280 
bis 370 
Zeile 300 


Zeile 310 

Zeile 320 

Zeile 330 
Zeile 350 

Zeile 340 


ausgegeben. 

unterdrückt hierbei die Ausgabe des jeweiligen Gliedes, falls 
der Faktor vor y gleich 0 ist, da das Glied in diesem Fall nicht 
vorhanden ist. 

druckt ein + aus, wenn dies erforderlich ist. Ein — wird 
automatisch ausgegeben, 
gibt den Faktor aus. 

Falls der Exponent von y gleich 0 ist, erfolgt mittels 

sofort die Ausgabe des nächsten Gliedes. 

druckt dagegen ein zusätzliches x aus, wenn der Exponent 

gleich 1 ist. Trifft beides nicht zu, wird über 

bzw. 
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Zeile 360 ein x und der dazugehörige Exponent ausgegeben. Der 
beschriebene Vorgang wird in den 

Zeilen 380 

bis 670 für die 1. bis3. Ableitung wiederholt; und schließlich kann das 
Programm über die Abfrage in 
Zeile 680 neu gestartet werden. 

Programmlänge: etwa 2,3 K. 


10 REM Programm polablelt 

20 REM Mit diesem Programm koennen von ganzen ratlon 
alen Funktionen die ersten drei Ableitungen berechne 
t werden 

30 MODE 2 : PRINT’Untersucht wird eine ganze rationa 
le Funktion der Form” 

40 PRINT:PRINT"y = al*x~nl + a2*x Ä n2 + ... + am*x A nm 

M 

50 PRINT:INPUT "Wie viele Glieder hat die Funktion “ 
; m 

60 DIM ao(m),n<m) 

70 PRINT:PRINT"Geben Sie die einzelnen Faktoren ein. 

M 

80 FOR 1=1 TO m 

90 LOCATE 1,10 : PRINT ”a”;l;" : ”;:INPUT ;a0(l):PRI 
NT TAB (30); "n“;l5” : ";:INPUT n0(l) 

100 LOCATE 1,10: PRINT“ 

N 

110 NEXT 1 

120 PRINT:INPUT ‘Eingaben alle l.o. < /n> ”;a* 

130 IF a*="” THEN 150 

140 IF a*="n* THEN RUN 10 ELSE 120 

150 DIM al(m),nl(m) : REM 1.Ableitung 
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160 DIM a2(m),n2<m> : REM 2.Ableitung 
170 DIM a3(m),n3<m) : REM 3.Ableitung 
180 FOR 1=1 TO m 
190 a1(i) = a0(1) * nO(l) 

200 nl<1) = nO<1) - 1 
210 a2(1) = a1(i) * nl(I) 

220 n2(1) = nl(1) - 1 
230 a3<i) = a2<i> * n2(i) 

240 n3<1) = n2<1) - 1 
250 NEXT 1 

260 REM konstanter Faktor wird automatisch beruecksl 
chtlgt 

270 REM Ausgabe der Gleichungen mit den Unbekannten 

280 CL8 : PRINT’Glelchung H 

290 FOR 1=1 TO m 

300 IF aO(1) = 0 THEN 370 

310 IF aO(1) > 0 THEN PRINT" + 

320 PRINT aOC1); 

330 IF nO< I) = 0 THEN 370 
340 IF nO(1) <> 1 THEN 360 
350 PRINT’x ";:GOTO 370 
360 PRINT "x*";nO<1>;" 

370 NEXT i 

380 PRINT:PRINT:PRINT’1.Ableltung": 

390 FOR 1=1 TO m 

400 IF al<1) = 0 THEN 470 

410 IF a1( i ) > 0 THEN PRINT" + "; 

420 PRINT a 1 < 1); 

430 IF nl<i) = 0 THEN 470 
440 IF nl<I) <> 1 THEN 460 
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450 PRINT“x “;: GOTO 470 
460 PRINT“x A “;nl<l>;" “; 

470 NEXT 1 

480 PRINT:PRINT:PRINT“2.Ableitung" 

490 FOR 1=1 TO m 

500 IF a2<1) = 0 THEN 570 

510 IF a2(l> > 0 THEN PRINT“ + “J 

520 PRINT a2< 1 > > 

530 IF n2<l> = 0 THEN 570 
540 IF n2<1> <> 1 THEN 560 
550 PRINT'x “»:GOTO 570 
560 PRINT"x Ä “5n2< 1 >»" 

570 NEXT 1 

580 PRINT:PRINT:PRINT“3.Ableltung“ 

590 FOR 1=1 TO m 

600 IF a3(1) = 0 THEN 670 

610 IF a3(I) > 0 THEN PRINT“ + "J 

620 PRINT a3<1)5 

630 IF n3<1) = 0 THEN 670 

640 IF n3(i) <> 1 THEN 660 

650 PRINT"x *;:GOTO 670 

660 PRINT'x*“;n3<1);“ “J 

670 NEXT 1 

680 PRINT:PRINT:PRINT:PRINT"Sol 1 noch eine Gleichung 
untersucht werden ( /n> “;:INPUT a$ 

690 IF a* = "" THEN RUN 10 

700 IF a* = “n“ THEN 710 ELSE 680 

710 END 
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2.2 Nullstellenbestimmung nach Newton 

Grundlagen 

Genau wie die in Abschnitt 1.5 behandelte Regula Falsi handelt es sich 
beim Newton-Verfahren (Bild 2.1) um ein Näherungsverfahren zur 
Nullstellenbestimmung. Es hat den Vorteil, daß nur ein Schätzwert 
benötigt wird, der sowohl rechts als auch links von der Nullstelle liegen 
kann. Allerdings muß die Ableitung der Funktion bekannt sein. 

Der Funktionspunkt für den ersten Schätzwert x { ist P, (x, ,y l ). In 
diesem Punkt legt man eine Tangente an die Kurve und bestimmt den 
Schnittpunkt der Tangente mit der x-Achse. 

Die Tangente im Punkt P, hat die gleiche Steigung m wie die Funktion 
selbst in P ,: diese ist aber identisch mit der 1. Ableitung der Funktion 
f'(.x j) für den .v-Wert x t . Somit ist man in der Lage, die Gleichung der 
Tangente aufzustellen (vgl. Abschnitt 1.5): 

= /»=/'(*,) ( 2 . 2 . 1 ) 

x — x t 

x = beliebiger .x-Wert 
y = dazugehöriger r-Wert 

Für den Schnittpunkt der Tangente mit der x-Achse x = x s gilt 
wiederum: y — 0. 

Daraus folgt: 


,y s = x. 


Ti 

ff c,l 


( 2 . 2 . 2 ) 


Analog zu Abschnitt 1.5 wird hierfür der Funktionswert der zu 
untersuchenden Gleichung berechnet, die zu x s gehörige Ableitung der 
Funktion bestimmt und damit dann über Gleichung (2.2.2) der nächste 
Wert für x s berechnet. Dieses Verfahren wird so oft angewendet, bis die 
erforderliche Genauigkeit erreicht ist; die berechneten Werte sind immer 
wieder einzusetzen in die Formel 
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Beispiel 

Von der Funktion 

y = x 3 - Ix + 12 

ist die Nullstelle in der Nähe von .v, = — 3 zu bestimmen, und zwar mit 
einer Genauigkeit von 0,001 bezogen auf den _f-Wert. 

f'(x) = 3 x 2 — 1 
/(*,) = -3 + 7*3 + 12 = 6 
/'(.v,) = 3*3 2 - 7 = 20 
In Gleichung (2.2.3) eingesetzt, ergibt sich 


Hierfür berechnet sich als Funktionswert J (x 2 )= —0,837, so daß ein 
neuer Rechenschritt notwendig wird. 


f'(x 2 ) = 25,67 


*3 


-3,3 


0,837 

25,67 


= -3,267 


Für diesen Wert ist der Funktionswert f(x 3 ) = —6,35*10 4 so daß die 
Genauigkeit ausreicht. 


Literaturhinweise 

G: [6], S. 196; G: [4], S.312 bis 314; G,W: [21], S.74ff. 


Programmbeschreibung 

Zur Nullstellenbestimmung nach Newton müssen sowohl die 
Funktionsgleichung als auch die 1. Ableitung dem Pro¬ 
gramm bekannt sein; beide sind in 
Zeile 60 bzw. 

Zeile 90 abgelegt. Durch 
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Zeile 110 

und 120 werden Schätzwert und Genauigkeitsanforderung dem Pro¬ 
gramm mitgeteilt; die eigentliche Rechnung wird in den 

Zeilen 140 

bis 200 durchgeführt. 

Zeile 220 dient der Ausgabe, und über 
Zeile 230 

bis 250 kann eine weitere Nullstelle gesucht werden. 

10 REM Programm Nullnewton 

20 REM Programm zur Nullstellenbestlmmung mit dem Ve 
rfahren nach Newton 
30 GOTO 110 

40 REM In der Zelle 60 muss die Funktlonsglelchung, 
ln der Zelle 90 die 1.Ableitung abgespeichert sein 

50 REM Funktlonsglelchung 

60 y = x*x*x - 4*x*x - ll*x + 30 

70 RETURN 

80 REM 1.Able1tung 
90 y = 3*x*x - 8*x - 11 
100 RETURN 

110 CL3 : PRINT’Mlt welcher Genauigkeit sollten) die 
Nullstelle(n) berechnet werden ? Geben Si 

e die maximale Abweichung el 
n .";:INPUT abw 

120 CLS : PRINT“Geben Sie einen Schaetzwert fuer die 
Nullstelle ein ."»:INPUT x 

130 REM so, jetzt geht’s los (oder wieder los) 

140 GOSUB 60 : REM Funktionswert fuer x 
150 IF ABS(y) <= ABS(abw) THEN 210 
160 yl = y 
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170 GOSUB 90 : REM 1. Ableitung fuer x 
180 yi 1 = y 

190 x = x - yi/yil : REM neuer Naeherungswert 

200 GOTO 130 

210 REM wow, geschafft! 

220 PRINT“xO = ";x,"y = “;y 

230 PRINT:PRINT“Noch eine Nullstelle bei derselben G 
leichung ( /n> :INPUT a$ 

240 IF a$ = ““ THEN 120 

250 IF a$ = “n“ THEN END ELSE 230 

2.3 Extremwerte von Funktionen 

Grundlagen 

Die Extremwerte einer Funktion (Minimum bzw. Maximum) sind 
dadurch gekennzeichnet, daß die Funktion in diesem Punkt waagerecht 
verluft (Bild 2.2 und 2.3). 

Legt man im Minimum bzw. Maximum eine Tangente an die Kurve, so 
ist diese eine Parallele zur x-Achse; die Steigung der Tangente ist gleich 0. 
Da Steigung von Tangente und Funktion per Definition aber identisch 
sind, gilt, daß auch die Steigung der Funktion gleich 0 ist. 

Dies heißt aber mit anderen Worten, daß die 1. Ableitung in diesem 
Punkt gleich 0 ist (s. auch [6], S. 207 ff). 

Um den Extremwert einer Funktion zu bestimmen, muß man also 
lediglich eine Nullstelle der 1. Ableitung suchen. Dies reicht jedoch nicht 
aus, um darüber zu entscheiden, ob ein Minimum oder ein Maximum der 
Funktion vorliegt. Um hierüber eine Aussage treffen zu können, muß 
man sich den Verlauf der 2. Ableitung vergegenwärtigen, d. h. den Verlauf 
der Steigung der 1. Ableitung (Bild 2.4). 

Im Bild 2.4 sind mit steigendem x-Wert Tangenten an eine Kurve mit 
Maximum gelegt. Wie man sieht, nimmt die Steigung der Tangenten, also 
auch die der Funktion, mit steigendem x-Wert bei Annäherung an das 
Maximum ab; die Tangenten werden immer flacher, bis im Maximum 
eine waagerechte Tangente vorliegt. Mit anderen Worten heißt das, daß 






Bild 2.2 Funktion mit Minimum 



Bild 2.3 Funktion mit Maximum 
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Bild 2.4 Tangenten verlauf bei einer Kurve mit Maximum 


die 1. Ableitung der Funktion immer kleiner wird, bis sie im Maximum zu 
0 wird, um dann in den negativen Bereich zu gehen (die Tangenten kippen 
zur anderen Seite). Wenn die 1. Ableitung jedoch fallt (kleiner wird), heißt 
das, daß ihre Steigung negativ ist oder daß die 2. Ableitung kleiner 0 ist. 

Zusammenfassend kann man sagen, daß folgende Bedingungen erfüllt 
sein müssen, wenn ein Maximum vorliegt: 

Maximum: 1. Ableitung = 0 
2. Ableitung < 0 

In gleicher Weise kann man die Bedingungen für ein Minimum beweisen: 

Minimum: 1. Ableitung = 0 
2. Ableitung > 0 
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Bild 2.5 Funktion mit Rechts- und Linkskrümmung 


Um den Krümmungsverlauf einer Kurve zu untersuchen, muß man 
ähnliche Betrachtungen anstellen: 

Die Funktion in Bild 2.5 beschreibt zunächst eine Rechtskurve, die 
Steigung der Tangente bzw. Kurve oder der Wert der 1. Ableitung nimmt 
ab. Von x w an aufwärts jedoch beschreibt die Funktion eine Linkskurve, 
d.h., ihr Drehsinn ist dann nach links gerichtet. lnx w geht die 1. Ableitung 
also von einer fallenden in eine steigende Bewegung über; die Steigung der 
I . Ableitung (gleichbedeutend mit der 2. Ableitung) gehl demnach aus 
dem negativen Bereich (= fallende Bewegung) in den positiven Bereich 
(= steigende Bewegung) über, in x w geht sie durch 0. Für den Wechsel von 
einer Rechtskrümmung in eine Linkskrümmung (= Wendestelle) gilt 
also, daß die 2. Ableitung eine Nullstelle haben muß. Weiterhin wurde 
festgestellt, daß die 2. Ableitung aus dem negativen in den positiven 
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Bereich übergeht, sie ist also steigend. Dies ist aber wiederum 
gleichbedeutend mit der Aussage, daß die Steigung der 2. Ableitung 
(= 3. Ableitung) positiv ist. 

Zusammenfassend lauten die Bedingungen für eine Rechts/Links- 
Wendestelle: 


R/L-Wendestelle: 2. Ableitung = 0 
3. Ableitung > 0 

Einen analogen Beweis kann man für eine Links/Rechts-Wendestelle 
führen: 

L/R-Wendestelle: 2. Ableitung = 0 
3. Ableitung < 0 

Beispiel 

Von der Funktion 

y = * 3 + 3x 2 - 9.x + 5 

sollen Extremwerte und Wendestellen bestimmt werden. Um diese 
Aufgabe lösen zu können, werden zunächst die ersten drei Ableitungen 
der Funktion ermittelt (vgl. Abschnitt 2.1) 

y' = 3.x 2 + 6x -9 
y" = 6.x + 6 
/" = 6 


Extremwerte 

Bedingung für die Extremwerte war: y' = 0. Also werden die Nullstellen 
der 1. Ableitung gesucht (vgl. Abschnitt 1.7): 

0 = 3.x 2 + 6.x - 9 

Normalform: 0 = jc 2 4- 2x — 3 
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Bei x = 1 und x = — 3 liegen also Extremstellen vor. Um zu prüfen, ob ein 
Minimum oder ein Maximum vorhanden ist, werden für diese x-Werte die 
jeweiligen Werte der 2. Ableitung berechnet. 

x= 1 y"= 12 

x = — 3 y" = —12 

fürx = 1 ist / = 0 und y" > 0; somit liegt ein Minimum vor. Fürx = — 3 
dagegen ist y' = 0 und y" < 0, womit nachgewiesen ist, daß es sich an 
dieser Stelle um ein Minimum der Funktion handelt. 

Wendestellen 

Wenn eine Wendestelle vorliegt, muß die 2. Ableitung zu 0 werden. 

0 = 6x + 6 
x — — 1 

Bei x = — 1 hat die Funktion also eine Wendestelle. Für die Entscheidung 
Links/Rechts- oder Rechts/Links-Wendestelle wird der Wert der 3. Ablei¬ 
tung für x = — 1 untersucht. 

x = -1 y'" = 6 

Da y'" größer als 0 ist, handelt es sich um eine Rechts/Links-Wendestelle. 

In Bild 2.6 sind die Funktion und alle drei Ableitungen dargestellt, um 
die in diesem Abschnitt angesprochenen Zusammenhänge zu verdeut¬ 
lichen. 


Allgemeine Literaturhinweise zu Kapitel 2 

G, W: [4]; G, W: [6], S. 166ff.; G, W: [22], S. 252 ff.; G, W: [24], S. 98 ff.; 
W: [20] 
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Bild 2.6 Darstellung der Funktion y — x 3 + ix 2 — 9x + 5 mit ihren 
Ableitungen 
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Zeilen 

bis 

Zeilen 

bis 

Zeilen 

bis 

Zeilen 

bis 

Zeilen 

bis 


Zeilen 

bis 

Zeile 

Zeilen 

bis 

Zeilen 

bis 


Program mhesch reibung 

Mit Hilfe des Programms Extremwerte können alle Null¬ 
stellen, alle Minima/Maxima und alle Wendestellen einer 
Funktion bestimmt werden; hierzu werden die Nullstellen 
der Funktion bzw. der entsprechenden Ableitungen mit 
Hilfe des in Abschnitt 1.6 beschriebenen Programms Null¬ 
rechner ermittelt. Demzufolge müssen die Funktionsglei¬ 
chung und die ersten drei Ableitungen dem Programm 
vorliegen; sie sind in den 
40 

150 abgelegt. Die 

160 

270 stellen das Hauptprogramm dar, wobei in den 

170 

210 die einzelnen zur Berechnung notwendigen Unterprogram¬ 
me aufgerufen werden. In den 

220 

270 ist die Abfrage abgelegt, ob eine weitere Funktion unter¬ 
sucht werden soll. In den 

670 

870 wird vom Rechner auf die Notwendigkeit hingewiesen, daß 
alle drei Ableitungen dem Programm vorliegen müssen. 
Falls das noch nicht der Fall ist, hat man hier die 
Möglichkeit, sie noch einzugeben. Das Programm stoppt in 
diesem Fall; nach Eingabe der Gleichungen wird es mit 
GOTO 860 wieder gestartet. 

Die Eingabe von Untersuchungsgrenzen und Genauigkeit 
der Nullstellenbestimmungen erfolgt in den 

280 

350 Die Anzahl der Rechenschritte wird in 

320 hierbei auf 100 festgelegt (vgl. Abschnitt 1.6). Die 

360 

440 stellen die Benutzerführung vor der Nullstellenbestimmung 
dar; die eigenliche Rechnung findet in den 

880 

1130 statt. Dieser Programmabschnitt entspricht dem Programm 
Nullrechner; lediglich die Berechnung des Funktionswertes 
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Zeile 

Zeile 

Zeile 

Zeile 

Zeile 


Zeile 


Zeile 


Zeile 

Zeile 


Zeile 

Zeile 


mußte etwas komplizierter gestaltet werden, da das gleiche 
Unterprogramm auch zur Bestimmung der Nullstellen der 
Ableitungen benutzt wird; somit sind die Sprungadresse für 
das GOSUB aber variabel. 

Hier wurde das Problem so gelöst, daß zunächst die 
Variable wo in 

410 auf 1 gesetzt wird und hiermit dann in 
940 sowie in 
1000 

1070 und 

1090 die richtige Sprungadresse für die Funktionsgleichung 
ausgesucht wird. 

Wenn die Rechnung zu den Extremwerten kommt, wird wo 
in 

510 auf 2 gesetzt und dann in den gleichen Zeilen die Gleichung 
der 1. Ableitung zur Nullstellenbestimmung benutzt. 

Das Entsprechende geschieht für die Wendestellen; hier 
wird wo in 

610 auf 3 gesetzt und damit dann für die 3. Ableitung die 
Nullstellen berechnet. 

Ist eine Nullstelle gefunden, wird sie ausgegeben und die 
Rechnung fortgesetzt. Wenn die Intervallgrenze erreicht ist, 
erfolgt der Sprung zurück ins Hauptprogramm, wo dann 
mit dem nächsten Programmabschnitt weitergemacht wird. 
Wenn die Extremstellen berechnet werden, wird neben den 
x-Werten auch noch der Wert der 1. Ableitung über 
1150 für diesen x-Wert ausgegeben; dieser muß gleich 0 sein, 
wenn tatsächlich eine Nullstelle vorliegt. In 
1160 wird außerdem noch der Wert der 2. Ableitung ausgegeben, 
damit man entscheiden kann, ob ein Maximum oder ein 
Minimum vorliegt. Die Berechnung der Wendestellen 
erfolgt analog zu der Bestimmung der Extremwerte. 

Bei der Ausgabe der jeweils berechneten Werte wartet das 
Programm so lange, bis eine beliebige Taste gedrückt wird. 
Dies wird erreicht durch die INKEY-Funktion in 
1120 Vorher wird in 

1110 der entsprechende Text ausgegeben. 

Genau der gleiche Effekt wird erreicht durch Aufrufen einer 
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entsprechenden Maschinenspracheroutine; dies geschieht 
mit CALL &BB06. 

Dieser Befehl ist zwar etwas kürzer zu tippen; die Adresse 
stimmt jedoch nur für den CPC464. Um eine möglichst 
einfache Übertragung der Programme zu gewährleisten, 
wird in diesem Buch nur mit der INKEY-Funktion 
gearbeitet. 

Programmlänge: etwa 4,6K. 


10 REM Programm extremwerte 

20 REM Mit diesem Programm werden die Extremwerte vo 
n beliebigen Funktionen berechnet (einschliesslich N 
ullstellen) 

30 GOTO 160 

40 REM FunktIonsgleichung 
50 y = x*x*x + 3*x*x - 9*x + 5 
60 RETURN 

70 REM 1.Ableitung 
80 y= 3*x*x + 6*x - 9 
90 RETURN 

100 REM 2.Ableitung 
110 y= 6*x + 6 
120 RETURN 

130 REM 3.Ableitung 
140 y=6 
150 RETURN 

160 REM hier geht’s los 

170 GOTO 670 : REM Eingabe der Funktion und der Able 
itungen 

180 GOSUB 280 : REM Untersuchungsgrenzen und Genauig 
ke 11 
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190 GOSUB 360 : REM Nullstellen 
200 GOSUB 470 : REM Extremstellen 
210 GOSUB 570 : REM Wendestellen 

220 CLS : PRINT"So1 len von einer weiteren Funktion d 
ie Extremwerte bestimmt werden ( /n) ";:INPUT a$ 

230 IF a$ = “ “ THEN RUN 
240 IF a$<>-n- THEN 220 

250 CLS : LOCATE 20,10 : PRINT'Na, denn man Tschuess 

260 LOCATE 1,24 
270 END 

280 REM UP Eingabe von Untersuchungsgrenzen und Gena 
uigke it 

290 CLS : PRINT’Geben Sie das Intervall fuer die Un 
tersuchungen ein.“ 

300 PRINT:PRINT TAB(20);“xmin ";:INPUT xs 

310 PRINT:PRINT TAB(20);"xmax “;:INPUT xe 

320 swO = (xe-xs)/100 : REM Schrittweite 

330 PRINT:PRINT"Mit welcher Genauigkeit sollen die 

Extremwerte berechnet werden ?“ 

340 PRINT:PRINT"Geben Sie die maximale Abweichung ei 
n. “J:INPUT abw 
350 RETURN 

360 REM UP Nul1ste1lenbestimmung 

370 CLS : INPUT“Sollen die Nullstellen der Funktion 
bestimmt werden ( /n)“;a$ 

380 IF a$ = THEN 400 

390 IF a$ = “n“ THEN 460 ELSE 360 

400 b$ = “Nullstelle" : c$ = “y =“ 

410 wo = 1 

420 CLS : PRINT“Nu 11ste11en der Funktion“ 
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430 PRINT 


440 PRINT 

450 GOSUB 880 : REM Nu 11ste1lenberechnung 
460 RETURN 

470 REM UP Extremwertbestimmung 

480 CLS : PRINT"Sollen die Extremwerte der Funktion 
bestimmt werden ( /n) “;:INPUT a$ 

490 IF a$="" THEN 510 

500 IF a$=’n“ THEN 560 ELSE 470 

510 b$=’Extremwert" : c$ = ’y' =" : wo = 2 

520 CLS : PRINT"Extremwerte der Funktion’ 

530 PRINT" _ 


540 PRINT 

550 GOSUB 880 : REM Extremwertberechnung 
560 RETURN 

570 REM UP Wendeste 1lenbestimmung 

580 CLS : PRINT’Sollen die Wendestellen der Funktion 
bestimmt werden ( /n) ";:INPUT a$ 

590 IF a$ = ” THEN 610 

600 IF a«="n" THEN 660 ELSE 570 

610 b$ = "Wendestel len“ : c$ = “y" =“ : wo = 3 

620 CLS : PRINT’Wendeste11en der Funktion’ 

630 PRINT’ _ 


640 PRINT 

650 GOSUB 880 : REM Wendeste11enbestimmung 
660 RETURN 

670 REM Funktionseingabe fuer beliebige Funktionen 
680 CLS 
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690 PRINT"Funktionsgleichung und je nach Anwendungsf 
all 1.-3.Ableitung muessen dem Programm vorli 

egen, und zwar die” 

700 PRINTrPRINT TABC20);"- Funktionsgleichung in Zei 
le 50" 

710 PRINT:PRINT TAB(20);”- 1. Ableitung in Zeile 80” 

720 PRINTrPRINT TAB(20);"- 2. Ableitung in Zeile 110 

730 PRINTrPRINT TAB<20);"- 3. Ableitung in Zeile 140 

740 PRINT” 


750 PRINTrPRINT TAB(20);”Das ist der Fall";TAB(60);" 
CD" 

760 PRINTrPRINT TAB(20);"Das muss noch gemacht werde 
n”;TAB(60);“(2)" 

770 INPUT a 

780 IF a<1 OR a>2 THEN 670 

790 IF a=1 THEN 870 

800 LOCATE 1,14 r PRINT " 

810 LOCATE 1,16 r PRINT” 

820 PRINT“ 

830 LOCATE 1,14 

840 PRINTrPRINT'Geben Sie die Gleichungen jeweils in 
der Form y=.... in den oben angegeben Zeilen e 

in. Starten Sie das Programm 
anschliessend wieder mit goto 860” 

850 STOP 

860 REM weiter geht’s 





Differentialrechnung 


69 


870 GOTO 180 

880 REM Programm NullRechner 

890 REM Die Nullstelle wird berechnet, indem die Kur 
ve mit vorgegebener Schritt- weite abgetastet wird, 
bei Vorzeichenwechsel die Sc 
hrittweite halbiert wird etc. 

900 ON ERROR GOTO 910 
910 RESUME NEXT 
920 Zaehler = 0 
930 x = xs 

940 ON wo GOSUB 50,80,110 
950 y 1 =y 

960 k=2 : sw = swO 

970 Zaehler = Zaehler + 1 

980 x=x+sw 

990 IF x>= xe THEN 1110 
1000 ON wo GOSUB 50,80,110 

1010 y2=y : IF ABS(y) < ABS(abw) THEN 1050 
1020 IF SGN(yl)=SGN<y2) THEN 1030 ELSE 1040 

1030 y1=y2 : GOTO 980 

1040 sw=-sw/k : y1=y2 : GOTO 980 

1050 PRINT:PRINT b$;Zaehler;" = ";x, 

1060 GOSUB 50 : PRINT'y = ";y 
1070 ON wo GOSUB 1130,1140,1180 
1080 x=x+sw0 

1090 ON wo GOSUB 50,80,110 

1100 sw = swO : yl=y : GOTO 970 

1110 PRINT:PRINT'Weiter mit beliebiger Taste" 

1120 IF INKEYS="" THEN 1120 
1130 RETURN 
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1140 REM UP Ausgabeprogramm fuer Extremwertsuche 
1150 GOSUB 80 : PRINT TAB(12);'y' = "ly, 

1160 GOSUB 110 : PRINT'y” = “Jy 
1170 RETURN 

1180 REM UP Ausgabeproramm fuer Wendestellensuche 
1190 GOSUB 110 : PRINT TAB<13);-y’' = 'ly. 

1200 GOSUB 140 : PRINT’y"’ = "ly 
1210 RETURN 
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3.1 Integration nach der Trapezregel 

Grundlagen 

Die Bestimmung eines Flächeninhaltes ist nur lösbar, wenn die Fläche 
durch Geraden begrenzt wird bzw. wenn geometrisch genau definierte 
Formen wie Kreise, Elipsen usw. vorliegen. Schwieriger wird es bei der in 
Bild 3.1 dargestellten Figur: 

Es handelt sich hier um die in der Technik sehr häufig vorkommende 
Aufgabe, die Fläche zwischen Kurve und .Y-Achse innerhalb der Grenzen 
x A und x E zu bestimmen. Diese Fläche wird auch als bestimmtes Integral 
bezeichnet. Ist die Funktion der Kurve bekannt, wird man normalerweise 
in der Lage sein, die Kurve zu integrieren und somit die Fläche zu 
berechnen (vgl. z. B. [6], S. 261 ff). 

Es besteht jedoch auch die Möglichkeit, numerisch zu integrieren, d.h. 
Näherungsverfahren anzuwenden. Solche Verfahren muß man auch dann 
benutzen, wenn die Funktion nicht integrierbar ist. 

Will man die Integration nach der sogenannten Trapezregel vorneh¬ 
men, teilt man die zu berechnende Fläche in n gleich breite Streifen. 
Die Breite eines solchen Streifens ist 

h = Xa ~ X - e (3.1.1) 

n 

Verbindet man die Funktionswerte an den einzelnen Stellen .y ( durch 
Geraden, so berechnet sich die Fläche eines einzelnen Streifens zu 


_ Ti+ 1 + Ti . h 


(3.1.2) 




Bild 3.2 Numerische Integration nach der Trapezregel 
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Die Summe aller einzelnen Streifen ergibt aber die gesamte Fläche: 






4 + y t + y 2 + ... +y 


— + 2 ! ) 


(3.1.3) 


Gleichung (3.1.3) stellt die gesuchte Trapezregel dar. Je höher man die 
Anzahl n der Streifen wählt, desto genauer wird die Rechnung, da dann 
durch die schmäler werdenden Streifen die Abweichungen vom tatsäch¬ 
lichen Verlauf der Funktion abnehmen. 


Beispiel 

Die Fläche zwischen der Funktion 

y = 21 x 

und der x-Achse soll in den Grenzen x A = 2 und x E = 4 bestimmt werden. 
Als Streifenanzahl wird n = 10 gewählt. 

Zur Lösung stellt man folgende Tabelle auf: 


1 

-V 

Ti =/(*,) = 2/*,- 

0 

2,0 

1 

1 

2,2 

0,909090909 

2 

2,4 

0,833333333 

3 

2,6 

0,769230769 

4 

2,8 

0,714285714 

5 

3,0 

0,666666667 

6 

3,2 

0,625 

7 

3,4 

0,588235294 

8 

3,6 

0,555555556 

9 

3,8 

0,526315789 

10 

4,0 

0,5 
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Als Streifenbreite h ergibt sich 

h = "Tr = 0,2 

Somit berechnet sich die Fläche zu 

A = h(y 0 /2 + £ V,- +>’io/2) 

i = 1 

= 1,38754281 

Die Fläche unter der Kurve beträgt also 1,39 Flächeneinheiten. 

Literaturhinweise 

G: [4], S. 304; G: [5], S. 356 bis 358; G: [22], S. 399 bis 400 


Programmbeschreihung 


Die zu untersuchende Funktionsgleichung ist wieder in 
Zeile 50 abgelegt. Über 

Zeile 90 bzw. das Unterprogramm in den 

Zeilen 1000 bis 

bis 1050 erfolgt die Eingabe des Intervalles, in dem die Funktion 
integriert werden soll sowie die Vorgabe der Schrittweite. 
Die Berechnung des Integrals wird in 
Zeile 110 bzw. in 

Zeile 2100 

bis 2220 vorgenommen, wobei die Sehnen-Trapez-Formel zum 

Tragen kommt. 

Das Ergebnis wird als Variable integ gespeichert. In 
Zeile 130 erfolgt über das Unterprogramm in den 
Zeilen 1300 

bis 1390 die Ausgabe des Integrals sowie der Intervallgrenzen. 


Länge des Programms: etwa 1,3 K. 
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10 REM Programm Itrapez 

20 REM Integration einer vorgegebenen Funktion mit d 
er Sehnen-Trapez-Forme1 
30 GOTO 70 

40 REM Zelle 50 = Funktlonsglelchung 

50 y=l/x 

60 RETURN 

70 REM los geht's 

80 MODE 2 

90 GOSUB 1000 : REM Startwerteingabe etc. 

100 CLS 

110 LOCATE 20,10 : PRINT"Ich arbeite, Bitte warten." 
120 GOSUB 2100 : REM Integration nach der Sehnen-Tra 
pez-Formel 

130 GOSUB 1300 : REM Ausgabe 
140 END 

1000 REM UP Startwerteingabe etc. 

1010 CLS 

1020 PRINT“In welchen Grenzen soll die Funktion Inte 
grlert werden ?" 

1030 PRINT:PRINT"von xmln = "J:INPUT xmin:PRINT"bis 
xmax = ";:INPUT xmax 

1040 PRINT:PRINT"Mlt welcher Schrittweite soll der R 

echner Vorgehen ";:INPUT h 

1050 RETURN 

1300 REM UP ausgabe 

1310 CLS 

1320 PRINT a* : PRINT"_" : P 

RINT 

1330 PRINT"Integral : "Jlnteg 
1340 PRINT:PRINT"xmln : ";xmin 
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1350 PRINT:PRINT"xmax : ";xmax 

1360 LOCATE 1,15 : PRINT"Noch eine Berechnung ( /n>" 
;:INPUT a$ 

1370 IF a* = " " THEN RUN 

1380 IF a*=“n‘ THEN 1390 ELSE 1360 

1390 RETURN 

2100 REM UP Integration nach der Sehnen-Trapez-Forme 
1 

2110 a$=“Sehnen-Trapez-Formel" 

2120 lnteg=0 : x=xmln 

2130 GOSUB 50 : REM Funktionsgleichung 

2140 lnteg = lnteg + y/2 

2150 x=xmax : GOSUB 50 

2160 lnteg = lnteg + y/2 

2170 FOR x=xmin+h TO xmax-h/2 STEP h 

2180 GOSUB 50 

2190 lnteg = lnteg +y 

2200 NEXT x 

2210 lnteg = integ*h 

2220 RETURN 


3.2 Integration nach der Simpsonschen Regel 

Grundlagen 

Wie auch die Trapezregel zählt man die Simpsonsche Regel zu den 
Näherungsverfahren, mit deren Hilfe man Funktionen numerisch 
integrieren kann; gegenüber dieser besteht der Unterschied, daß die 
y-Punkte nicht durch eine Gerade (s. Bild 3.2), sondern durch eine 
Parabel verbunden werden. Außerdem muß die Anzahl der Streifen 
gerade sein. 
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Bild 3.3 Integration nach der Simpsonschen Regel 


Will man die Simpsonsche Formel anwenden, faßt man die Streifen aus 
Bild 3.2 zu Doppelstreifen zusammen. Des weiteren verschiebt man das 
Koordinatensystem so, daß die y-Achse durch den Punkt P x geht. 

Die gesamte Fläche A 0 + wird als A 02 bezeichnet, d.h., die Fläche 
zwischen den Punkten x 0 , y 0 , x 2 und y 2 . 

Nach Bild 3.3 kann man für die drei Punkte P 0 , P, und P 2 durch die 
spezielle Wahl des Koordinatensystems auch schreiben: 

Po = (x 0 ,y 0 ) = (~K y 0 ) 

Pi=(x l ,yi)= (O.yj) 

Pz = (x 2 ,y 2 ) = (+h, y 2 ) 

Gesucht wird jetzt die Parabel, die durch alle drei Punkte geht. Hierzu 
werden folgende Gleichungen aufgestellt: 
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y = a 2 x 2 + a t x + a 0 

(allgemeine Gleichung) 

(3.2.1) 

y 0 =a 2 h 2 -a t h + o 0 

(x= -h) 

(3.2.2) 

y i =a 0 

(x = 0) 

(3.2.3) 

y 2 = a 2 h 2 +a { h + a 0 

(x = h) 

(3.2.4) 

Aus Gleichung 3.2.3 folgt direkt: 



a 0 =yi 

(3.2.5) 


Addition von Gleichung (3.2.2) und Gleichung (3.2.4) ergibt a 2 : 

y 0 + yi = a 2 lj2 - a \ h + Ji + a 2 h 2 + a { h + a 0 
y 0 +y 2 = 2a 2 h 2 + 2y l 


a 2 


y 0 + y2~ iy 1 

2h 2 


(3.2.6) 


Subtraktion von Gleichung 3.2.2 von Gleichung 3.2.4 ergibt o,: 

y 2 ~ yo~ a ih 2 + a x h + y t - a 2 h 2 + a { h - y x 

y 2 -yo = 2 a i h 


«i 


-G .'o 
2h 


(3.2.7) 


A 02 ist aber wiederum nichts anderes als das bestimmte Integral der 
Parabel in Gleichung (3.2.1) in den Grenzen —h und h. 

(Für die folgenden Rechenschritte bzw. Schreibweisen vgl. z.B. [22], 
S. 357 ff) 


a o 2 = J (a 2 x 2 + a,x + a 0 )dx 


1 3 1 J 

= h a 2 x + y x + a 0 x 


a 2 h 3 a,h 2 , a 7 h 3 a,h 2 
= 3+^2~ + a ° h+ ^3 - J~ + “o 1 ' 


-2°?f + 2a 0 l, 
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^02 — 


200 - 2 .)-! + v 2 )h } 

2-3h 2 + yl " 

O'o — 2 >’i + y 2 )h 


3 


+ 2 y x h 


— j O'o + 4y ( + y 2 ) 


Führ! man die gleiche Rechnung für die angrenzende Fläche A 24 durch 
(das Koordinatensystem legt man durch x 3 ), so erhält man die gleiche 
Formel, jedoch sind die Indizes von y jeweils um 2 erhöht. 

-424 = 3 0'2 + 4j- 3 +>’4) 

Das Gleiche gilt für die dann angrenzenden Flächen: 

Ae = 3 0’4 + 4y 5 + y 6 ) usw. 

Die Gesamtfläche unter der Kurve (s. Bild 3.1) errechnet sich als Summe 
der Einzelflächen: 


A — A 02 + A 24 + A 46 + ... + A„_ 2 „ 

Setzt man für die einzelnen Flächen die entsprechenden Gleichungen ein, 
so erhält man: 

A — 2 [O’o + 4>-i + y 2 ) + (y 2 + 4y 3 + y*) +... 

+ O'n- 2 + 4>’„_, + y „)] 

Diese Gleichung liefert nach weiteren Umformungen die Simpsonsche 
Formel: 


h 

A = 3 O’o + 4 y, + ly 2 + 4y 3 + 2\\ + ... + 2y n _ 2 + 4 y n _ t + y n ) 

(3.2.8) 

h 

A — 3 O’o + y n + 40 ’1 + y 3 + ...y„- 1 ) + 20 ' 2 + y 4 + • ■ • + + y n - 2 )] 
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Beispiel 

Um den Unterschied zu verdeutlichen, der zwischen der Trapezregel und 
der Simpsonschen Regel besteht, soll die gleiche Aufgabe wie in 
Abschnitt 3.1 gelöst werden. Es geht also um die Integration der Funktion 

y = 2/X 

in den Grenzen x = 2 und x = 4. Als Anzahl der Streifen wird wieder 10 
gewählt; die Streifenbreite ist somit wieder 0,2. 

Der genaue Wert der Fläche ist 1,38629436 Flächeneinheiten (aus einer 
geschlossenen Integration berechnet). Zur Anwendung der Simpsonschen 
Regel stellt man zweckmäßigerweise folgende Tabelle auf: 


D 

m 

Ti = /(*,■) = 2 Ix 

0 

2,0 

1 



i 

2,2 


0,909090909 


2 

2,4 



0,833333333 

3 

2,6 


0,769230769 


4 

2,8 



0,714285714 

5 

3,0 


0,666666667 


6 

3,2 



0,625 

7 

3,4 


0,588235294 



3,6 



0,555555556 


3,8 


0,526315789 



4,0 

0,5 





1,5 

3,459539428 

2,728174603 



= 1« 

-I* 

= 1^ 


Die Fläche A ist laut Gleichung (3.2.8): 

A = ^=-(1,5 + 4*3,459539428 + 2*2,728174603) 
A = 1,38630046 Flächeneinheiten 
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Trapezregel: 

A = 1,38754281 Flächeneinheiten 

Genauer Wert: 


A = 1,38629436 Flächeneinheiten 

In diesem Fall ist die Simpsonsche Regel also zwei Dezimalstellen genauer 
als die Trapezregel. Diese Aussage kann man verallgemeinern; die 
Simpsonregel ist immer genauer (bzw. in Grenzfällen gleich genau) als die 
Trapezregel. Die Ursache hierfür ist die bessere Anpassungsmöglichkeit 
an die integrierende Funktion durch die Parabel. 


Literaturhinweise 

G: [4], S. 304 bis 305; G: [5], S. 356 bis 360; G: [22], S. 399 bis 402 


Programmbesch reibung 

Das Programm Isimpson ist genauso aufgebaut wie das in Abschnitt 3.1 
beschriebene Programm Itrapez; selbst die Zeilennummern sind gleich. 
Aus diesem Grund sei hier auf eine detaillierte Programmbeschreibung 
verzichtet. Der einzige Unterschied besteht darin, daß im Unterpro¬ 
gramm nach Zeile 1000 nicht die Sehnen-Trapez-Formel, sondern die 
Simpsonsche Regel angewendet wird. 

Programmlänge: etwa 1,4K. 


10 REM Programm Isimpson 

20 REM Integration einer vorgegebenen Funktion mit d 
er Simpson'schen Formel 
30 GOTO 70 

40 REM Zelle 50 = Funktionsgleichung 
50 y=SQR(1+EXP< 2*x)) 

60 RETURN 

70 REM los geht’s 

80 MODE 2 
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90 GOSUB 1000 : REM Startwerteingabe etc. 

100 CLS 

110 LOCATE 20,10 : PRINT’Ich arbeite. Bitte warten." 
120 GOSUB 1100 : REM Berechnung nach Simpson 
130 GOSUB 1300 : REM Ausgabe 
140 END 

1000 REM UP Startwerteingabe etc. 

1010 CLS 

1020 PRINT“In welchen Grenzen soll die Funktion Inte 
griert werden ?" 

1030 PRINT:PRINT*von xmln = ";:INPUT xmln:PRINT"bls 
xmax = INPUT xmax 

1040 PRINT:PRINT"Mlt welcher Schrittweite soll der R 
echner Vorgehen INPUT h 
1050 RETURN 

1100 REM UP Simpson-Berechnung 
1110 a*="Slmpson'sche Formel“ 

1120 integ = 0 : x=xmln 

1130 GOSUB 50 : REM Funktion 

1140 integ = integ + y 

1150 x=xmax : GOSUB 50 

1160 Integ = integ + y 

1170 FOR x=xmln+h TO xmax STEP h*2 

1180 GOSUB 50 

1190 Integ = Integ + y*4 

1200 NEXT x 

1210 FOR x=xmln+2*h TO xmax-h STEP 2*h 

1220 GOSUB 50 

1230 integ = integ + 2*y 

1240 NEXT x 

1250 Integ = integ*h/3 
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1260 RETURN 

1300 REM UP ausgabe 

1310 CLS 

1320 PRINT a* : PRINT"_" : P 

RINT 

1330 PRINT"Integral : "Jinteg 
1340 PRINT:PRINT"xmln : "Jxmln 
1350 PRINT:PRINT"xmax : "Jxmax 

1360 LOCATE 1,15 : PRINT'Noch eine Berechnung < /n>" 
;:INPUT a$ 

1370 IF a*=“" THEN RUN 

1380 IF a$=“n“ THEN 1390 ELSE 1360 

1390 RETURN 





4 

Vektoralgebra 


4.1 Grundrechenarten 

Grundlagen 

Ein Skalar ist eine Größe, die durch eine einzige Angabe eindeutig 
beschrieben ist, wie z.B. Masse, Zeit oder auch Länge. Ein Vektor 
dagegen braucht zur eindeutigen Bestimmung einen Skalar, der den 
Betrag des Vektors angibt, und eine Richtungsangabe, in die der Vektor 
wirkt. Beispiele hierfür sind Kräfte, Momente oder Geschwindigkeiten. 
Nicht genau definiert werden muß der genaue Ort, an dem der Vektor 
beginnt; man kann Vektoren beliebig parallel verschieben. 

In Bild 4.1 sind die beiden Vektoren a und S in ein Koordinatensystem 
eingezeichnet. Um sie eindeutig zu bestimmen, gibt es zwei Möglichkei¬ 
ten: Entweder gibt man den jeweiligen Betrag a = \ä\ bzw. b = \b\ an 
sowie den Winkel a„ bzw. a b , den sie mit der x-Achse einschließen, oder 
aber man gibt die jeweiligen Beträge in x- und y-Richtung für jeden 
Vektor an. Bei der zuerst genannten Möglichkeit spricht man von 
Polarkoordinaten, bei der zweiten von rechtwinkligen Koordinaten. Will 
man den Vektor a z.B. in Polarkoordinaten angeben, so lautet die 
Angabe: 

a = |a| (sina„ + cosa„) (4.1.1) 

In rechtwinkligen Koordinaten lautet die Angabe dagegen: 

ä = aj+aj (4.1.2) 

aj = Vektorkomponente bezüglich der x-Achse 
aj = Vektorkomponente bezüglich der j -Achse 

Polar- und rechtwinklige Koordinaten können ineinander umgerechnet 
werden: 
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a = 

\3\ 

= Ja 2 x + a 2 

(4.1.3) 

<*a = 

arctan ( ^ | 

(4.1.4) 




«x = 

\3\ 

* cos o. a 

(4.1.5) 

“y = 

l«l 

* sina u 

(4.1.6) 


Im Folgenden wird hauptsächlich mit rechtwinkligen Koordinaten 
gerechnet. 

Verschiebt man den Vektor b derart, daß sein Beginn an das Ende von 
Vektor a kommt, so hat man die Möglichkeit geschaffen, beide zu 
addieren. Beim sogenannten Kräfteparallelogramm, kurz Krafteck 
genannt, zieht man einfach eine Linie vom Beginn des Vektors ä zum 
Ende des Vektors b und erhält so den resultierenden Vektor a + b. 



Bild 4.1 Zwei Vektoren in einem Koordinatensystem 
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X 

Bild 4.2 Addition zweier Vektoren 


Um dies auch rechnerisch zu erfassen, braucht man lediglich die 
jeweiligen .v- und /-Werte beider Vektoren zu addieren, um den 
resultierenden Vektor zu erhalten. 

a + f> = (a x + b x ) i + ( a y + b y )j (4.1.7) 

Will man zwei Vektoren dagegen voneinander subtrahieren, zieht man die 
beiden x- und /-Koordinaten voneinander ab. 

a — b = (a x — b x ) t + ( a y - b y )j (4.1.8) 

Will man einen Vektor mit einem Skalar, d.h. mit einer Zahl, 
multiplizieren (z. B.: Die Kraft ist zu verdoppeln), so multipliziert man 
lediglich den Betrag; der Winkel bleibt erhalten. 

2a = 2 |ä| (sina„ + cosa,) (4.1.9) 

Wenn man in rechtwinkligen Koordinaten rechnet, sind demnach die 
Beträge in w- und in /-Richtung mit dem Skalar zu multiplizieren: 

2a = 2aJ+ 2aJ (4.1.10) 
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Beispiele 

1. Der Vektor ä = 3/ + 4/ ist in Polarkoordinaten anzugeben. 

a=\ä\ = Ja\ + a) 

= V^+4 2 
= 5 

oc 0 = arctan 

= arctan 
= 53,1° 

2. Die Vektoren a = 3i + 4/sowie S = 5i + 10 j sind zu addieren und zu 
subtrahieren. 



a + £ = (a x + b x ) i + (a y + b y )j 
= (3 + 5 )f+ (4+10)/ 

= 8T+ 14/ 

ä-S=(a x - b x ) i + (a y - b y )J 
— (3 — 5) /*+ (4-10)/ 

= — 2f-6/ 

3. Der Vektor a = 3/ + 4 j ist mit 2,5 zu multiplizieren. 

2,5a = 2,5 aj + 2,5 aj 
= 2,5*3 1 + 2,5*4/ 

= 7,5 1+ 10/ 


Anmerkung 

Das Programm zu dem gesamten Kapitel sowie die Literaturhinweise 
finden sich im Anschluß an Abschnitt 4.4.3. 
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4.2 Skalarprodukt 

Grundlagen 

Vektoren lassen sich nicht nur durch die beschriebenen Grundrechenar¬ 
ten miteinander verknüpfen, sondern genügen noch weiteren Rechenge¬ 
setzen. Eines dieser Gesetze beschreibt die Bildung des sogenannten 
skalaren Produktes. Das skalare Produkt selbst ist ein Skalar, die 
Richtungsangabe fällt also weg. Berechnet wird es, indem die jeweiligen 
Koordinaten der beiden Vektoren multipliziert werden und hiervon die 
Summe gebildet wird. 

a - S=a x - b x + a y - b y (4.2.1) 

Kennzeichen für das skalare Produkt ist der Multiplikationspunkt 
zwischen den beiden Vektoren. Entstanden ist diese Vorschrift, weil eine 
Beschreibung der Arbeit nötig wurde; diese hängt ab von einer Kraft F 
und dem Weg s, über den die Kraft insgesamt ansteht. 

W=Fs (4.2.2) 

Die Arbeit ist also das Produkt aus den beiden Vektoren Kraft und Weg; 
da sie selbst aber unabhängig von einer Richtung ist, ist sie kein Vektor, 
sondern ein Skalar. 


Beispiel 

Gegeben ist die Kraft F= F x i + Fyj mit F x = 300N und F y = 1000N. 

Gesucht ist die Arbeit, die auf dem Weg von P l = (5 m, 2 m) nach 
P 2 = (10m, lm) verrichtet wird. 

Zunächst wird der Weg .? bestimmt: 

s = R, >7 

= (10 m - 5m)f + (1 m - 2 m)/ 

= 5 m i — 1 mj 
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Somit errechnet sich die Arbeit wie folgt: 

W=Fs 

= F x ■ s x + F y ■ s y 
= 300N • 5m - 1000N • 1 m 
= 500 Nm 

4.3 Vektorprodukt 

Grundlagen 

Das Vektorprodukt der beiden Vektoren a und b beschreibt die Fläche des 
Parallelogramms, das zwischen den beiden Vektoren aufgespannt ist 
(Bild 4.3). 

Kennzeichen der Schreibweise des Vektorproduktes ist das Kreuz 
zwischen zwei Vektoren. Das Ergebnis ist ein neuer Vektor. 

Der Betrag von Ä ist der Flächeninhalt des Parallelogramms. 
Weiterhin ist festgelegt, daß die Richtung von A senkrecht ist zu der 
Richtung von ä und der Richtung von fi\ Ä steht also senkrecht auf dem 
Parallelogramm. 

Zur Berechnung des Vektorproduktes wird angenommen, daß der 
Vektor in einem dreidimensionalen Koordinatensystem vorliegt. 

ä = a x i 4- a y j + a.Ü (4.3.1) 

6 =bj+ hj+bj (4.3.2) 

aj, bj = Vektorkomponente in x-Richtung 
a y l b y J — dito in y-Richtung 
a z i <, b z R = dito in z-Richtung 

Das Vektorprodukt ist jetzt der Wert der folgenden Determinante (zur 
Berechnung von Determinanten vgl. Kapitel 6). 

• r I 

a x a y a. 
b x b y b z 


ax b = 


(4.3.3) 
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Beispiel 


Wie groß ist das Parallelogramm zwischen den Vektoren 

ä=3i + 2j + 0i< und B = 5i + 2j + oBl 

Beide Vektoren haben keine z-Komponente, liegen also in einem 
zweidimensionalen Koordinatensystem. Trotzdem muß nach Gleichung 
(4.3.3) ein dreidimensionales betrachtet werden. 


A = a ■ B = 


i J B 
3 2 0 
5 2 0 


= 0 / + o /— aB 


Der Flächenvektor Ä hat also eine alleinige Komponente in z-Richtung, 
steht somit senkrecht zu a und B. 


\Ä | — \J A\ + Ay + A* 

(vgl. Abschnitt 4.1) 
= v /Ö y To y +4 1 
= 4 


Die Fläche hat eine Größe von 4 Flächeneinheiten. 
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Literaturhinweise 

G: [19], S. 141 ff; G,W: [1], S.67ff; G,W: [5], S. 790ff; G,W: [7], 
S. 170ff; W: [16], S.9ff. 


Programmbeschreibung 

Alle in Kapitel 4 beschriebenen Rechenarten sind im Programm 
Vektoralgebra zusammengefaßt. Das Programm ist hierbei so aufgebaut, 
daß auf dem Bildschirm zunächst ein Menü erscheint; hieraus wählt man 
über eine Code-Zahl die gewünschte Rechenoperation aus. Nach erfolgter 
Ausführung kehrt das Programm dann jeweils zum Menü zurück. 

In 

Zeile 30 erfolgt zunächst die Eingabe der Dimension des zu 
berechnenden Vektorraumes; entsprechend erfolgt in 
Zeile 40 die Dimensionierung der Felder v und fl, in denen die 
Vektoren abgespeichert werden. 

Zeile 50 

bis 180 stellt das erwähnte Menü mit Eingabe der Code-Zahl dar, 
aufgrund der 

Zeile 190 das entsprechende Unterprogramm auswählt. 

Zeile 200 bringt das Programm wieder an den Anfang der Menüaus¬ 
gabe zurück, nachdem die Rechnung durchgeführt wurde. 

Mittels des Unterprogramms in den 

Zeilen 300 

bis 460 wird der erste Vektor eingegeben. 

Dies geschieht zunächst durch das Unterprogramm in 

Zeile 2000 

bis 2150 mit dem die Beträge des Vektors in jc-, y-, z- usw. Richtung 
nach v\ gespeichert werden. In den 

Zeilen 330 

bis 350 wird der Vektor dann nach v gespeichert. Dieser erste 
Programmabschnitt muß in jedem Fall durchlaufen wer¬ 
den, wenn mit Vektoren gerechnet werden soll, da alle 
weiteren Operationen mit dem in v gespeicherten Vektor 
ablaufen. 
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Zeilen 

bis 


Zeile 

bis 


Zeile 


Zeile 

bis 

Zeilen 

bis 


Zeile 

bis 


Zeile 


Auch das Ergebnis aller Rechenschritte wird jeweils nach v 
gespeichert, sofern es einen Vektor darstellt; auf diese Art 
und Weise sind Kettenrechnungen leicht möglich. 

Das etwas umständliche Verfahren der Eingabe wurde 
gewählt, da auch bei den eigentlichen Rechenoperationen 
z.T. Vektoren eingegeben werden müssen. Diese können 
hier über das gleiche Unterprogramm eingegeben werden, 
lediglich die weitere Behandlung ist dann von Fall zu Fall 
unterschiedlich. 

Dies kann man bereits am Beispiel der 

410 

460 sehen, die das Unterprogramm für die Addition und die 
Subtraktion eines Vektors darstellen. Auch hier wird das 
Unterprogramm in 

2000 

2150 für die Eingabe des neuen Vektors benutzt; dieser wird dann 
zu dem in v befindlichen Vektor addiert bzw. von diesem 
abgezogen. Der Ergebnisvektor wird dann wieder nach v 
gespeichert. 

Bei der Faktormultiplikation wird der Faktor in 

520 eingegeben und jeder Betrag von v gemäß der Rechenregel 
mit diesem Faktor multipliziert. Auch hier bleibt der 
Ergebnisvektor in v. Bei der Bildung des Skalarproduktes in 

600 

700 wird zunächst der Multiplikationsvektor eingegeben und 
dann in den 

630 

650 das Skalarprodukt gebildet. Dies wird dann sofort ausgege¬ 
ben, da es keinen Vektor darstellt und somit nicht nach v 
gespeichert werden kann. Durch Betätigen einer beliebigen 
Taste springt man dann zum Hauptprogramm zurück. Zur 
Bildung des Vektorproduktes in 

800 

910 wird zunächst geprüft, ob die betrachteten Vektoren 
dreidimensional sind, da nur dann das Vektorprodukt 
errechnet wird. Ist das nicht der Fall, wird in 

820 der entsprechende Hinweis gegeben und dann zurück zum 
Hauptmenü gesprungen. 
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Zeilen 

bis 

Zeilen 

bis 

Zeile 


Zeilen 

bis 


Zeilen 

bis 


Andernfalls wird der Multiplikationsvektor eingegeben und 
in den 
870 

900 das Vektorprodukt berechnet und nach v gespeichert. In den 

1000 

1370 werden Polar- und rechtwinklige Koordinaten ineinander 
umgerechnet; nachdem man in 

1060 eingegeben hat, was gemacht werden soll, fragt das 
Programm die notwendigen Eingaben ab und gibt die neu 
berechneten Koordinaten dann aus. 

Das vorletzte Unterprogramm in den 

1400 

1510 dient dazu, den in v befindlichen Vektor auszugeben. Hat 
man also eine Rechenoperation durchgeführt, deren Ergeb¬ 
nis ein neuer Vektor ist und das nach v gespeichert ist, wird 
der Vektor nicht automatisch ausgegeben, sondern man 
muß die Ausgabe erst über das Menü anwählen. Auch dies 
erscheint im ersten Moment sicherlich etwas umständlich, 
erleichtert jedoch Kettenrechnungen und dient zudem der 
Übersichtlichkeit. Mit dem Unterprogramm in den 

1600 

1690 wird schließlich entsprechend Abschnitt 4.1 der Betrag des 
in v befindlichen Vektors berechnet. 

Das Programm benutzt bei den Berechnungen ausschließ¬ 
lich rechtwinklige Koordinaten, da diese programmtech¬ 
nisch leichter zu verarbeiten sind. Fehleingaben werden 
durch die Dialogführung jedoch ausgeschlossen. 

Länge des Programms: etwa 5,1 K. 
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10 REM Gesamtprogramm Vektoralgebra 
20 MODE 2 

30 CLS : INPUT"Wievlele Dimensionen hat der Vektorra 
um";d1 mens Ion 

40 DIM v(d1 mens Ion),v1(dl mens Ion) 

50 CLS : REM Hauptmenue 

60 PRINT"Was soll gemacht werden ?":PRINT“_ 


70 PRINT:PRINT"- der 1.Vektor soll eingegeben werden 
“;TAB(70);"(D" 

80 PRINT:PRINT"- ein Vektor soll addiert werden";TAB 
<70>;"<2>" 

90 PRINT:PRINT"- ein Vektor soll subtrahiert werden" 
;TAB(70);"< 3)“ 

100 PRINT:PRINT"- es soll mit einem Faktor multlpllz 
iert werden“,TAB(70)>"(4)" 

110 PRINT:PRINT"- es soll ein Skalarprodukt gebildet 
werden">TAB(70)J"(5)" 

120 PRINT:PRINT"- es soll ein Vektorprodukt gebildet 
werden";TAB(70);" (6) “ 

130 PRINT:PRINT"- Polar- bzw. rechtwinklige Koordlna 
ten sollen umgerechnet werden"»TAB(70);"(7)" 

140 PRINT:PRINT“- der berechnete Vektor soll ausgege 
ben werden“;TAB(70);"(8)" 

150 PRINT:PRINT"- der Betrag des berechneten Vektors 
soll bestimmt werden"5TAB(70)»"(9 ) " 

160 PRINT:PRINT"- Schluss";TAB(70);“<10)" 

170 INPUT was 

180 IF was<1 OR was>10 THEN 50 

190 ON was GOSUB 300,400,400,500,600,800,1000,1400,1 






96 


Vektoralgebra 


600,1700 
200 GOTO 50 

300 REM Eingabe des 1.Vektors 
310 ERASE v : DIM v<dlmenslon> 

320 GOSUB 2000 : REM Eingabe nach vl 

330 FOR k=l TO dimension : REM Umspeicherung nach v 

340 v<k)=vl<k) 

350 NEXT k 
360 RETURN 

400 REM UP/HP fuer + und - 

410 GOSUB 2000 : REM Eingabe des neuen Vektors nach 
VI 

420 FOR k=1 TO dimension 

430 IF was=2 THEN v(k)=v(k) + vl(k) 

440 IF was=3 THEN v(k)=v<k> - vl<k> 

450 NEXT k 
460 RETURN 

500 REM UP/HP Faktormultiplikation 
510 CLS 

520 PRINT’Faktor";:INPUT fak 
530 FOR k=l TO dimension 
540 v<k)=v(k) * fak 
550 NEXT k 
560 RETURN 

600 REM HP/UP Skalarprodukt 

610 GOSUB 2000 : REM Eingabe des Multiplikationsvekt 
ors 

620 sp=0 

630 FOR k=l TO dimension 
640 sp=sp + v<k>*vl<k> 

650 NEXT k 
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660 CLS 

670 PRINT'Skalarprodukt : ">sp 

680 LOCATE 1,23 : PRINT"Weiter mit beliebiger Taste" 
690 IF INKEY*="" THEN 690 
700 RETURN 

800 REM Vektorprodukt 

810 IF dimenslon = 3 THEN 860 

820 CLS : PRINT“Vektorprodukte koennen nur von dreld 
lmenslonalen Vektoren berechnet werden (die z-K 

omponente muss u.U. mit 0 ei 
ngegeben werden)." 

830 LOCATE 1,23 : PRINT’Welter mit beliebiger Taste" 

840 IF INKEY$ ="" THEN 840 
850 RETURN 

860 GOSUB 2000 : REM Eingabe des 2.Vektors 
870 xb = v(2)*vl(3) - vl(2)*v<3) 

880 yb = -v <1>*v1< 3) + vl<l)*v(3) 

890 zb = v<1)*v1< 2) - v1<1)*v< 2) 

900 v(1) = xb : v<2)=yb : v<3) = zb 
910 RETURN 

1000 REM UP/HP Umrechnung von Polar- und rechtwinkli 
gen Koordinaten 
1010 DEG 
1020 CLS 

1030 PRINT'Was soll in was umgerechnet werden ?" 

1040 PRINT:PRINT"- Polarkoordinaten ln rechtwinklige 
Koordinaten"»TAB(55); " (1) " 

1050 PRINT:PRINT"- rechtwinklige Koordinaten ln Pola 
rkoordinaten*;TAB(55);"(2)" 

1060 INPUT wa3 

1070 IF was<1 OR was>2 THEN 1000 
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1080 IF was = 2 THEN GOTO 1250 

1090 REM Umrechnung von Polarkoordinaten ln rechtwln 
klige Koordinaten 

1100 CLS : PRINT"Betrag des Vektors : ";:INPUT betra 
9 

1110 PRINT:INPUT"Winkel zur x-Achse (In Grad) : ";wi 
nke 1 

1120 x=betrag*SIN<winke1) 

1130 y=betrag*COS(wlnkel) 

1140 rcos=COS(wlnkel) 

1150 GOSUB 1170 : REM Ausgabe 
1160 RETURN 
1170 REM Ausgabe-UP 
1180 CLS 

1190 PRINT'x : "»x, “y : ";y 

1200 PRINT:PRINT"Betrag : "Ibetrag,"Winkel : "Jwlnke 
1 *"grad","(eventuel1 nicht richtig, da tan nicht ein 
deut lg)" 

1210 PRINT:PRINT"Richtungscoslnus : "ircos 

1220 LOCATE 1,23 : PRINT'Welter mit beliebiger Taste 

1230 IF INKEY*="* THEN 1230 
1240 RETURN 

1250 REM Umrechnung von rechtwinkligen Koordinaten 1 
n Polarkoordinaten 

1260 CLS : PRINT"x-Betrag : ";:INPUT x 
1270 PRINT:PRINT"y-Betrag : ";:INPUT y 
1280 betrag = SQR(x*x + y*y> 

1290 Winkel = ATN(y/x) 

1300 rcos = x/betrag 

1310 GOSUB 1170 : REM Ausgabe 
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1320 RETURN 

1400 REM Ausgabe des Vektors 
1410 CLS 

1420 FOR k=1 TO dimenslon 

1430 IF k=1 THEN LOCATE 12,k : PRINT"x : 

1440 IF k=2 THEN LOCATE 12,k : PRINT'y : “J 

1450 IF k=3 THEN LOCATE 12,k : PRINT"z : ", 

1460 IF k>3 THEN PRINT*Dimens 1on "Jk;": 

1470 PRINT v(k) 

1480 NEXT k 

1490 LOCATE 1,23 : PRINT”We1ter mit beliebiger Taste 

1500 IF INKEY*="" THEN 1500 
1510 RETURN 

1600 REM Betragsberechnung eines Vektors 
1610 sum=0 

1620 FOR k=l TO dimenslon 
1630 sum=sum + v<k>*v<k> 

1640 NEXT k 

1650 betrag=SQR<sum) 

1660 CLS : PRINT’Betrag : ’Jbetrag 

1670 LOCATE 1,23 : PRINT'Welter mit beliebiger Taste 

1680 IF INKEY* = ” THEN 1680 
1690 RETURN 

1700 REM UP/HP Schluss machen 
1710 CLS 

1720 LOCATE 20,10 : PRINT’Na gut. Denn man Tschuess. 

N 

1730 LOCATE 1,23 
1740 END 
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2000 REM Eingabe des Vektors ln das Feld vl 
2010 ERASE vl : DIM v1(dimens Ion) 

2020 CLS : PRINT'Geben Sie den Vektor ein. Es werden 
rechtwinklige Koordinaten benoetigt, d.h., die Bet 
raege des Vektors ln x-, y-, 
z-, etc. Richtung":PRINT"_ 


2030 FOR k=1 TO dimenslon 

2040 LOCATE 1,6 

2050 IF k=1 THEN PRINT"x : "; 

2060 IF k=2 THEN PRINT"y : "; 

2070 IF k=3 THEN PRINT"z : "J 

2080 IF k>3 THEN PRINT’Dimension "JkJ" : "5 

2090 INPUT vl(k) 

2100 LOCATE 1,6 : PRINT" 

2110 NEXT k 

2120 LOCATE 1,23 : PRINT’EIngaben alle i.o. < /n>";: 
INPUT a* 

2130 IF a*="" THEN 2150 

2140 IF a*="n" THEN 2000 ELSE 2120 

2150 RETURN 
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Komplexe Zahlen 


Grundlagen 

In Naturwissenschaft und Technik ist es bei vielen Gelegenheiten 
wünschenswert, Gleichungen der Form x 2 = — 1 zu lösen. Dies ist in der 
Menge der reellen Zahlen nicht möglich. Um sich aus dieser Verlegenheit 
zu helfen, hat man die imaginäre Einheit i definiert: 

i 2 = -l (5.1) 

Durch diese Definition entsteht eine ganze Klasse neuer Zahlen, die 
komplexen Zahlen. Zum Rechnen mit i verfährt man grundsätzlich wie 
bei reellen Zahlen; so gilt z.B. 

i + i =2i 3,5 i — i = 2,5 i 

i° = 1 i 1 = i 

i 2 =-l i 3 = i - i 2 = — 1 

Hieraus kann man folgende Regeln herleiten: 

l=i 4 " (5.2) i = i 4n+1 (5.3) 

— 1 = i 4 " + 2 (5.4) -i = i 4n+3 (5.5) 

wobei n ein Element der natürlichen Zahlen einschließlich 0 sein kann. 

Eine komplexe Zahl Z besteht aus einem reellen und einem 
sogenannten imaginären Anteil; veranschaulichen kann man sich dies mit 
der Gaußschen Zahlenebene, wo in der Waagerechten der reelle Teil von 
Z und in der Senkrechten der imaginäre Teil aufgetragen wird (Bild 5.1). 
Definiert ist die komplexe Zahl als Punkt in dieser Ebene: 

Z = a + bi (5.6) 

a = Real teil von Z 
bi = Imaginärteil von Z 
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s 



Bild 5.1 Die komplexe Zahl in der Gaußschen Zahlenebene 


Verschwindet bi, so hat man eine reelle Zahl vorliegen, die Zahl befindet 
sich auf der waagerechten Achse. Verschwindet dagegen a, so hat man 
eine sogenannte imaginäre Zahl; diese liegt auf der senkrechten Achse. 
Eine «größer als» oder «kleiner als» Beziehung läßt sich bei komplexen 
Zahlen nicht definieren; dagegen lassen sich die vier Grundrechenarten 
anwenden. Hierbei wird auch die Verwandschaft mit den in Kapitel 4 
besprochenen Vektoren deutlich. 

Addition 

Zwei komplexe Zahlen werden addiert, indem man jeweils die Realteile 
und die Imaginärteile addiert. 

(a, +b 1 \) + (a 2 + b 2 i) = (a, +a 2 ) + (/>, +b 2 ) i (5.7) 

Dies entspricht genau der Addition zweier Vektoren; das gilt auch für die 
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Subtraktion 

Zwei komplexe Zahlen werden subtrahiert, indem man jeweils die 
Realteile und die Imaginärteile voneinander subtrahiert. 

(«, + b x i) - (a 2 + * 2 ») = («i - a 2 ) + (b x - b 2 )i (5.8) 


Multiplikation 

Bei der Multiplikation zweier komplexer Zahlen werden beide gliedweise 
miteinander multipliziert. 

(a x + b x i)*(a 2 +/> 2 >) = a x a 2 + a x b 2 i + a 2 b x \ + b x b 2 i 2 
Da i 2 aber gleich —1 ist, kann man hierfür auch schreiben: 

(a x + b x i) * (a 2 + b x i) = (a x a 2 - b x h 2 ) + (a x b 2 + a 2 b x ) '\ (5.9) 


Division 

Will man die Zahl a x + b x i durch die Zahl a 2 + b 2 i dividieren, so 
erweitert man den Bruch zunächst mit dem Konjugierkomplex a 2 - b 2 \ 

a x +b x \ _ (a x +b x i)(a 2 -b 2 i) 
a 2 + b 2 i “ (a 2 +b 2 \)(a 2 - b 2 \) 

Jetzt wendet man auf Zähler und Nenner jeweils die Multiplikationsregel 
an: 

Zähler: ( a x + b x i) ( a 2 - b 2 i) = (a x a 2 + b x b 2 ) — (a x b 2 — a 2 b x )i 
Nenner: ( a 2 +b 2 i) (a 2 -b 2 i) = (a\ + b\) - (a 2 b 2 - a 2 b 2 )'\ = a\ + b\ 
Somit erhält man als Rechenregel: 

a x +b x \ a x a 2 + b x b 2 , -a x b 2 + a 2 b x . 
a 2 + b 2 i « 1+^2 fl 2 + Ä 2 

Komplexe Zahlen kann man genau wie Vektoren auch in einer Polarform 
angeben. Hierfür dient die sogenannte Zeigerdarstellung (Bild 5.2). 

Die Polarform lautet dann: 


Z = r(sin c^> + icoscp) 


(5.11) 
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Bild 5.2 Zeigerdarstellung der komplexen Zahl Z 


Die bisher verwendete Normalform Z = a + b i und die Polarform lassen 


sich nach den gleichen Gesetzen wie bei Vektoren ineinander umrechnen 
(vgl. Abschnitt 4.1) 


(5.12) 

/ a\ 

(p = arctan 1 - 1 

(5.13) 

a = r ■ cos (p 

(5.14) 

b = r • sin cp 

(5.15) 

Beispiel 


Auf die beiden komplexen Zahlen 


Z = 5 + 7i 

Z= -3 - 5i 

und 


sind die vier Grundrechenarten anzuwenden; des weiteren ist Z in die 
Polarform umzuwandeln. 
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(5 + 7i) + (- 
(5 + 7i) — ( — 
(5 + 7i) * ( — 

(5 + 7i)/( — 


3 — 5i) = 2 + 2i 
3 — 5i) = 8 + 12i 

3 — 5i) = (—15 + 35) + (— 25 — 21)i = 20 - 46i 

-15 - 35 25 - 21 . 50 4 . 

3 31 — 9 + 25 + 9 + 25 1 ~ _ 34 + 34 1 

Z, = 5 + 7i 


r = J25 + 49 = 8,602 


cp = arctan 



35,54° 


Die Polarform von Z, lautet also: 

Z ( = 8,602 (sin 35,54° + icos 35,54°) 


Literaturhinweise 

G: [19], S. 16 bis 20; G,W: [5], S. 400 bis 408; G,W: [13], S. 130 bis 146; 
G, W; [22], S. 151 ff; W: [19], S. 127ff; W: [23] 


Programmbeschreibung 

Die Struktur des Programms komplex ist mit dem in Kapitel 4 
beschriebenen Programm Vektoralgebra identisch. Neben den u.U. 
anderen Rechenvorschriften besteht ein wesentlicher Unterschied zwi¬ 
schen Programmen eigentlich nur darin, daß das Programm komplex nur 
zweidimensionale Rechnungen durchführen muß, während die Anzahl 
der Dimensionen des Vektorraumes im Programm Vektoralgebra frei 
wählbar ist. Das Programm hat demzufolge auch nur eine Länge von 
etwa 3,3 K, während Vektoralgebra etwa 5,1 K an Speicherplatz benötigt. 

In den 

Zeilen 50 

bis 180 ist das Menü abgelegt, über das die verschiedenen Rechen¬ 
arten angewählt werden können, die mit komplexen Zahlen 
möglich sind. Über die 

Zeilen 300 

bis 330 wird die erste komplexe Zahl eingegeben; die 
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Zeilen 400 

bis 440 werden bei Addition/Subtraktion und die 
Zeilen 500 

bis 550 bei Multiplikation aufgerufen. 


Zeile 600 

bis 660 ist schließlich zuständig für die Division von komplexen 
Zahlen. In den 

Zeilen 700 

bis 1070 werden Polarform und Normalform ineinander umgerech¬ 
net; dieser Programmabschnitt ist praktisch identisch mit 
dem entsprechenden Abschnitt im Programm Vektoral¬ 
gebra. Über die 

Zeilen 1100 

bis 1160 kann zum Schluß noch die komplexe Zahl, die das Ergebnis 
der Rechnungen darstellt, ausgegeben werden. 


10 REM Programm komplex 

20 REM Rechenoperationen mit komplexen Zahlen 
30 MODE 2 

40 REM faellt aus wegen Bodennebel 
50 CLS : REM Hauptmenue 

60 PRINT’Was soll gemacht werden*; :PRINT"_ 


70 PRINT:PRINT“- die 1.komplexe Zahll soll eingegebe 
n werden"J TAB(60);" < 1)" 

80 PRINT:PRINT"- eine komplexe Zahl soll addiert wer 
den";TAB(60);*(2)" 

90 PRINT:PRINT"- eine komplexe Zahl soll subtrahiert 
werden";TAB<60>;"<3>" 

100 PRINT:PRINT"- es soll mit einer komplexen Zahl m 
ultlpliziert werden";TAB<60);"<4)" 
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110 PRINT:PRINT"- es soll durch eine komplexen Zahl 
dividiert werden" ; TAB<60);"<5) * 

120 PRINT:PRINT"- Polarform und Normalform sollen um 
gerechnet werden";TAB<60);“(6)* 

130 PRINT:PRINT"- die berechnete Zahl soll ausgegebe 
n werden*»TAB<60)*"<7>" 

140 PRINT:PRINT"- Schluss";TAB(60);"(8) " 

150 INPUT was 

160 IF was <1 OR was>8 THEN 50 

170 0N was GOSUB 300,400,400,500,600,700,1100,1200 

180 GOTO 50 

300 REM HP/UP Eingabe der l.Zahl 
310 GOSUB 1500 : REM Eingabe nch .2 
320 x1=x2 : yl=y2 
330 RETURN 

400 REM HP/UP fuer + und - 
410 GOSUB 1500 : REM Eingabe der 2. Zahl 
420 IF was = 2 THEN xl=xl+x2 : yl=yl+y2 
430 IF was = 3 THEN xl=xl-x2 : y1=y1-y2 
440 RETURN 

500 REM HP/UP Multiplikation 

510 GOSUB 1500 : REM MultlpllkatoreIngabe 

520 x3 = xl*x2 - yl*y2 

530 y3 = xl*y2 + yl*x2 

540 xl=x3 : yl=y3 

550 RETURN 

600 REM HP/UP Division 

610 GOSUB 1500 : REM Eingabe des Divisors 
615 y2=-y2 

620 x3=xl*x2 - + y1*y2 
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630 y3 = xl*y2 + yl*x2 
640 nenner = x2*x2 + y2*y2 
650 xl=x3/nenner : y1=y3/nenner 
660 RETURN 

700 REM Umrechnungen von Polarform und Normal form 
710 CLS 

720 PRINT’Bel einer Umrechnung von der Normalform ln 
die Polarform muessen die Faktoren der Normalform 
dem Programm bereits ln xl 
und yl vorliegen.* 

730 PRINT:PRINT’Bei der umgekehrten Rechnung wird dl 
e berechnete Normalform nach xl und yl gespelch 

ert.":PRINT*_ 


740 PRINT:PRINT“Was soll umgerechnet werden ?“ 

750 PRINTsPRINTsPRINT"- Normalform ln Polarform < 
1 >* 

760 PRINTsPRINT*- Polarform ln Normalform (2)" 

770 INPUT was 

780 IF was=1 THEN 900 

800 REM Umrechnung Polarform ln Normal form 
810 CLS 
820 DEG 

830 INPUT’Betrag : “»betrag 
840 PRINT:INPUT’Wlnkel : *;wlnkel 
850 xl=betrag*COS<wlnkel) 

860 yl=betrag*SIN<wlnkel) 

870 GOSUB 1000 : REM Ausgabe 
880 RETURN 
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900 REM Umrechnung Normal form ln Polarform 
910 DEG 

940 betrag = SQR(xl*xl + yl*yl) 

950 Winkel = ATN(yl/xl> 

960 GOSUB 1000 : REM Ausgabe 
970 RETURN 

1000 REM Ausgabe fuer Umrechnungen Polar ln Normalfo 
rm und umgekehrt 
1010 CLS 

1020 PRINT-Realte11 : *;x1,■Imaglnaerte 11 : " J y 1 
1030 PRINT:PRINT“Betrag : “ibetrag,"Phi : “»Winkel»" 
grad" 

1040 PRINT:PRINT“(Phi u.U. nicht richtig, da atan nl 
cht eindeutig)“ 

1050 LOCATE 1,23 : PRINT“Welter mit beliebiger Taste 

M 

1060 IF INKEY*=““ THEN 1060 
1070 RETURN 

1100 REM HP/UP Ausgabe einer komplexen Zahl 
1110 CLS 

1120 PRINT’Realte 11 : “Jxl 

1130 PRINT’Imaglnaerte 11 : “,yl 

1140 LOCATE 1,23 : PRINT’Weiter mit beliebiger Taste 

1150 IF INKEY$="“ THEN 1150 
1160 RETURN 

1200 REM HP/UP Schluss 
1210 CLS 

1220 LOCATE 30,10 : PRINT“Auf Wiedersehen“ 

1230 LOCATE 1,23 
1240 END 
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1500 REM UP Eingabe einer komplexen Zahl 
1510 CLS 

1520 INPUT"Realteil : "1x2 
1530 INPUT"Imaglnaerte11 : ";y2 
1540 RETURN 
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6.1 Eigenschaften von Determinanten 

Grundlagen 

Eine zweireihige Determinante ist dadurch definiert, daß dem quadrati¬ 
schen Zahlenschema 

«11 «12 
«21 «22 

der Zahlenwert 

= «1 1 «22 — «12 «2 1 ( 6 . 1 . 1 ) 

«21 «22 

zugeordnet wird. a ik sind die Elemente der Determinante, wobei i die Zeile 
und k die Spalte angibt. Die Linie öj , a 22 nennt man die Hauptdiagonale, 
die Linie a l2 a 2l Nebendiagonale. 

Aus obiger Definition kann man die Richtigkeit der folgenden Regeln 
herleiten: 


Regel 1 


Der Wert einer Determinante bleibt unverändert, wenn man an der 
Hauptdiagonalen spiegelt. 


«11 

«12 


«11 

«21 

«21 

«2 2 


«12 

«22 


«22 «12 «21 


Das Spiegeln an der Hauptdiagonalen ist auch bekannt als Stürzen; 
Zeilen und Spalten werden hierbei einfach vertauscht. Somit ist auch 
bewiesen, daß alle für Zeilen hergeleiteten Eigenschaften auch für Spalten 
gelten und umgekehrt. 
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Regel 2 

Eine Determinante wird mit einem Faktor multipliziert, indem man die 
Elemente einer beliebigen Zeile oder Spalte mit ihm multipliziert. 

k * {a 11 Ö 22 — «21 « 12 ) = ^«i i «22 — ka 2l a l2 
ka n a l2 
ka 2l a 22 



Regel 3 

Der Wert einer Determinante bleibt unverändert, wenn man zu einer Zeile 
ein beliebiges Vielfaches einer anderen Zeile addiert. 


«ii + ta 2l 
«21 


a, 2 +ta 22 / \ / \ 

= (a n + ta 2l )a 22 -a 21 (a l2 + ta 22 ) 

«22 

= «1 1 «22 " t " *«21 «22 — «21 «12 — *«22 «21 


— « 11^22 «21 «12 


Regel 4 

Der Wert einer Determinante ist gleich 0, wenn eine Zeile ein Vielfaches 
einer anderen Zeile ist. 

= «11 £«12 - ka u a, 2 = 0 

Wenn alle Elemente der Determinante gleich 0 sind, kann man ebenfalls 
obige Regel anwenden. 


«11 «12 
kau ka l2 


Regel 5 

Eine Determinante, die in der Hauptdiagonalen nur Einsen, sonst nur 
Nullen enthält, hat den Wert 1. 

= 1*1 - 0*0 = 1 

Alle bisher hergeleiteten Regeln gelten nur für zweireihige Determinanten. 


1 1 
0 1 
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Genauso gibt es aber auch drei- und mehrreihige Determinanten: 


11 

«12 

«13 

21 

«2 2 

«2 3 


= dreireihige Determinante 


«31 «32 a 33 


Bei drei- bzw. mehrreihigen Determinanten hat man sogenannte 
Unterdeterminanten definiert, die jeweils einem Element zugeordnet sind; 
die zu einem Element a ik gehörige Unterdeterminante erhält man durch 
Streichen der /-ten Zeile und Ä-ten Spalte. Diese Definition wird benutzt, 
um den Wert einer mehrreihigen Determinante zu bestimmen. Hierzu 
werden nacheinander alle Elemente einer beliebigen Zeile i oder Spalte k 
genommen, das Produkt dieser Elemente mit der jeweiligen Unterdeter¬ 
minante gebildet, wobei das Vorzeichen bei + anfangt und dann bei 
jedem Glied wechselt, und schließlich die Summe gebildet. 


«11 

«12 

«13 


«2 2 

«23 


«21 

«23 

«21 

«22 

«23 

= «11 



- «12 



«31 

«32 

«33 


«32 

«33 


«31 

«33 


+ «i. 


Dieses etwas komplizierte Verfahren nennt man Entwickeln einer 
Determinante nach der Spalte k bzw. Zeile i; als Ergebnis erhält man 
Determinanten, die jeweils eine Reihe weniger besitzen als die Ursprungs¬ 
determinante. Man muß also so oft entwickeln, bis nur noch zweireihige 
Determinanten vorhanden sind, von denen man nach (6.1.1) den Wert 
bestimmen kann. Ein anderes Verfahren der Eigenwertbestimmung ist in 
Abschnitt 6.2 beschrieben. 


Regel 6 

Alle für 2reihige Determinanten aufgestellten Sätze gelten auch für die 
drei- und mehrreihigen Determinanten. Für den Beweis dieses Satzes sei 
auf die Literatur [9] verwiesen. 
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Regel 7 

Hat eine Determinante oberhalb (und/oder unterhalb) der Hauptdiago¬ 
nalen nur Nullen, so ist ihr Wert gleich dem Produkt der diagonalen 
Elemente. 


«11 

«12 

«13 


«22 

«23 


«12 

«13 


«12 

«13 

0 

«22 

«2 3 

= «11 



-0 



+ 0 



0 

0 

«33 


0 

«33 


0 

«3 3 


«22 

«23 


— «n(«22«33 0 * «23) — «1 1 «22«33 


Die Entwicklung wird zweckmäßigerweise nach der 1. Spalte oder 3. Zeile 
vorgenommen, alle anderen Entwicklungen führen aber zum gleichen 
Ergebnis. 


Beispiel 

Der Wert folgender Determinante soll berechnet werden: 


2 4-2 


-1 5 


6 5 


6 -1 

6-1 5 

= 2 


- 4 


- 2 


3-2 6 


-2 6 


3 6 


3 -2 


= 2( —1 *6 + 2*5) — 4(6*6 - 3*5) - 2( —2*6 +1*3) 
= 2*4 - 4*21 + 2*9 
= -58 


Literaturhinweise 

G, W: [5], S. 116ff; G,W: [9], S.81 ff; G,W: [13], S.281 ff; G,W: [16], 
S. 87 ff; G, W: [19], S. 37 ff; G,W: [22], S.84IT. 

6.2 Eigenwertbestimmung 

Grundlagen 

Das Entwicklungsverfahren zur Bestimmung des Eigenwertes einer 
Determinante wird leicht unübersichtlich, wenn es sich um mehrreihige 
Determinanten handelt; hier wendet man besser den sogenannten 
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Gaußschen Algorithmus an. Die Bestimmung des Eigenwertes geht dann 
folgendermaßen vor sich: 



«u 

«12 •• 

• «Im 

Determinante = 

«21 

«22 •• 

• «2m 


«nl 

«n2 

«nm 


Die erste Zeile wird mit — « u /«,, multipliziert und zur k-ten Zeile addiert 
(k = 2,3,...,«) (1 fehlt!). Der Wert der Determinante bleibt hierbei nach 
Regel 3, Abschnitt 6.1, unverändert. Das Ergebnis ist folgende Determi¬ 
nante: 


«n a n ■■■ a im 

0 «2 2 ••• a 2m 

o a'„ 2 ... a' nm 

Anschließend nimmt man die zweite Zeile mit —a' k2 /a' 22 mal und addiert 
sie zur &-ten Gleichung (k = 1,3,.... n) (jetzt fehlt 2!). Dieses Verfahren 
wendet man «-mal an. Das Ergebnis ist die Determinante 

«*! 0 ... 0 
0 a* 22 ... 0 

0 0 a* 

Der Wert dieser Determinante ist aber nach Regel 7, Abschnitt 6.1 gleich 
dem Produkt der Elemente der Hauptdiagonalen. 


Beispiel 

Zur Verdeutlichung soll die gleiche Determinante wie im Abschnitt 6.1 
durchgerechnet werden. 


2 4-2 

6-1 5 

3 -2 6 


n = 3 
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1. Schritt: 


2. Schritt: 



fu = 

_^2 


«n 


2 

4 

_ 2 

0 

-13 

11 

3 

-2 

6 


a k\ _ 



«11 


2 

4 

_ 2 

0 

-13 

11 

0 

-8 

9 


a k2 _ 



«22 

ö 

2 

0 

18 


13 

0 

-13 

11 

0 

-8 

9 


«lc2 _ 

a 


«22 

a 

2 

0 

18 

13 

0 

-13 

11 

0 

0 

29 

13 


k = 1: -=?£=- 


k = 3: -^1= - 


6 

2 


= -3 


3 

2 


4 _ 4 
-13 ~ 13 


-8 _ 8 
- l3 ~ _ 13 


An dieser Stelle könnte man bereits aufhören, da die Bedingungen für die 
Regel 7 erfüllt sind. Zur Verdeutlichung wird jedoch der dritte Schritt 
auch noch durchgeführt. 


11 

13 

29 

13 


18 

29 


0 

0 


0 0 
-13 11 


0 


29 

13 


3. Schritt: 


k = 1: 
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«23 


11 

, 

— 

29 

«33 


1 3 

0 

0 


-13 

0 



0 0 fl 


143 

29 


Der Eigenwert der Determinante ist hiernach also 


2 * ( — 13) ♦ 


29 

13 


= -58 


Wie man sieht, ist dieses Verfahren für dreireihige Determinanten sehr 
viel komplizierter als das Entwicklungsverfahren; da es aber einem ganz 
klaren Schema folgt, läßt es sich leicht auf einem Rechner anwenden. 


Literaturhinweise 

G, W: [5], S. 116ff; G,W: [9], S. 81 ff; G,W: [13], S.281 ff; G,W: [16], 
S. 87 ff; G, W: [19], S. 37 ff; G,W: [22], S. 84 ff. 


Program mhesch reibung 

Auch für das Programm determinante wurde die in den 
Kapiteln 4 und 5 beschriebene Programmstruktur beibehal¬ 
ten, obwohl es hier vielleicht aufgrund der geringen Anzahl 
von Rechenoperationen nicht notwendig gewesen wäre; es 
wird jedoch auf diese Weise eine bessere Übersichtlichkeit 
erreicht. 

Nach Eingabe der Determinantengröße in 
Zeile 30 und der Dimensionierung des Hauptfeldes dadd und des 
Berechnungsfeldes d (vgl. Kapitel 4 und 5) in 
Zeile 40 erfolgt in 
Zeile 50 

bis 200 die Ausgabe des Menüs. Im Gegensatz zu den erwähnten 
Programmen erfolgt die Eingabe der Codezahlen hier nicht 
direkt; vielmehr werden Zeichen eingegeben, die in logischer 
Verbindung zu der durchzuführenden Rechnung stehen. 
Diese Zeichen werden in den 
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Zeilen 

bis 

Zeile 


Zeilen 

bis 


Zeile 

bis 


Zeile 

Zeile 


Zeile 

bis 


Zeilen 

bis 


Zeile 

bis 

Zeilen 

bis 


140 

180 dann in Codezahlen umgewandelt, so daß 
190 das entsprechende Unterprogramm auswählen kann. Da 
sich die Unterprogramme für Eingabe und Ausgabe ähneln, 
werden sie im Unterprogramm in den 
300 

400 zusammen ge faßt; wo Unterscheidungen notwendig sind, 
geschieht dies aufgrund der Codezahl. Entsprechend erfolgt 
in 
500 

560 die Multiplikation mit einem Faktor. Obwohl hier prinzi¬ 
piell eine beliebige Spalte oder Zeile der Determinante mit 
dem Faktor multipliziert werden kann, wählt das Pro¬ 
gramm immer die erste Spalte; soll dies anders sein, müssen 
die entsprechenden Indizes in 
540 geändert werden. Ab 

700 erfolgt die Eigenwertbestimmung der Determinante, die in 
dadd gespeichert ist. Um diese für weitere Rechnungen noch 
zur Verfügung zu haben, erfolgt zunächst in 
710 

750 die Umspeicherung nach d. 

Die Berechnung der Lösungsdeterminanten, in der nur noch 
in der Hauptdiagonalen Werte stehen, erfolgt in den 

1100 

1200 Die Lösung erhält man durch Multiplikation der einzelnen 
Faktoren der Hauptdiagonalen. 

Dies wird in 

1220 

1260 durchgeführt. Die 
800 

1040 stellen Unterprogramme dar, die von den anderen Teilpro¬ 
grammen benutzt werden; ihre Bedeutung kann den 
Kommentarzeilen entnommen werden. 

Programmlänge: etwa 2,7K. 
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10 REM HP Determinantenrechnungen 
20 MODE 2 

30 PRINT’Wievleie Zeilen bzw. Spalten hat (haben) di 
e Determinante(n) “>:INPUT danz 
40 DIM d(danz,danz),dadd(danz,danz) 

50 REM Hauptmenue 
60 CLS 

70 PRINT’Was soll gemacht werden ?“ : PRINT "_ 


80 PRINT:PRINT“- die 1.Determinante soll eingegeben 
werden“;TAB(70)>"<1)" 

90 PRINT:PRINT“- es soll mit einem Faktor multlpllzl 
ert werden";TAB<70>;"<*)" 

100 PRINT:PRINT“- die berechnete Determinante soll a 
usgegeben werden";TAB<70)>"(a)" 

110 PRINT:PRINT“- der Eigenwert der berechneten Dete 
rmlnante soll bestimmt werden";TAB<70)(e>" 

120 PRINT:PRINT"- es soll endlich Schluss gemacht we 
rden";TAB(70)*"(s)“ 

130 INPUT a* 

140 IF a*="1" THEN was =1 

150 IF a*="*" THEN was = 2 

160 IF a*="a" THEN was = 3 

170 IF a$="s" THEN was = 4 

180 IF a$="e“ THEN was = 5 

190 ON was GOSUB 300,500,300,600,700 

200 GOTO 50 

300 REM HP/UP fuer 1 und a 
310 IF was = 3 THEN 330 

320 GOSUB 900 : REM Eingabe der Determinante 
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330 CLS 

340 FOR k=1 TO danz 

350 FOR 1=1 TO danz 

360 ON was GOSUB 810,1320,830 

370 NEXT 1 

380 NEXT k 

390 IF was <> 3 THEN 420 

400 LOCATE 1,23 ; PRINT'Welter mit beliebiger Taste“ 
410 IF INKEY*=““ THEN 410 
420 RETURN 

500 REM UP/HP Multiplikation mit einem Faktor 
510 CLS 

520 INPUT “Faktor : “;fak 

530 FOR k=1 TO danz 

540 dadd(k,1) = dadd(k,1)*fak 

550 NEXT k 

560 RETURN 

600 REM UP/HP Ende 

610 CLS 

620 LOCATE 20,10 : PRINT“Ich hab auch keine Lust meh 
r“ 

630 LOCATE 1,23 
640 END 

700 REM UP/HP Eigenwertbestimmung 
710 FOR k=1 TO danz 
720 FOR 1=1 TO danz 
730 d(k,l) = dadd(k, 1) 

740 NEXT 1 
750 NEXT k 

760 GOSUB 1100 : REM Eigenwertbestimmung 
770 CLS : PRINT'Elgenwert : “;dew 
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780 RETURN 

800 REM UP's fuer 1 und a 

810 dadd(k,l> = d<k,l> : REM Umspeicherung 
820 RETURN 

830 LOCATE 1*10-9 , k : PRINT dadd(k,l> 

840 RETURN 

900 REM UP Determinanteneingabe 
910 WINDOW #1,1,80,5, 10 
920 CLS 

930 PRINT'Geben 31e die Determinante zeilenweise ein 

M 

940 FOR k=1 TO danz 

950 LOCATE 1,3 : PRINT k;".te Ze1le":PRINT 
960 CLS#1 

970 FOR 1=1 TO danz 

980 PRINT#1,"d<k1;"> : -;:INPUT#1, d(k,l) 

990 NEXT 1 
1000 NEXT k 

1010 LOCATE 1,20 : INPUT'Elngabe l.o ( /n> -Ja* 

1020 IF a* = "" THEN 1040 

1030 IF a$="n“ THEN 900 ELSE 1010 

1040 RETURN 

1100 REM UP Eigenwertbestimmung 

1110 REM verwendet wird der Gauss’sche Algorithmus 

1120 FOR m=1 TO danz 

1130 FOR k=1 TO danz 

1140 IF k=m THEN 1190 

1150 a1=-d(k,m)/d(m,m) 

1160 FOR 1=1 TO danz 

1170 d<k,1>=d(k,1) + d(m,1)*al 

1180 NEXT 1 
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1190 NEXT k 
1200 NEXT m 

1210 REM Loesungsmatrlx Ist fertig 

1220 REM dew wird durch Multiplikation der Hauptdlag 

onalen bestimmt 

1230 LET dew=1 

1240 FOR 1=1 TO danz 

1250 LET dew=dew*d(1,1) 

1260 NEXT 1 
1270 CLS 

1280 PRINT-Elgenwert : "Jdew 

1290 LOCATE 1,23 : PRINT’We1ter mit beliebiger Taste 

N 

1300 IF INKEY*=“" THEN 1300 
1310 RETURN 
1320 REM 
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Grundlagen 

Mit den Determinanten verwandt sind die Matrizen; auch hier handelt es 
sich um Zahlenschemata, die jedoch nicht mehr zwangsläufig quadratisch 
sein müssen. Das Produkt aus der Zeilenanzahl i und der Spaltenanzahl k 
nennt man auch Grad der Matrix. 

grad = (/ * k) 


Gleich sind zwei Matrizen nur dann, wenn sie in allen ihren Koeffizienten 
übereinstimmen. Im Gegensatz zu Determinanten lassen sich Matrizen 
z. B. durch Addition miteinander verknüpfen: 



«12 

«m2 



... a Xn 

\ / “ 

b i2 



+ : 



••• U mn 

/ 

b m 2 


+ h x , 

«12 + ^12 •• 

■ “ln 

+ b] 

+ (’ml 

a m2 + b m 2 ■ ■ 

amn 

+ b, 



(7.1) 


Matrizen werden also koeffizientenweise addiert; das gleiche gilt für die 
Subtraktion. Daraus folgt auch, daß Matrizen nur dann addiert bzw. 
subtrahiert werden können, wenn sie in Zeilen- und Spaltenzahl 
übereinstimmen. Ähnliches gilt für die Multiplikation einer Matrix mit 
einem Faktor: 



(7.2) 
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Auch die Faktormultiplikation wird demnach koeffizientenweise durch¬ 
geführt. 

Komplizierter wird es, wenn zwei Matrizen multipliziert werden sollen. 
Dies ist nur möglich, wenn die Spaltenanzahl k der ersten Matrix mit der 
Zeilenanzahl i der zweiten Matrix übereinstimmt. Wenn r die Zeilenan¬ 
zahl der ersten Matrix und s die Spaltenanzahl der zweiten Matrix 
bezeichnet, so ist das Ergebnis eine (r*.?)-Matrix. 



Das Glied c,, der Ergebnismatrix erhält man, indem man jeweils das 
Produkt bildet aus den Koeffizienten der ersten Zeile von Matrix a mit 
den entsprechenden Koeffizienten der ersten Spalte von Matrix b; c,, ist 
die Summe aller Produkte. 

Für c 12 multipliziert man die Elemente der ersten Zeile von Matrix a 
mit den Elementen der zweiten Spalte von Matrix b und bildet die Summe. 
Entsprechend verfährt man mit den übrigen Gliedern der Ergebnisma¬ 
trix; für c iA z.B. nimmt man die Elemente der dritten Zeile von Matrix a 
und die der vierten Spalte von Matrix b. 

Zusammen gefaßt bedeutet dies: 

Gm =Yj a n» b w (7-3) 

t = 1 

m = 1,..., r 
n = 1 ,..., s 
i = k\ 

r = Zeilenanzahl Matrix a 
k = Spaltenanzahl Matrix a 
i = Zeilenanzahl Matrix b 
s = Spaltenanzahl Matrix b 

Für die Herleitung von Regel (7.3) sei auf die Literatur [22] verwiesen. 

Genau wie eine Determinante läßt sich eine Matrix an der Hauptdiago¬ 
nalen spiegeln; dies nennt man hier auch transponieren. Aus einer (i * k)- 
Matrix entsteht eine (k * /)-Matrix; Zeilen und Spalten werden vertauscht. 
Per Definition heißt eine quadratische Matrix, deren Determinante den 
Wert 0 hat (vgl. Abschnitt 6.1), singulär, anderenfalls heißt sie regulär. 
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Für reguläre Matrizen gibt es neben den beschriebenen Verknüpfungen 
noch eine weitere Umformung, die des Öfteren in der Technik 
angewendet werden muß, das sogenannte Invertieren. 

Für die zu a inverse Matrix a~' gilt: 

a*a l =a l *a= E (7.4) 

E ist die Einheitsmatrix, sie enthält in der Hauptdiagonalen nur Einsen, 
alle anderen Glieder sind 0. Zur Aufstellung der inversen Matrix schreibt 
man zunächst eine gleichrangige Einheitsmatrix E neben die umzufor¬ 
mende Matrix a. Man betrachte dann beide Matrizen als Determinanten 
und wendet die elementaren Umformungen nach Kapitel 6 an, bis die 
Determinante der Matrix a eine Einheitsdeterminante darstellt. Zweck¬ 
mäßigerweise verwendet man hierzu den Gaußschen Algorithmus wie in 
Abschnitt 6.2, bis nur noch in der Hauptdiagonalen Werte stehen. Diese 
zieht man dann durch Division vor die Determinante, so daß die 
Einheitsdeterminante entsteht. Alle Rechenoperationen, die man an a 
ausführt, führt man parallel dazu auch an E aus. Wenn a zu E geworden 
ist, steht an der Stelle, an der die Einheitsmatrix stand, die Inverse a l . 
Prüfen kann man dies, indem man a und a~' miteinander multipliziert; 
das Ergebnis muß nach Gleichung (7.4) die Einheitsmatrix E sein. 


Beispiele 


a = 


a + b = 



b = 


16\ 

13 


a — b = 

11/ 


( 9 

15 21 

3 * a = I 

12 

6 27 

' 

^12 

9 24 
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Matrizenmultiplikation a*b 
r=k=i=s=3 

c u =I«iAi = 3*1+5*8 + 7*3 = 64 
c 'i2 = = 3*2 + 5*0 + 7*7 = 55 

c 13 = £a 1( £, 3 = 3*9 + 5*4 + 7*3 = 68 
c 21 = 'Z a 2 ,b, l =4*1 +2*8 + 9*3 = 47 
c 22 = Y, a 2tb,2 — 4*2 + 2*0 + 9*7 = 71 
c 23 = X a 2i*t3 = 4*9 + 2*4 + 9*3 = 71 
c 31 = £a 3 Ai = 4*1 + 3*8 + 8*3 = 52 
c 32 = £«3,6,2 = 4*2 + 3*0 + 8*7 = 64 
Gs = Z«3,^3 = 4*9 + 3*4 + 8*3 = 72 

bedeutet hier £ \ 


Die Ergebnismatrix lautet also: 


'64 

55 

68 

47 

71 

71 

.52 

64 

72 


Invertieren der Matrix a 


Da der Gaußsche Algorithmus bereits in Abschnitt 6.2 ausführlich 
behandelt wurde, sind hier nur noch die Ergebnisse der einzelnen Schritte 
angegeben. 



r 5 7 X 

I 0 -4,667 -0,333 I 
\0 -3,667 -1,333/ 


E = 


/! 0 0 
I 0 1 0 

VO 0 1. 

( 1 

-1,333 

\-1,333 


0 

1 

0 




Matrizen 


127 




E = 


Multipliziert man a mit a ', so erhält man 


-1,244 

1,244 

-0,311 

-0,286 

-0,786 

1 / 

-0,733 

-1,267 

2,067 

0,267 

-0,267 

0,067 

0,267 

0,733 

-0,933 


a*a 1 = 


'1,001 

0,002 

0,002 

0,001 

1,002 

0,002 

,0,001 

0,002 

1,002 


Die Abweichungen von der Einheitsmatrix sind auf Rundungsfehler 
zurückzuführen, da nur mit drei Dezimalstellen gerechnet wurde. 


Literaturhinweise 

G, W: [5], S. 537 bis 543; G,W: [9], S. 55 ff; G,W: [22], S. 236 ff; G,W: 
[25]; W: [1],; W: [16], S. 57 ff. 

Program mbesch reibung 

Das Programm Matrizenprogramme ist mit etwa 8K das längste in 
diesem Buch abgedruckte; Ursache sind die vielen und zum Teil 
komplizierten Rechenmöglichkeiten, die mit Matrizen möglich sind. Das 
Programm kann mit Matrizen rechnen, die aus bis zu zehn Zeilen bzw. 
Spalten bestehen; soll diese Zahl erhöht werden, müssen die dementspre¬ 
chenden Dimensionierungen geändert werden. 

Zeile 50 

bis 200 enthält wiederum das Hauptmenü mit Eingabe der Code¬ 
zahl und entsprechender Auswahl des Unterprogramms. 

Zeile 300 

bis 400 enthält unter Benutzung der 
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Zeilen 2100 

bis 2280 die Eingabe der ersten Matrix; die Gründe für diese 
Vorgehensweise wurden bereits in Kapitel 4 erläutert. 

Zeile 500 

bis 650 ist zuständig für Addition und Subtraktion, wobei wie 
übrigens bei allen anderen Eingaben von Matrizen vor 
Eingabe der zu addierenden bzw. subtrahierenden Matrix 
das Eingabefeld gelöscht und neu dimensioniert wird. Da 
nur Matrizen mit gleicher Spalten- und Zeilenzahl addiert 
und subtrahiert werden können, erfolgt die entsprechende 
Überprüfung in 

Zeile 530 

Durch die 


Zeilen 700 
bis 820 

Zeile 900 
bis 980 


Zeilen 1000 
bis 1250 

Zeile 1040 


Zeile 1220 


kann die Matrix mata, in der alle Ergebnisse gespeichert 
werden, ausgegeben werden. Mit dem Unterprogramm in 

wird die in mata befindliche Matrix mit einem Faktor 
multipliziert, wobei jedes Element mit diesem Faktor 
malgenommen wird, ln den 

ist dagegen die Multiplikation mit einer Matrix beschrieben. 

Entsprechend der Rechenvorschrift wird in 

geprüft, ob die Spaltenzahl von Matrix a mit der Zeilenzahl 

von Matrix b übereinstimmt; ist das nicht der Fall, erfolgt 

die entsprechende Ausgabe in 

und der Rücksprung zum Menü. 


Für die eigentliche Rechnung wird zunächst in den 

Zeilen 1060 

bis 1220 die Ergebnismatrix matc berechnet; diese wird in den 
Zeilen 1130 

bis 1190 nach mata gespeichert, damit sie mit dem vorher genannten 
Unterprogramm ausgegeben werden kann. Die 

Zeilen 1300 

bis 1420 transponieren mata, wobei die neue Matrix zunächst in 
matb gebildet wird und dann nach mata zurückgespeichert 
wird. In den 
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Zeilen 1500 


bis 

1860 

Zeilen 1520 

Zeile 

2370 

bis 

2430 

Zeile 

2500 

bis 

2670 


Zeile 

1550 

bis 

1740 

Zeile 

1750 

bis 

1790 

Zeile 

1810 


Zeilen 2000 

bis 

2050 

Zeile 

2300 

bis 

2670 


wird die Inversion der Matrix a durchgeführt. Hierzu wird 
in 

zunächst geprüft, ob es sich um eine quadratische Matrix 
handelt. Ist das der Fall, wird mata in 

zunächst in das Feld d gespeichert. Dies wird durchgeführt, 
weil anschließend zwischen 

die Eigenwertbestimmung von d durchgeführt wird. Hierzu 
wird das Programm aus Kapitel 6 benutzt, so daß die 
Variablennamen angepaßt werden müssen. 

Nur wenn der Eigenwert der Determinante von Matrix a 
ungleich 0 ist, wird in 

die inverse Matrix in matb berechnet und in 


nach mata zurückgespeichert. Ist der Eigenwert dagegen 
gleich 0, erfolgt in 

die Ausgabe, daß es sich um eine singuläre Matrix handelt 
und dementsprechend der Rücksprung ins Menü. 

Die Eigenwertbestimmung der Determinante einer Matrix 
kann über Codezahl 9 und dementsprechend über die 

und 


auch alleine angewählt werden. 

Bei diesem Programm handelt es sich, wie schon gesagt, um 
das gleiche Programm wie in Kapitel 6; für die Erläuterung 
des Ablaufes sei auf die Programmdokumentation dort 
verwiesen. 
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10 REM HP Matr1zenprogramme 
20 MODE 2 

30 DIM mata<10,10),matb(10,10) 

50 CLS : REM Hauptmenue 

60 PRINT’Was soll gemacht werden ? " : PRINT“ 


70 PRINT:PRINT“- Eingabe der 1.Matr1x";TAB<50);“<1>“ 

80 PRINT:PRINT"- Addition einer Matr1x“;TAB(50) i "<2) 

90 PRINT:PRINT“- Subtraktion einer Matrlx“;TAB(50) i * 
< 3) “ 

100 PRINT:PRINT‘- Multiplikation mit einem Faktor";T 
AB(50)J"C 4)" 

110 PRINT:PRINT“- Multiplikation mit einer Matrlx";T 
AB(50)J“<5)“ 

120 PRINT:PRINT“- Transponieren einer Matrlx";TAB<50 
);"( 6 )“ 

130 PRINT:PRINT‘- Invertieren einer Matrlx";TAB<50); 
"(7)" 

140 PRINT:PRINT"- Ausgabe der berechneten Matrlx'JTA 
B<50);"(8)“ 

150 PRINT:PRINT"- Eigenwertbestimmung einer Matrix"; 
TAB(50)»“<9) " 

160 PRINT:PRINT"- Schluss ";TAB<50);"C10)" 

170 INPUT was 

180 IF was<1 OR was)10 THEN 50 

190 ON was GOSUB 300,500,500,900,1000,1300,1500,700, 
2000,1900 
200 GOTO 50 
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300 REM HP/UP Eingabe der 1.Matrix 
310 ERASE matb : DIM matb(10,10) 

320 GOSUB 2100 : REM Eingabe Zellen- und Spaltenzahl 

330 GOSUB 2150 : REM Eingabe der Matrix 

340 FOR k=1 TO bzelle 

350 FOR 1=1 TO bspalte 

360 mata(k,1)=matb<k,1) 

370 NEXT 1 
380 NEXT k 

390 azelle = bzelle : aspalte = bspalte 
400 RETURN 

500 REM HP/UP + und - 

510 ERASE matb : DIM matb<10,10> 

520 GOSUB 2100 

530 IF bze 1 le<>aze 1 le OR bspalteOaspalte THEN 620 
540 GOSUB 2150 : REM Matrixeingabe 
550 FOR k=1 TO bzelle 
560 FOR 1=1 TO bspalte 

570 IF was=2 THEN mata(k,i>=mata(k,i) + matb(k,i) 

580 IF was = 3 THEN mata<k,1)=mata(k,1> - matb<k,l> 

590 NEXT 1 
600 NEXT k 
610 GOTO 650 

620 CLS : PRINT'Es koennen nur Matrixen mit gleicher 
Spalten- und Zeilenzahl addiert bzw. subtrahi 

ert werden.* 

630 PRINT:PRINT“Welter mit beliebiger Taste* 

640 IF INKEY* = * * THEN 640 
650 RETURN 

700 REM HP/UP Ausgabe 
710 CLS 
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720 a2=0 

730 FOR k=l TO aspalte 
740 FOR 1=1 TO azeIle 

750 LOCATE <k-a2*6>* 12-11,1 : PRINT mata(i.k) 

760 NEXT 1 

770 IF k=6 THEN GOSUB 800 
780 IF k=6 THEN a2=a2+l : CLS 
790 NEXT k 

800 LOCATE 1,23 : PRINT’Welter mit beliebiger Taste" 
810 IF INKEY*="" THEN 810 
820 RETURN 

900 REM HP/UP Faktormultipplikatlon 
910 CLS 

920 INPUT"Faktor : * i fak : CLS 

930 FOR k=1 TO azeIle 

940 FOR 1=1 TO aspalte 

950 mata(k,1)=mata<k,1)*fak 

960 NEXT 1 

970 NEXT k 

980 RETURN 

1000 REM HP/UP Multiplikation mit einer Matrix 
1010 ERASE matb : DIM matb(lO.lO) 

1020 GOSUB 2100 : REM eIngabe neue Matrix 
1030 GOSUB 2150 : REM dito 
1040 IF aspalteObze Ile THEN 1220 
1050 DIM matc(azeile,bspalte) 

1060 FOR k=l TO azeIle 
1070 FOR 1=1 TO bspalte 
1080 FOR n=l TO aspalte 

1090 matcCk,1) = mata(k,n> * matb(n,l) + matc(k,l) 
1100 NEXT n 
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1110 NEXT 1 
1120 NEXT k 

1130 ERASE mata : DIM mata(10,10) : REM Nullsetzen v 
on mata 

1140 FOR k=1 TO azeile 
1150 FOR 1=1 TO bspalte 
1160 mata(k,i) = matc<k,l) 

1170 NEXT 1 
1180 NEXT k 
1190 aspalte=bspalte 
1200 ERASE matc 
1210 GOTO 1250 

1220 CLS : PRINT'Matr1zenmultlpllkatIonen sind nur m 
oeglich, wenn die Spaltenanzahl von Matrix A = der Z 
ellenanzahl von Matrix B ist 

1230 L.0CATE 1,23 : PRINT"We 1 ter mit beliebiger Taste 

M 

1240 IF INKEY*="" THEN 1240 
1250 RETURN 

1300 REM HP/UP Transponieren 
1310 ERASE matb : DIM matb<10,10> 

1320 FOR k=1 TO azeile 
1330 FOR 1=1 TO aspalte 

1340 matb<l,k> = mata(l,k> : REM zunaechst Spelcheru 
ng nach mat b 
1350 NEXT 1 
1360 NEXT k 

1370 FOR k=1 TO azeile 
1380 FOR 1=1 TO aspalte 
1390 mata(l,k> = matb(k,l) 
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1400 NEXT i 
1410 NEXT k 
1420 RETURN 

1500 REM HP/UP Inversion einer Matrix 
1510 ERASE matb : DIM matb(10,10> 

1520 IF aspalte <> azelle THEN 1840 

1530 GOSUB 2300 : REM Eigenwertbestimmung mit Pruefu 
ng auf quadratische Matrix und Umspeicherung 
1540 IF dew = 0 THEN 1810 

1550 FOR k=l TO aspalte : REM Aufbau einer Elnheltsm 
atrix in matb 
1560 matb(k,k)=l 
1570 NEXT k 

1580 FOR m=l TO aspalte 
1590 FOR k=l TO aspalte 
1600 IF k=m THEN 1660 
1610 al= -mata<k,m)/mata(m,m) 

1620 FOR i=1 TO aspalte 

1630 mata<k,l> = mata<k,l> + mata(m,l)*al 
1640 matb<k,l) = matb(k,i> + matb(m,i)*al 
1650 NEXT 1 
1660 NEXT k 
1670 NEXT m 

1680 FOR k=l TO aspalte 
1690 al=mata<k,k) 

1700 FOR 1=1 TO aspalte 
1710 mata<k,1)=mata(k,1)/al 
1720 matb(k , 1)=matb(k,1)/a1 
1730 NEXT 1 
1740 NEXT k 

1750 FOR k=l TO aspalte 
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1760 FOR 1=1 TO aspalte 

1770 mata<k,1>=matb<k,1> : REM Umspeicherung der lnv 

ersen Matrix 

1780 NEXT 1 

1790 NEXT k 

1800 GOTO 1860 

1810 CLS : PRINT“Dle Matrix Ist slngulaer. Eine Inve 
rtlerung Ist sinnlos.“: LOCATE 1,23 : PRINT“Welter m 
lt beliebiger Taste“ 

1820 IF INKEY*="" THEN 1820 
1830 GOTO 1860 

1840 CLS : PRINT“Invertlerungen koennen nur von quad 
ratischen Matrizen durchgefuehrt werden": LOCATE 1,2 
3 : PRINT’We1ter mit bellebl 
ger Taste“ 

1850 IF INKEY* = ’" THEN 1850 
1860 RETURN 
1900 REM UP/HP ende 
1910 END 

2000 REM HP/UP Eigenwertbestimmung mit Pruefung auf 
quadratische Matrix und Umspeicherung 
2010 GOSUB 2300 : REM ab da passlert's 
2020 CLS : PRINT’Elgenwert : “Jdew 

2030 LOCATE 1,23 : PRINT’Welter mit beliebiger Taste 

M 

2040 IF INKEY*=““ THEN 2040 
2050 RETURN 

2100 REM UP Eingabe von bmat 
2110 CLS 

2120 INPUT'Wlev leie Zellen hat die Matrix“Jbze Ile 
2130 PRINT:INPUT’Wlev leie Spalten hat die Matrix “;b 
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spalte 
2140 RETURN 

2150 REM Eingabe der Matrix 
2160 WINDOW #1,1,80,8,18 
2170 FOR k=1 TO bzelle 
2180 CLS #1 

2190 LOCATE 1,6 : PRINT k;‘.te Zelle" 

2200 FOR 1=1 TO bspalte 

2210 PRINT#1,"matk»","»1»") : ";: INPUT#1,matb(k, 1 

> 

2220 NEXT 1 
2230 NEXT k 

2240 LOCATE 1,23 : PRINT’Alle Eingaben l.o. < /n> *; 
:INPUT a* 

2250 IF a* = " " THEN 2280 

2260 IF a*="n" THEN GOSUB 2100 ELSE 2240 
2270 GOSUB 2150 
2280 RETURN 

2300 REM Pruefung auf quadratische Matrix fuer Eigen 

Wertbestimmung ln 2500 

2310 IF aspalte = azeile THEN 2370 

2320 CLS 

2330 PRINT’ElgenwertbestImmungen sind nur von Determ 
lnanten quadratischer Matrizen moegllch“ 

2340 LOCATE 1,23 : PRINT’We1ter mit beliebiger Taste 

w 

2350 IF INKEY$="" THEN 2350 
2360 GOTO 50 

2370 REM Umspeicherungen fuer die Elgenwertbest1mmun 
gen ab 2500 

2380 DIM d<aspalte,aspalte) 
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2390 FOR k=l TO aspalte 
2400 FOR 1=1 TO aspalte 
2410 dCk,1)=mata(k, 1) 

2420 NEXT 1 
2430 NEXT k 
2440 danz=aspalte 

2450 GOSUB 2500 : REM Eigenwertbestimmung 
2460 ERASE d 
2470 RETURN 

2500 REM Eigenwertbestimmung von Determinanten quadr 
atischer Matrizen 

2510 REM Benutzt wird der Gauss'sche Algorithmus 

2520 FOR m=l TO danz 

2530 FOR k=1 TO danz 

2540 IF k=m THEN 2590 

2550 al= -d<k,m)/d<m,m> 

2560 FOR 1=1 TO danz 

2570 d<k,1)=d(k,1) + d<m,l)*al 

2580 NEXT 1 

2590 NEXT k 

2600 NEXT m 

2610 REM Loesungsdetermlnante Ist fertig 

2620 REM Eigenwertbestimmung durch Multiplikation de 

r Hauptdiagonalen 

2630 dew=l 

2640 FOR 1=1 TO danz 
2650 dew = dew * d(i,l) 

2660 NEXT 1 
2670 RETURN 
10000 CLS 

10010 FOR 1 = 1 TO aze Ile 
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10020 FOR k=1 TO azelle 

10030 LOCATE k*10-9,1 : PRINT matb<l,k> 

10040 NEXT k 

10050 NEXT 1 
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8.1 Lösung mit dem Gaußschen Algorithmus 

Grundlagen 

Eine Anzahl von n Gleichungen, in denen insgesamt m Unbekannte 

Vorkommen, heißen dann linear, wenn alle Unbekannten in der ersten 

Potenz auftreten. Lösen läßt sich dieses System, wenn 

□ die Anzahl der Gleichungen gleich der Anzahl der Unbekannten ist 
(m = ri) 

□ die Gleichungen linear unabhängig sind (zwei Gleichungen sind dann 
linear abhängig, wenn die eine durch Multiplikation mit einem Faktor 
in die andere überfuhrt werden kann) 

□ sich die Gleichungen nicht widersprechen 

Ein solches Gleichungssystem hat die Form 


«11*1 

+ «12*2 + • ' 

• + «1 m *m = h 

«21*1 

+ «22*2 + • • 

■ + «2m*m = b. 

«»1 *1 

+ «n2 *2 + • • 

■ + «nm *m = b, 


X,... x m = Unbekannte 

Indizes von a : 1. Index = 1... n = Zahl der Unbekannten 

2. Index = 1... m = Zahl der Gleichungen 

Wenn die Glieder b i ,b 2 ,..., 6„alle gleich 0 sind, heißt ein solches System 
homogen. In diesem Fall läßt es sich mit der vorliegenden Methode nicht 
lösen. 

Die Lösung des Gleichungssystems erfolgt hier mit Hilfe des 
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Gaußschen Algorithmus. Dabei wird vorausgesetzt, daß das System 
lösbar ist, eine Prüfung auf Lösbarkeit erfolgt nicht. 

Der prinzipielle Lösungsablauf ist wie in Abschnitt 6.2 bei der 
Bestimmung des Eigenwertes einer Determinante; da jedoch in der 
Ausführung einige Unterschiede bestehen, wird der Ablauf hier erneut in 
allen Einzelheiten erläutert. 

Die Lösung geht folgendermaßen vor sich: 

Die 1. Gleichung wird mit — a ki /a l , multipliziert und zur £-ten Gleichung 
addiert, wobei k von 2 bis n läuft (man beachte: 1 fehlt!). 

Das Ergebnis sieht dann folgendermaßen aus: 


«n*i +a l2 x 2 + .. 

■ + «1 m *m = />! 

«2 2*2 + •• 

■ + «2m*m = h 'i 

«n2 *2 + • • 

■ + a 'nm X m ~ K 


Wie man sieht, existiert nur noch in der 1. Gleichung ein Glied mit der 
Unbekannten x,. 

Im zweiten Schritt läßt man alle Glieder mit x 2 bis auf «22*2 
verschwinden. Hierzu nimmt man die 2. Gleichung mit — a' k2 /a 22 ma l ur| d 
addiert sie zur &-ten Gleichung (k = 1,3,..., n). Man beachte, daß 
diesmal bei k die 2 fehlt, es handelt sich um den 2. Schritt. 

Mit dem neuen Gleichungssystem verfährt man dann wieder entspre¬ 
chend, bis schließlich n Gleichungen vorliegen, die jeweils nur noch eine 
Unbekannte enthalten. Diese lassen sich dann leicht lösen. 


Gleichung I: 
Gleichung II: 
Gleichung III: 

1. Schritt 


Beispiel 

2 *! + * 2 - * 3 = - 2 

x, — 3 x 2 — 2x 3 — —5 
— x, — 4 x 2 + x 3 = 7 


«n «n 2 

Gleichung I multipliziert mit —1/2 und zu Gleichung 11 addiert, ergibt als 
neue Gleichung II': 
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Gleichung II': 


■ 3,5x 2 — 1,5*3 = 


k = 3: 


^31 


= 1/2 


Jetzt wiederum Gleichung I multipliziert mit 1/2 und addiert zu 
Gleichung III ergibt als Gleichung III': 

Gleichung III': —3,5*+ 0,5* = 6 

Als neues Gleichungssystem ergibt sich somit: 

Gleichung I: 2*, + * 2 — * 3 = — 2 

Gleichung II': —3,5* 2 — 1,5* 3 = —4 

Gleichung III': — 3,5* 2 + 0,5* 3 = 6 


2. Schritt: 


k = 1: = 

a,, a. 


1 


1 

3,5 


*22 “22 2,5 

Gleichung II' multipliziert mit 1/3,5 und addiert zu Gleichung I ergibt: 


Gleichung I': 


-3 5 11 

2, '-0*- = -w 


k= 3: = 


-3,5 

-3,5 


= -1 


Gleichung III": 


2*3 = 10 


Damit erhält man nach dem zweiten Schritt als neues Gleichungssystem: 

5 11 

Gleichung I': 2*, — — * 3 


Gleichung II': 
Gleichung III": 


3,5 V3 3,5 
— 3,5* 2 — 1 , 5*3 = —4 
2*3 = 10 


An dieser Stelle würde man normalerweise mit dem Einsetzverfahren 
weitermachen; da das Programm jedoch den Gaußschen Algorithmus bis 
zum Ende durchführt, sei dies auch hier zur Verdeutlichung gemacht. 
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3. Schritt 


k = 1: 

Gleichung I": 

k = 2: 


Gleichung II": 


a k3 _ _ 4>3 . 
«33 «33 

2x , =4 

«> i - 3 _ _ «23 
«33 _ «33 


5 


3,5 
2 “ 


= 5/7 


-1,5 

2 


= 3/4 


-3,5*2 = 3,5 


Das zusammengefaßte, endgültige Gleichungssystem lautet also: 


Gleichung I": 2.x:, =4 

Gleichungll": — 3,5* 2 = 3,5 

Gleichung III": lx 3 = 10 

Daraus folgen die einzelnen Unbekannten als: 

= 2 
x 2 = -1 
*3 = 5 

Setzt man diese Werte zur Probe in Gleichung I ein, so erhält man: 


2jc, + x 2 - x 3 = -2 
2*2 + (-1) - 5 = -2 

— 2=—2 w.d.R.d.R.b.w. 


(womit die Richtigkeit der Rechnung bewiesen wäre) 


Literaturhinweise 

G, W: [5], S. 116 ff; G,W: [9], S. 81 ff; G,W: [13], S. 281 ff; G,W: [16], 
S. 87 ff; G, W: [19], S. 37 ff; G,W: [22], S. 84 ff. 
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Programmbeschreibung 

Bei der Eingabe des Gleichungssystems ist es wichtig, daß die 
Gleichungen in der Form vorliegen, wie sie unter Grundlagen angegeben 
ist; das Restglied, d.h. der Faktor ohne Unbekannte, steht rechts vom 
Gleichheitszeichen. Falls eine Unbekannte in einer Gleichung nicht 
vorkommt, muß an der entsprechenden Stelle der Faktor 0 angegeben 
werden. Im übrigen läuft das Programm wie folgt ab: 

Zeile 20 Umschaltung auf 80 Zeichen pro Zeile 
Zeilen 30 

bis 260 Eingabe des Gleichungssystems Besonderheit hier: 

Zeile 40 Das Feld wird für die einzelnen Faktoren erst dimensioniert, 
nachdem die Anzahl der Unbekannten eingegeben worden 
ist. Damit geht einerseits nicht zuviel Speicherplatz verloren, 
andererseits lassen sich aber Gleichungen mit beliebig vielen 
Unbekannten lösen. 

Zeile 170 Die vorherige Zeile wird gelöscht, der übrige Text bleibt 
stehen. 

Zeilen 280 

bis 370 Lösung des Gleichungssystems 

Am Ende ist in jeder Gleichung neben dem Restglied nur noch 
ein Glied vorhanden, und zwar in der 1. Gleichung das 
1. Glied, in der 2. Gleichung das 2. Glied usw. 

Zeile 420 Berechnung der jeweiligen Unbekannten durch Division des 
Restgliedes durch das jeweilige noch vorhandene Glied 
Benötigter Speicherplatz: etwa 1,6 K. 
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10 REM Programm 1 ln 2.0 
20 MODE 2 

30 PRINT "Anzahl der Unbekannten :";:INPUT n 
40 DIM a<n,n+l),x<n) 

50 CLS 

60 PRINT "Geben Sie die Faktoren der einzelnen Gl 
elchungen ein." 

70 PRINT:PRINT"Geben Sie alle Werte einer Gleichung 
einschliesslich des Restglledes ein, bevor Sie zur n 
aechsten uebergehen." 

80 PRINT:PRINT"Beruecksichtigen Sie fehlende Unbekan 
nte ln einer Gleichung durch Eingabe von 0." 

90 LOCATE 1,25:PRINT“Weiter mit beliebiger Taste" 

100 IF INKEY*="" THEN 100 
110 LOCATE 1,25:PRINT SPC<27) 

120 LOCATE 1,12 : PRINT"A 1,3 = 1.Gleichung, 3.Glied 
etc" 

130 PRINT:PRINT"Restglled = rechts vom Gleich 
heitszeichen !!!!!“ 

140 FOR 1 = 1 TO n 
150 FOR k = 1 TO n 

160 LOCATE 5,18 : PRINT "A" i 1J","iki:INPUT a<l,k> 

170 LOCATE 5,18 : PRINT SPC<30> 

180 NEXT k 

190 LOCATE 5,18 : PRINT "B";1»:INPUT a(l,n+l) 

200 LOCATE 5,18 : PRINT SPC<30) 

210 NEXT 1 

220 PRINT:PRINT:PRINT "War die Eingabe l.o. ( /n> 

:INPUT a$ 

230 IF a* = "" THEN 260 
240 IF a$ = "n" THEN 50 
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250 GOTO 220 
260 CLS 

270 PRINT’Ich arbeite. Bitte einen Augenblick Geduld 

M 

280 FOR k = 1 TO n 
290 FOR 1 = 1 TO n 
300 IF i=k THEN GOTO 360 
310 f = -a<l,k)/a<k,k> 

320 FOR 1=1 TO n+1 

330 b=a<k,1) * f 

340 a<l,l> = a< 1,1) + b 

350 NEXT 1 

360 NEXT 1 

370 NEXT k 

380 CLS 

390 PRINT “We've got lt ! ! ! ■ 

400 PRINT 

410 FOR 1=1 TO n 

420 x<l> = a(1,n+1) /a( 1 , 1) 

430 print "x*;i;* = ";x(i> 

440 NEXT 1 

450 LOCATE 1,23 : PRINT’Soll noch ein weiteres Gleic 
hungssystem geloest werden ( /n) “;:INPUT a$ 

460 IF a* = ” THEN RUN 20 
470 IF a$ = "n" THEN 490 
480 GOTO 450 
490 MODE 1 

500 LOCATE 3,10 : PRINT’Schoenen Feierabend" : LOCAT 
E 1,24 
510 END 
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8.2 Determinantenverfahren 

Im Gegensatz zu dem vorstehend beschriebenen Programm erfolgt hier 
eine Prüfung auf Lösbarkeit des Systems. Dazu wird zunächst die 
Koeffizientendeterminante D des Gleichungssystems berechnet, d.h. die 
Determinante, deren einzelne Glieder aus den Faktoren a,, bis a nm 
bestehen. Lösbar ist das System, wenn /)<) 0 ist. 

Daß D zu 0 wird, kann verschiedene Ursachen haben; diese werden 
nicht näher untersucht. 

Die eigentliche Lösung erfolgt jetzt nach der Cramerschen Regel. 
Hierzu wird nacheinander jeweils eine Spalte der Koeffizientendetermi¬ 
nante ersetzt durch die Faktoren h x bis h n und dann diese neue 
Determinante durch D geteilt; der Quotient stellt die jeweilige Glei¬ 
chungslösung dar. 


Beispiel 

Im Folgenden wird das in Abschnitt 8.1 angegebene Beispiel erneut 
durchgerechnet. Für die Berechnung des Eigenwertes der Determinanten 
sei auf Kapitel 6 verwiesen. 


Gleichung I: 


2*i 

+ *2 

Gleichung II: 


*i 

- 3.y 2 

Gleichung III 


-*i 

- 4*2 


2 

1 

-1 


D = 

1 

-3 

-2 

= —: 


-1 

-4 

1 



-2 

1 

-1 


Z), = 

-5 

-3 

-2 

= 28 


7 

-4 

1 



2 

-2 

-1 


d 2 = 

1 

-5 

-2 

= 14 


-1 

7 

1 



■ *3= -2 

-2.v 3 = -5 

+ *3=7 


14< > 0 es gibt also eine 

eindeutige Lösung 

— 28 

x>=DJD=- i- = 2 


14 

*2 = D 2 /D = — = - 1 
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2 1 -2 

D 3 = 1-3 -5 =-70 x 3 = D i /D=- 1 4 ° = 5 

-1-4 7 

Beide Verfahren kommen also zu der gleichen Lösung. 

Literaturhinweise 

G: [5], S. 264 bis 267; G: [6], S. 76ff; G: [8], S.134 bis 137; G: [25], 
S. 58 ff; G, W: [1 ], S. 75 ff; G,W:[9], S. 71 bis 81; G,W:[22], S.72ff. 

Programmbeschreibung 

Für die Eingabe des Gleichungssystems gilt das bereits in 
Abschnitt 8.1 Gesagte. 

Zeile 20 

bzw. 

Zeile 670 

bis 910 Eingabe des Gleichungssystems. 

Anzahl der Gleichungen = Uzahl 
Abspeicherung im Feld gs 

Zeile 40 Dimensionierung der Determinante d , deren Glieder im 
Laufe der nachfolgenden Berechnungen unterschiedliche 
Werte annehmen. 

Zeile 50 

bis 90 Zunächst werden die Faktoren der Koeffizientendetermi- 

nante nach d gespeichert. 

Zeile 100 

bzw. 

Zeile 520 

bis 660 Eigenwertbestimmung von d. 

Der Eigenwert liegt anschließend in der Variablen dew vor. 
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Zeile 120 

bzw. 

Zeile 320 

bis 510 Eigenes Ausgabeprogramm für den Fall, daß die Koeffi¬ 
zientendeterminante gleich 0 ist, da das Gleichungssystem 
nicht lösbar ist. Die Werte der einzelnen Zählerdeterminan¬ 
ten werden ausgegeben. 

Zeile 130 

bzw. 

Zeile 930 

bis 1080 Berechnung der Zählerdeterminante für den «Normalfall». 

Abspeicherung im Feld x(i). 

Zeilen 150 

bis 170 Berechnung der einzelnen Unbekannten 

Zeilen 180 

bis 230 Lösungsausgabe 

Benötigter Speicherplatz: etwa 3 K. 


10 REM Programm 1 ln det 

20 GOSUB 670 : REM Eingabe des Glelchungssystems 
30 danz = Uzahl 
40 DIM d(danz,danz) 

50 FOR k = 1 TO danz : REM Umspeicherung der Koefflz 
lentendetermlnante 
60 FOR 1 = 1 TO danz 
70 d(k,1) = gs<k,1) 

80 NEXT 1 
90 NEXT k 

100 GOSUB 520 : REM Eigenwertbestimmung 
110 kd = dew 
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120 IF kd = 0 THEN 320 : REM eigenes Ausgabeprogramm 
130 GOSUB 930 : REM Berechnung der einzelnen Zaehler 
determInanten 

140 REM Berechnung der einzelnen Zaehlerdetermlnante 
n 

150 FOR 1 = 1 TO uzahl 
160 x<l> = x<l> / kd 
170 NEXT 1 

180 REM Loesungsausgabe 
190 CL3 

200 PRINT’So,Jetzt haben wlr’s geschafft!" : PRINT 
210 FOR 1 = 1 TO uzahl 
220 PRiNT"x";i;" = ";x<i> 

230 NEXT 1 

240 PRINT:PRINT:PRINT"So 11 noch ein Gleichungssystem 
geloest werden < /n> ";:INPUT a* 

250 IF a* = "" THEN RUN 
260 IF a* = "n" THEN 280 
270 GOTO 240 

280 CLS : LOCATE 15,10 : PRINT'Na gut. Wer nicht wll 
1, der hat schon." 

290 LOCATE 30,13 : PRINT"Tschuess" 

300 LOCATE 1,23 
310 END 

320 REM Ausgabe fuer Koeffizientendeterminante = 0 
330 CLS 

340 PRINT'Koefflzlentendetermlnante = 0 " 

350 PRINT 

360 PRINT"Moegllche Ursachen" 

370 PRINT" - lineare Abhaenglgke1t von mindestens 2 
Glelchungen" 
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380 PRINT “ (wenn mindestens eine Zaehlerdetermlnan 
te = 0)* 

390 PRINT 

400 PRINT " - Widerspruch lm Glelchungssystem" 

410 PRINT " (wenn alle Zaehlerdetermlnanten <> 0)" 
420 PRINT 

430 PRINT" - homogenes Gleichungssystem" 

440 GOSÜB 930 : REM Berechnung der einzelnen Zaehler 

determlnanten 

450 PRINT:PRINT:PRINT 

460 PRINT "Werte der einzelnen Zaehlerdetermlnanten’ 
470 PRINT 

480 FOR 1 = 1 TO uzahl 
490 PRINT x(1) 

500 NEXT 1 
510 GOTO 240 

520 REM UP Eigenwertbestimmung von Determinanten 

530 FOR m = 1 TO danz 

540 FOR k = 1 TO danz 

550 IF k = m THEN GOTO 600 

560 al = -d(k,m)/d(m,ra) 

570 FOR 1 = 1 TO danz 

580 d(k,1) = d(k, 1) + d(m,l)*al 

590 NEXT 1 

600 NEXT k 

610 NEXT m 

620 dew = 1 

630 FOR 1 = 1 TO danz 

640 dew = dew*d(1,1) 

650 NEXT 1 
660 RETURN 



Lineare Gleichungssysteme 


151 


670 REM Eingabe des Glelchungssystems 
680 CLS 

690 PRINT"Wlevleie Unbekannte hat das Gleichungssyst 
em " > 

700 INPUT uzahl 

710 DIM gs(uzahl,uzahl+l>,x(uzahl> 

720 PRINT 

730 PRINT"A3 = Faktor vor x3, A2 = Faktor vor x2 etc 

740 PRINT 

750 PRINT"fehlende Glieder durch 0 beruecksichtigen 
! ! ! " 

760 FOR k = 1 TO uzahl 

770 LOCATE 1,10 : PRINT k;"te Gleichung" 

780 FOR 1 = 1 TO uzahl 

790 LOCATE 1,12 : PRINT"A’;1;" : "J 

800 INPUT qs(k,1) 

810 LOCATE 1,12 : PRINT" 

820 NEXT 1 

830 LOCATE 1,12 :PRINT "Y : "5 
840 INPUT gs(k,uzahl+l) 

850 NEXT k 
860 PRINT 

870 PRINT:PRINT " Eingabe o.k. ( /n> "5: INPUT a* 

880 IF a* = "" THEN 910 
890 IF a* = ’n" THEN RUN 
900 GOTO 870 
910 PRINT 
920 RETURN 

930 REM UP zaehlerdetermInantenberechnung 



152 


Lineare Gleichungssysteme 


940 FOR n = 1 TO uzahl 
950 FOR k = 1 TO danz 
960 IF k = n THEN GOTO 1000 
970 FOR 1 = 1 TO danz 
980 LET d(1,k> = gs (1,k) 

990 NEXT 1 
1000 NEXT k 

1010 FOR 1 = 1 TO danz 
1020 d<l,n> = gs(l,danz+l> 
1030 NEXT 1 
1040 GOSUB 530 

1050 REM Eigenwertbestimmung 
1060 x(n) = dew 
1070 NEXT n 
1080 RETURN 
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Wenn der Wert einer Größe gemessen wird, ist dieser Meßwert 
normalerweise mit Fehlern behaftet. Man unterscheidet hier systemati¬ 
sche Fehler wie Fehler der Analysenmethode usw. von zufälligen Fehlern. 
Systematische Fehler sind grundsätzlich feststellbar und können demzu¬ 
folge berücksichtigt werden, zufällige Fehler können jedoch nur mit Hilfe 
der Wahrscheinlichkeitsrechnung erfaßt und korrigiert werden. 

Soll eine Größe jc bestimmt werden, wird man normalerweise die 
Einzelmessung «-mal wiederholen. Der wahrscheinliche Wert der Größe 
ist dann der Mittelwert. 


- _ *1 + X 2 + ■ ■ • + *n 

n 


~ I 

»i-i 


(9.1) 


x ist ein Schätzwert für den wahren Wert von x. Daß für die Bildung von x 
das arithmetische Mittel und nicht z. B. das geometrische herangezogen 
werden muß, kann über die Methode der minimalen Fehlerquadrate 
nachgewiesen werden, die von Gauß zum erstenmal zur Bahnberechnung 
von Planeten angewendet wurde. 

v t bezeichne die jeweilige Abweichung der Meßwerte x t vom tatsächli¬ 
chen Wert x. Die Forderung nach Gauß ist dann, daß die Summe der 
Fehlerquadrate zu einem Minimum werden soll: 


£ vf = Minimum! 

i = 1 

bzw. in der Schreibweise nach Gauß 


[r>r>] = Minimum! 

Der quadratische Fehler wird genommen, damit sich positive und 
negative Fehler nicht aufheben; außerdem werden «Ausreißer» verstärkt 
berücksichtigt. 
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Ersetzt man v f durch x — x h so ergibt sich: 

M = (x- x,) 2 + (x - x 2 ) 2 + ... + (x - x n ) 2 

Daß diese Summe zu einem Minimum werden soll, heißt, daß die erste 
Ableitung zu 0 werden muß (vgl. Kapitel 2): 

= 2(x — x,) + 2(x — x 2 ) + • • • + 2(x — x„) = 0 

x - x, + x - x 2 + ... + x - x n = 0 
X, + x 2 + ... + x„ = nx 

Xi Xj “1" -Y„ _ . . . 

—---- = x = arithm. Mittel 

n 

Es reicht jedoch häufig nicht aus, den Mittelwert x zu kennen, sondern 
man muß ein Maß dafür haben, wie wahrscheinlich dieser Wert ist bzw. in 
welcher Bandbreite der wahre Wert x zu suchen ist. Hierfür kennt man 
zunächst die Varianz s 2 , die folgendermaßen definiert ist: 

s 2 = ——r X (Xj — x) 2 (9.2) 

n — t i = i 

Die positive Wurzel aus s 2 ist die Standardabweichung s, auch mittlerer 
Fehler der Einzelmessung oder empirische Streuung genannt. 


Wenn der Meßwert angegeben wird, lautet die Angabe dann: 

x± s 

wobei s sowohl in Prozent als auch in absoluten Einheiten angegeben 
werden kann. 

Da die Einzelmessungen mit Fehlern behaftet sind, muß auch der 
Mittelwert mit Fehlern behaftet sein. Ein Maß hierfür ist die Streuung des 
Mittelwertes; dieser Wert ist auch als Standardabweichung des Mittel¬ 
wertes bekannt: 
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Beispiel 

Die Masse eines Produktgefaßes wurde 8mal überprüft: 


1 

1 

2 

3 

4 

5 

6 

7 

8 

* [kg] 

20,34 

20,34 

20,38 

20,20 

20,31 

20,50 

20,33 

20,34 


Wie groß ist der Mittelwert und die Standardabweichung? 

_ 20,34 + 20,34 + 20,38 + 20,20 + 20,31 + 20,50 + 20,33 + 20,39 

* = - 8 - 

= 20,35 kg 


i 

i 

2 

3 

4 

5 

6 

7 

8 

x — x ( 

0,01 

0,01 

-0,03 

0,15 

0,04 

-0,15 

0,02 

-0.04 


/0,01 2 + 0,01 2 + 0,03 2 + 0,15 2 + 0,04 2 + 0,15 2 + 0,02 2 + 0,04 2 
V 8-1 

5 = 0,08 = 0,41 % 
x = 20,35 kg ± 0,41 % 


Literaturhinweise 

G, W: [2], S. 41 bis 49; G,W: [8]; G,W: [10], S. 188 bis 191; G,W: [14]; 
G, W: [17], S. 100 ff. 


Programmbeschreihung 

Das Programm Fehler kann bis zu 100 Meßwerte verarbei¬ 
ten; ist die Anzahl höher, muß x in 
Zeile 20 entsprechend höher dimensioniert werden. In den 
Zeilen 50 

bis 110 werden die Meßwerte eingegeben und die Anzahl festgestellt; 
beendet wird der Vorgang durch Eingabe von 333.33. In 

Zeile 130 

bis 180 wird der Mittelwert berechnet, in 
Zeile 190 
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bis 240 die empirische Varianz und in 

Zeile 260 als Wurzel aus der Varianz die empirische Streuung. 
Schließlich wird in 

Zeile 270 

bis 290 noch die Standardabweichung berechnet, und zwar sowohl in 
absoluten Einheiten als auch in %. 

Über die 

Zeile 300 

bis 360 werden alle berechneten Werte schließlich ausgegeben. 
Programmlänge: etwa 1,7 K. 


10 REM Programm Fehler 
20 DIM x<100) 

30 REM Dateneingabe 
40 CLS 

50 PRINT'Geben Sie die elnzelen Messwerte ein. Beend 
en Sie die Eingabe durch 333.33." 

60 WINDOW #1,1,80,3,25 
70 FOR 1=1 TO 100 

80 PRINT#1 ,"x"> 1 5": INPUT#1,x(1) 

90 IF x<1>=333.33 THEN 110 

100 NEXT 1 

110 anzahl = 1-1 

120 CLS 

130 REM Mittelwert 
140 sum=0 

150 FOR 1=1 TO anzahl 
160 sum=sum + x(l) 

170 NEXT 1 

180 mw=sum/anzahl 

190 REM Empirische Varianz 

200 sf2=0 
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210 FOR 1=1 TO anzahl 

220 sf2 = sf2 + (mw-x(1)>*(mw-x( 1) > 

230 NEXT 1 

240 var=sf2/(anzahl-1) 

250 REM Empirische Streuung 
260 str=SQR(var> 

270 REM Standardabweichungen 
280 stabw=SQR(sf2/(anzahl*(anzahl-1 
290 stabwrel=stabw/ABS(mw>*100 
300 REM Ausgabe 
310 CLS 

320 PRINT“Mittelwert : 

330 PRINT:PRINT“Standardabwelchung 
340 PRINT“ +-’ 

350 PRINT:PRINT"EmpIr. Streuung 
360 PRINT:PRINT’Empir. Varianz 
370 LOCATE 1,23 
380 END 


))) 


“; mw 

: +-“>stabw 

;stabwre1;" %“ 
: “Jstr 


; var 




10 

Ausgleichsrechnung 


10.1 Lineare Regression 

Grundlagen 

Ein oft auftauchendes Problem in nahezu allen Bereichen der Wirtschaft 
besteht darin, zu einer gegebenen Reihe von Meßwertpaaren den 
funktionalen Zusammenhang zu finden, der die Werte optimal be¬ 
schreibt, gleichzeitig aber Meßfehler möglichst eliminiert. Normalerweise 
geht man so vor, daß man einen bestimmten Funktionstyp auswählt, 
anhand der Meßwerte für diesen Typ entsprechende Parameter berechnet 
und dann schließlich mittels eines speziellen Kriteriums die Güte der 
Annäherung beurteilt. Führt man dies für verschiedene Funktionstypen 
durch, kann man anhand der Gütekriterien feststellen, welcher Funk¬ 
tionstyp am besten geeignet ist. 

Die einfachste Funktion, die für diese Ausgleichsrechnungen betrachtet 
wird, ist die Gerade. 

y = a 0 + a l x (10.1) 

In diesem Fall spricht man von linearer Regression. Hier geht es darum, 
die Faktoren a 0 unda, aus Gleichung (10.1) zu bestimmen. Hierfür gibt es 
diverse Methoden; am bekanntesten ist die Approximation durch 
orthogonale Funktionssysteme, vorgeschlagen von Fourier und Tscheby- 
scheff [3] (Seiten 526 bis 533) und die Methode der kleinsten 
Fehlerquadrate nach Gauß. Im Folgenden wird wie auch in Kapitel 9 das 
Gaußsche Verfahren verwendet. 

Die zentrale Forderung nach Gauß lautet, daß die Summe der 
Fehlerquadrate zu einem Minimum werden soll (vgl. Kapitel 9). Als 
Fehler v : ist hierbei jeweils der Abstand des Meßwertes y t zu dem 
berechneten Funktionswert/( jc,) zu verstehen, die beide zu dem Meßwert 
Xj gehören. 
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Z CF« -/(*.)) 2 = M = Minimum 

i = 1 

M = (öi*i +öo-Ji ) 2 + ( a 1 x 2 + a 0 - y 2 ) 2 + ... 

+ (a l x l +a 0 -y l ) 2 (10.2) 

Aus dieser Gleichung müssen «, und a Q so bestimmt werden, daß [m>] zu 
einem Minimum wird. Die 1. Ableitung muß also zu 0 werden; da a l und 
a 0 beides Veränderliche sind, muß nach diesen beiden Größen partiell 
abgeleitet werden. x f und y t sind als Konstanten zu betrachten (für den 
Begriff der partiellen Ableitung und für die verwendeten Ableitungsregeln 
siehe z.B. [22], S. 252 ff). 

= 2(a, -V, + a 0 -y x ) + 2 («, x 2 + a 0 - y 2 ) + ... 

+ 2(a 1 x„ + a o ->’ 1 ) = 0 (10.3) 

r/[yr] 

= 2 ( 0 !.v, +a 0 -y l )x l + 2(a, x 2 + a 0 -y 2 )x 2 

+ ... + 2(a,.v„ + a 0 ->'„).v„ = 0 (10.4) 

Aus Gleichung (10.3) folgt nach Umformung die Berechnungsgleichung 
für a 0 

öi(*i +x 2 + .. ■ + x,) + n*a 0 = y t + y 2 + ... + y t 


1 " n " 
n i = 1 “i=l 


(10.5) 


Im Folgenden werden aus Vereinfachungsgründen die Laufwerte für die 
Summationen weggelassen; alle Summen laufen von 1 bis n. 

Aus Gleichung (10.4) folgt nach Umformung: 

a x xj + a 0 x x - y, x, + a { x\ + a 0 x 2 - y 2 x 2 
+ ... + a t x 2 +a 0 x n -y„x„ =0 

«1 Z-v 2 + «oZ*i = Z*.* 

Wenn man in diese Gleichung die Gleichung (10.5) einsetzt, ergibt sich die 
Berechnungsformel für a,; 
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«, ( e «. ! - ( -^) = Zw, - 


Zw 


Z*iZyi 


Cl, = 


Z*?~ 


G >,) 2 


( 10 . 6 ) 


Als Gütekriterium wird häufig das Bestimmtheitsmaß R 2 verwendet [11 ], 
(S. 579 bis 591). Zu dessen Berechnung wird zunächst der Mittelwert y der 
j'-Werte gebildet und dann die Summe der quadratischen Abweichungen 
der Meßwerte y t vom Mittelwert und der berechneten Funktionswerte 
f(xj) vom Mittelwert ermittelt. 

SQMY = Z(f( x i) ~ y ) 2 (10.7) 

SQAY = X(y-y) 2 (10.8) 

Das Bestimmtheitsmaß ist definiert als der Quotient dieser beiden 
Summen. 


R 2 


SQMY 

SQAY 


(10.9) 


Für den Idealfall, daß die Ausgleichsgerade durch alle Meßpunkte geht, 
sind beide Summen gleich, R 2 ist demnach gleich 1. Je mehr R 2 von 1 
verschieden ist, desto schlechter ist die Annäherung. 

Eine weitere Möglichkeit für ein Gütekriterium ist die Bestimmung der 
Abweichung dy. Hier wird die quadratische Summe der Differenzen 
zwischen Meßwert^,- und Funktionswert f(x,) gebildet, diese Summe auf 
die Anzahl der Meßwerte bezogen und hieraus die Wurzel bestimmt. 


dy = 


pZcM 

n 


( 10 . 10 ) 


Wird dieser Wert auf den Mittelwert bezogen, erhält man die relative 
Abweichung dyrel 
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dyrel = 


^* 100 % 


( 10 . 11 ) 


dy und dyrel sind also mit den in Kapitel 9 besprochenen Standardabwei¬ 
chungen verwandt. 


Beispiel 

Für folgende Wertetabelle ist die optimale Gerade zu bestimmen. Des 
weiteren sind die Gütekriterien zu berechnen. 


I*. 


i 

l 

2 

3 

4 

5 

6 

X 

0,5 

1 

2 

3 

5 

7 

y 

0,1 

0,3 

0,7 

1,5 

2 

4 


18,5 X.v ; = 8,6 = 44,25 £*? 




44,25 - 


88,25 



0,568 


88,25 


Hiermit folgt aus Gleichung (10.5): 


a 0 = 


8,6 

~ 6 ~ 


0,568 

6 


• 18,5 = -0,318 


Die Gleichung der optimalen Geraden lautet also 


y = -0,318 + 0,568* 
Mittelwert: y = 8,6/6 = 1,433 
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1 

i 

2 

3 

4 

5 

6 

/(*,') 

-0,035 

0,250 

0,818 

1,386 

2,522 

3,659 

/(*«) - y 

-1,468 

-1,183 

-0,615 

-0,047 

1,089 

2,226 

yt-y 

-1,333 

-1,133 

-0,733 

0,067 

0,567 

2,567 

y t -/(*, ) 

0,135 

0,050 

-0,118 

0,114 

-0,522 

-0,341 


SQMY = ^(Rx i -y) 2 = 10,08 
sQAY = Z(yi-y) 2 = 10,51 
R 2 = 10,08/10,51 = 0,959 
IO ’i-A*i)) 2 = 0,439 



dyrel = 


0,0266 

1,433 


100% = 18,9% 


Das Programm zur linearen Regression ist unter Abschnitt 10.2 zu finden, 
da es mit den dortigen Programmen verknüpft wurde. 


Literaturhinweise 


G,W: [14], S. 80ff; G,W:[18]; G,W: [22], S. 340 bis 343. 
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10.2 Anpassung durch einfache Funktionen 


Grundlagen 

Wenn die Meßwertanpassung durch eine Gerade kein befriedigendes 
Ergebnis liefert, wird man normalerweise zunächst auf andere einfache 
Funktionen zurückgreifen; hier sind insbesondere die Logarithmusfunk¬ 
tion, die e-Funktion und die Potenzfunktion zu nennen. Man geht dann 
so vor, daß man die jeweilige Funktion durch Termersetzung in eine 
Geradengleichung überführt und die zu berechnenden Koeffizienten über 
Gleichung (10.5) und (10.6) ermittelt. Da durch die Termersetzung neue 
Variablen eingeführt werden, lautet die Geradengleichung 


p = a 0 + a x z 


( 10 . 12 ) 


wenn in der Ausgangsgleichung die Variablen .v und y heißen sollen. 
Die Faktoren a 0 und a x berechnen sich also nach 


Yj Z iPi “ 


I-.!/>, 


L-t 2 




(10.13) 


_ oiZfi 


(10.14) 


Die Gleichungen für die Gütekriterien (10.7) bis (10.11) können dagegen 
unverändert übernommen werden. 


Logarithmusfunktion 

Funktionsgleichung: y = a 0 + a { ln(.r) (10.15) 

Ersetzt man y durch p und ln (x) durch z, ergibt sich die Geradengleichung 
(10.12), so daß in Gleichung (10.13) und (10.14) eingesetzt werden kann. 
Beim Einsetzen wird die Termersetzung wieder rückgängig gemacht, d. h., 
in beiden Gleichungen wird p durch y und z durch ln (x) ersetzt. 
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Somit ergibt sich: 


5>,mW - 


(10.16) 


(10.17) 


Exponentialfunktion 
Funktionsgleichung: y = a 0 e“>* 

Diese Gleichung wird zunächst beidseitig logarithmiert: 

ln(j) = ln(a 0 ) + a, x 


(10.18) 


(10.19) 


Hier wird jetzt ln (y) durch p und x durch z ersetzt. Beim Einsetzen in 
Gleichung (10.14) muß statt a 0 zunächst ln(a 0 ) geschrieben werden, da 
Gleichung (10.19) sonst keine Geradengleichung darstellt. Beim Einset¬ 
zen wird die Termersetzung wieder rückgängig gemacht: 


E^lntVi) - ■ 

I >, 2 - 


£ x i I |n (>’■■) 


, n( „ o) ,IüM_£Ja bzw , 

n n 


a 0 = exp 


Z ln (Ti) q.I-Vj 


( 10 . 20 ) 


( 10 . 21 ) 


Funktionsgleichung: 


Potenzfunktion 
y = a 0 x a ' 


( 10 . 22 ) 


Logarithmieren führt zu 


ln(>0 = ln(a 0 ) + öi ln(x) 


( 10 . 23 ) 
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Termersetzung: ln(_y)=/> ln(x) = z 

Auch hier muß in Gleichung (10.14) zunächst ln(a 0 ) statt a 0 geschrieben 
werden. 


I,ln W ) <■»<*>) 


I*? 


l-v,r 


l„, ao ).I l £W_»lTl b2w . 

n n 


a 0 = exp 


X ln (v,) _ « 11 -T 


(10.24) 


(10.25) 


Beispiele 

Im Folgenden werden die Meßwerte aus Abschnitt 10.1 durch die 
genannten Funktionstypen angenähert und die dazugehörigen Gütekrite¬ 
rien berechnet. 

Logarithmusfunktion 

J» (*,.) = 13,07 £ ln (*i) = 4,653 

2>.-= 8,6 £(ln(jc,)) 2 = 8,545 

y = (!>’,)/»= 1-433 

13,07 _ 4,653.8T, 

Gleichung (10.16): a, =- TZm — = 1,297 

8,545 - 4 '“ 3 
6 

Gleichung (10.17): a 0 = ~ - —“ 97 . 4,653 = 0,428 
6 6 

SQMY = 8,30 SQAY = 10,51 

8 30 

Gleichung (10.9): R 2 = ^ = 0,79 

I -/(x,)) 2 = 2,225 
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Gleichung (10.10): 


Gleichung (10.11): 


dy = 


dyrel = 



0,609 

1,433 


= 0,609 


100% = 42,5% 


Die Gleichung lautet also y = 0,428 + 1,297 ln (x), die Anpassung ist 
jedoch völlig unzureichend. 


£>,.= 18,5 
£.v f = 88,25 
SQMY = 22,01 


Exponentialfunktion 


£lnta)= "1.378 
= 8 ^ 

SQAY = 10,51 


X A 'i ln (T’i) = 11,32 
j? = 1.433 
I(v i -/(-v.)) 2 = 3,202 



11,32 - 18 ’ 5 (“ 1 ’ 378 ) 

Gleichung (10.20): 

«. =-77TT2-= 0,499 

88,25 - 18 #- 
6 

Gleichung (10.21): 


Gleichung (10.9): 


Gleichung (10.10): 


Gleichung (10.11): 

dyrel= \4K ' 100% = 51 % 


0,171 


Die Gleichung _y = 0,171 exp(0,499.v) ist also schon besser als die vorher 
berechnete, aber immer noch wesentlich schlechter als die lineare 
Regression. 
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Potenzfunktion 


£ (ln (*,.)) (ln 00) = 5,607 
£ln(jc ( ) = 4,654 
I>i = 8,6 


£>00 = 
I(lnW ) 2 = 
y = 


-1,378 

8,545 

1,433 


5,607 - 


-1,378-4,654 


Gleichung (10.24): 


a, = 


8,545 - 


4,654 


= 1,353 


Gleichung (10.25): a 0 = expf —- - 1,353 

\ 6 6 

SQMY = 10,62 SQAY= 10,52 


= 0,278 


Gleichung (10.9): 


, 2 10,62 


R 2 = 


10,52 


= 1,01 


I(* -f(*i)) 2 = 0,296 

Gleichung (10.10): dy = — 0,22 

Gleichung (10.11): dyrel = ' 100% = 15,5% 

Funktionsgleichung: y = 0,278.x 1 - 353 

Mit der Anpassung durch eine Potenzfunktion erreicht man also schon 
eine fast optimale Annäherung an die Meßwerte. 


Literaturh in weise 

G, W: [14], S. 80 ff; G,W: [18]; G,W: [22], S. 340 bis 343 


Programmbeschreibung 

Da mit dem Programm Ausgleich vier verschiedene Anpas¬ 
sungsmöglichkeiten gegeben sind, für die jeweils die 
Gütekriterien bestimmt werden, hat auch dieses Programm 
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Zeile 30 

Zeile 40 

bis 130 

Zeile 130 

Zeile 140 

bis 250 


Zeilen 260 
bis 290 


Zeilen 300 
bis 330 

Zeilen 400 
bis 540 


Zeilen 1700 
bis 1880 


Zeilen 430 
bis 510 

Zeilen 2500 
bis 2810 
Zeilen 2000 
bis 2100 


eine Länge von ca. 8 K. Nach Dimensionierung der Felder 

für die x- und j'-Werte in 

deren Größe auf 100 festgesetzt ist, folgt in 

ein Menü, über das die gewünschte Anpassungsart ausge¬ 
wählt wird. Entsprechend der Codezahl verzweigt 
in die notwendigen Unterprogramme. Nach erfolgter 
Ausgleichung hat man über 

zunächst die Möglichkeit, mit der erhaltenen Gleichung 
Zwischenwerte zu berechnen; diese werden einzeln auf dem 
Bildschirm ausgegeben. In den 

ist die Möglichkeit des Rücksprungs ins Menü gegeben, 
damit für die gleichen Meßwerte eine neue Anpassungsart 
gewählt werden kann. Diese Möglichkeit wird man nutzen, 
wenn die Gütekriterien aussagen, daß die Anpassung nicht 
gut genug war. Mittels des 

kann das Programm schließlich neu gestartet werden. In 
den 

ist das Unterprogramm für die lineare Regression abgelegt. 
Es benutzt wie alle anderen entsprechenden Programmteile 
die 

zur Eingabe der Meßwerte; in der Variablenanza/z/ wird die 
Zahl der Meßwerte festgehalten. 

Die Berechnung der Regressionsfaktoren erfolgt dann in 
den 

Die nächsten beiden Schritte sind dann wieder gleich für alle 
Anpassungsarten; nämlich die Benutzung der 

für die Berechnung der Gütekriterien und der 

für die Ausgabe der Ergebnisse. Prinzipiell der gleiche 
Vorgang läuft ab in den 
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Zeilen 600 

bis 820 für die Anpassung durch eine Exponentialfunktion, in den 

Zeilen 1000 

bis 1220 für die Anpassung durch eine Logarithmusfunktion sowie in 
den 

Zeilen 1300 

bis 1520 für die Anpassung durch eine Potenzfunktion. In allen drei 
Teilprogrammen kommt jedoch hinzu, daß vor Durchfüh¬ 
rung der Rechnungen geprüft werden muß, ob die gewählte 
Anpassung überhaupt möglich ist (s. Grundlagen). Wenn 
das nicht der Fall sein sollte, erfolgt die entsprechende 
Ausgabe und der Rücksprung ins Menü, damit eine andere 
Anpassungsart gewählt werden kann. 


10 REM Programm Ausgleich 
20 MODE 2 

30 DIM x(100),y(100) 

40 CLS : REM Menue 

50 PRINT"Mlt welchem Funktionstyp sollen die Messwer 
te angepasst werden ?":PRINT"_____ 


60 PRINT:PRINT“- Gerade (lineare Regression) y = 

a + b*x ”;TAB(70); ’ (1) * 

70 PRINT:PRINT"- Exponentialfunktion y = 

a*exp(b*x)“> TAB(70);"< 2)“ 

80 PRINT:PRINT"- Logarithmusfunktion y = 

a + b*log<x)"J TAB(70);“(3)" 

90 PRINT:PRINT’- Potenzfunktion y = 

a*x A b (a>0)’;TAB(70);"(4)" 


100 PRINT:PRINT"- mit gar keiner";TAB(70)5"(5)" 
110 INPUT was 

120 IF was>5 OR was<l THEN 40 
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130 ON was GOSUB 400,600,1000,1300,1600 

140 CLS: PRINT“Soll die ermittelte Gleichung fuer In 

terpolatIonen verwendet werden < /n)“»:INPUT a* 

150 IF a$="“ THEN 170 

160 IF a$="n“ THEN 260 ELSE 140 

170 PRINT"Beenden Sie die Interpolationsrechnungen d 
urch Eingabe von x = 333.33“ 

180 PRINT:PRINT“x : ";:INPUT x 
190 IF was = 1 THEN y = a + b*x 

200 IF was = 2 THEN y = a*EXP(b*x) 

210 IF was = 3 THEN y = a + b*L0G<x> 

220 IF was = 4 THEN y = a*x A b 

230 IF x = 333.33 THEN 260 
240 PRINT"y = "Jy 
250 GOTO 180 
260 weiter = 0 

270 CLS : PRINT'Sollen die gleichen Messwerte fuer e 
Ine erneute Anpassung verwendet werden < /n)":INPUT 
a$ 

280 IF a*="" THEN weiter =1 : GOTO 40 
290 IF a*=“n“ THEN 300 ELSE 190 
300 CLS 

310 PRINT’Soll eine weitere Anpassung mit anderen Me 
sswerten durchgefuehrt werden < /n>“5:INPUT a$ 

320 IF a$=““ THEN RUN 

330 IF a*="n“ THEN 1600 ELSE 300 

340 END 

400 REM UP/HP lineare Regression 

410 al$ = "Lineare Regression y = a + b*x" 

420 GOSUB 1700 : REM Eingabe der Messwerte 
430 suml=0 : sum2=0 : sum3=0 : sum4=0 
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440 FOR k=1 TO anzahl 
450 suml=suml + x(k)*y(k> 

460 sum2=3um2 + x(k) 

470 sum3=sum3 + y(k) 

480 sum4=sum4 + x<k)*x(k) 

490 NEXT k 

500 b = (suml - <3um2*sum3/anzahl))/<sum4 - <sum2*su 
m2/anzahl)) 

510 a = sum3/anzahl - <b*sum2/anzahl> 

520 GOSUB 2500 : REM Guetekr1terlen 
530 GOSUB 2000 : REM Ausgabe 
540 RETURN 

600 REM Anpssung durch Exponentialfunktion 
610 GOSUB 1700 : REM Eingabe der Messwerte 
620 al* = “Exponentlalanpassung y = a*exp(b*x>“ 

630 suml=0 : sum2=0 : sum3=0 : sum4=0 
640 FOR k=l TO anzahl 
650 IF y<k)<=0 THEN 760 
660 sural=suml + x(k) 

670 sum2=sum2 + L0G(y(k>) 

680 sum3=sum3 + x(k)*L0G<y(k)) 

690 sum4=sum4 + x(k)*x(k) 

700 NEXT k 

710 b = (sum3 - (sum1*sum2/anzahl)) / (sum4 - (suml* 
sum1/anzahl)) 

720 a = EXP<(sum2/anzahl) - (b*sum1/anzahl)) 

730 GOSUB 2500 : REM Guetekr1terlen 
740 GOSUB 2000 : REM Ausgabe 
750 RETURN 
760 CLS 
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770 PRINT'Kurvenanpassungen ueber eine Exponentlalfu 
nktion sind nur moeglich, wenn alle y(l) groesser 
0 sind“ 

780 PRINT:PRINT’Hler ist aber y’;k;" = “;y(k> 

790 LOCATE 1,23 

800 PRINT"We1ter mit beliebiger Taste" 

810 IF INKEY$="" THEN 810 
820 RETURN 

1000 REM Kurvenanpassung per Logarithmusfunktion 
1010 GOSUB 1700 : REM Messpunkteingabe 
1020 al* = "Logarlthmlsche Funktionsanpassung y 

= a + b*log(x) 

1030 suml=0 : sum2=0 : sum3 = 0 : sum4=0 
1040 FOR k=1 TO anzahl 
1050 IF x<k)<= 0 THEN 1160 
1060 suml=suml + y<k)*LOG(x(k)> 

1070 sum2=sum2 + L0G(x<k)> 

1080 sum3=sum3 + LOG(x(k))*LOG(xCk)) 

1090 sum4=sum4 + y(k) 

1100 NEXT k 

1110 b = (suml - <sum2*sum4/anzahl)> / <sum3 - <sum2 
*sum2/anzahl)) 

1120 a = <sum4 - <b*sum2)>/anzahl 

1130 GOSUB 2500 : REM Guetekr1ter len 

1140 GOSUB 2000 : REM Ausgabe 

1150 RETURN 

1160 REM x<=0 

1170 CLS 

1180 PRINT’Elne 1ogar1thm1 sehe Kurvenanpassung Ist n 
lcht moeglich, wenn ein x-Wert kleiner oder 

gleich 0 Ist. Hier ist aber 
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x";k;" = ";x(k> 

1190 LOCATE 1,23 

1200 PRINT’Welter mit beliebiger Taste" 

1210 IF INKEY$="" THEN 1210 
1220 RETURN 

1300 REM Anpassung ueber eine Potenzfunktion 
1310 al* = "Anpassung ueber die Potenzfunktion 
y = a*x A b" 

1320 GOSUB 1700 : REM Messwerteingabe 
1330 suml=0 : sum2=0 : sum3=0 : sum4=0 
1340 FOR k=l TO anzahl 
1350 IF x<kX=0 OR y<kX = 0 THEN 1460 
1360 sum1=sum1 + L0G<x(k))*L0G(y<k>) 

1370 sum2=sum2 + L0G(x(k)> 

1380 sum3=sum3 + L0G<y<k)) 

1390 sum4=sum4 + L0G<x<k))*L0G<x(k)) 

1400 NEXT k 

1410 b = (suml - <sum2*sum3/anzahl)) / <sum4 - <sum2 
*sum2/anzahl)) 

1420 a = EXP(sum3/anzahl - (b*sum2/anzahl)) 

1430 GOSUB 2500 : REM Guetekrlterlen 
1440 GOSUB 2000 : REM Ausgabe 
1450 RETURN 
1460 CLS 

1470 PRINT'Bel einer Anpassung ueber eine Potenzfunk 
tlon muessen alle x(l) und alle yC1> groesser 0 se 
ln." 

1480 PRINT:PRINT“Hier Ist aber x";k;" = ";x(k>;" un 
d y";k;* = "}y<k> 

1490 LOCATE 1,23 

1500 PRINT"Weiter mit beliebiger Taste" 
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1510 IF INKEY*="“ THEN 1510 
1520 RETURN 

1600 REM Schlussprogramm 
1610 CLS 

1620 LOCATE 30,10 
1630 PRINT"Tschuess" 

1640 LOCATE 1,23 
1650 END 

1700 REM Eingabe der Messwerte 
1710 IF we1ter=1 THEN RETURN 
1720 anzahl=0 
1730 CLS 

1740 PRINT"Geben Sie die Messwerte eln.":PRINT:PRINT 

"Beenden Sie die Messwerteingabe durch Eingabe von 3 

/ 

33.33 als x-Wert" 

1750 LOCATE 1,6 
1760 anzahl=anzahl+l 

1770 PRINT“x"JanzahlINPUT x(anzahl) 

1780 IF x<anzahl)=333.33 THEN 1830 

1790 PRINT:PRINT“y“;anzahlINPUT y(anzahl) 

1800 LOCATE 1,6 : PRINT" 


1810 LOCATE 1,8 : PRINT" 


1820 GOTO 1750 

1830 anzahl=anzahl- 1 

1840 LOCATE 1,15 

1850 PRINT'Alle Eingaben i.o. 

1860 IF a*="" THEN 1880 


( /n) ";:INPUT a* 
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1870 IF a$="n" THEN 1700 ELSE 1840 
1880 RETURN 

2000 REM Ausgabeunterprogramm 

2010 CLS 

2020 PRINT al* 

2030 PRINT" 


2040 PRINT'a : ";a,,"b : ";b 
2050 PRINT:PRINT“r"‘2 : ";r2u 

2060 PRINT:PRINT"Abwelchung : ";abw 

2070 PRINT" = ";USING"###.##";ab 

wrel;:PRINT"\" 

2080 LOCATE 1,23 : PRINT’Welter mit beliebiger Taste 

M 

2090 IF INKEY*=“" THEN 2090 
2100 RETURN 

2500 REM Berechnung der Guetekrlterlen 
2510 REM Relative Quadratische Abweichungen 
2520 abw=0 

2530 FOR k=l TO anzahl 
2540 x=x(k) 

2550 IF was=l THEN yb = a + b*x 

2560 IF was=2 THEN yb = a * EXP(b*x) 

2570 IF was=3 THEN yb = a + b*LOG(x) 

2580 IF was=4 THEN yb = a * x*b 

2590 abw=abw + (yb-y(k)>*<yb-y(k)) 

2600 NEXT k 

2610 abw = SQRCabw/anzahl) 

2620 REM Best 1mmthe1tsmass nach Ullmann 
2630 sqay = 0 : sqmy = 0 
2640 suml=0 
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2650 FOR k=1 TO anzahl 
2660 3uml=suml+y<k) 

2670 NEXT k 

2680 mwy=suml/anzahl 

2690 FOR k=l TO anzahl 

2700 sqay = sqay + (y<k)-mwy)*<y<k)-mwy) : REM Quadr 
atIsche Abweichungen der Messwerte 
2710 x=x(k) 

2720 IF was=l THEN yb = a + b*x 

2730 IF was=2 THEN yb = a * EXP<b*x> 

2740 IF was=3 THEN yb = a + b*L0G<x> 

2750 IF was=4 THEN yb = a * x Ä b 

2760 sqmy = sqmy + <yb-mwy)*(yb-mwy) : REM Quadrat 
Ische Abweichungen der berechneten y-Werte 
2770 NEXT k 
2780 r2u = sqmy/sqay 

2790 REM relative Abweichung, bezogen auf den Mittel 
wert 

2800 abwrel = abw/mwy*100 
2810 RETURN 


10.3 Anpassung durch ein Polynom 

Grundlagen 

Wenn die Anpassung einer Meßwertreihe durch eine einfache Funktion 
nicht oder nur ungenügend gelingt, wird man als Kurventyp für die 
Anpassung ein Polynom wählen. Solch ein Polynom hat die allgemeine 
Funktionsgleichung 

y = a 0 + a x x 1 + a 2 x 2 + ... + a n x" 

n gibt hier den Grad des Polynoms an. 


( 10 . 26 ) 
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Die Ausgleichsforderung lautet wieder 

N 

[vv] = X U( x i) - yd = Minimum! 

i = 0 

V’, ist der Meßwert und / (a,) ist der berechnete Funktionswert, der jeweils 
zu x t gehört. Wenn z die Meßwertanzahl bezeichnet, läuft / im Gegensatz 
zu den beiden vorangegangenen Abschnitten nicht mehr von 1 bis z, 
sondern von 0 bis N, wobei N = z - 1 ist. Die Anzahl der Summations¬ 
glieder bleibt also gleich der Meßwertanzahl. 

Ausführlich geschrieben lautet die Ausgleichsforderung: 

N 

[vv]= £ Oo + «1 *,• + •• • +a„x?-y i ) 2 = Minimum! (10.27) 
; = o 

Gleichung (10.27) muß jetzt wie in Abschnitt 10.1 gezeigt partiell nach 
allen Variablen a 0 bis a„ abgeleitet werden; die so erhaltenen Gleichungen 
werden alle gleich 0 gesetzt, um die Bedingung für ein Minimum zu 
erfüllen. 

N 

2 £ (fl 0 + + ... +<*!*? - y t ) =0 

i = 0 
N 

2 I K + ö, *,■ + • • • + a, -v," - >',)*,. = 0 
1 = 0 

d[v v] i 

~rr~ = 2 I (a 0 + a, x, + ... + a n xr - yi )x" = 0 

Cia 0 i = 0 

Dieses Gleichungssystem kann man auch folgendermaßen schreiben (alle 
Summen laufen hier von 0 bis N): 

(TV + \)a 0 4- a t YjXi + ... + «, 

«oI-V + öi I-V 2 + • • • + a„2>r = X>', V 

öoI-y" + ö, Xv n+1 + ... + ö„£.vr +n = 

Aus diesem Gleichungssystem kann man jetzt die Koeffizienten a 0 bis a n 
berechnen und so das optimale Polynom aufstellen. 

In Kurzform geschrieben lautet das Gleichungssystem: 

1 a k ( £ *?•*/)= I yi-xl 0 = 0,1,...,«) (10.28) 

fc = 0 \i = 0 / i — 0 


d[vv] 

da 0 

d[vv] 

da 0 
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Beispiel 

Für die Meßwertreihe aus Abschnitt 10.1 soll das optimale Polynom 
1. Grades bestimmt werden. 

Meßwertanzahl = 6: N = 5 
n = 1 (Grad des Polynoms) 

Das Gleichungssystem lautet also (die Summen laufen von 0 bis 5): 
j = 0: 6« 0 + a, £>,- = X v f 

7 = 1: a 0 X-v, + «, X-V , 2 = X>'. *. 

Mit X^i = 18,5; X.v,- = 8,6; X*f = 88,25; Xji-^i = 44,25 ergibt sich 

6a 0 + 18,5a, = 8,6 
18,5« 0 + 88,25a, =44,25 

Hieraus errechnen sich a 0 und a l : 

a 0 = -0,318 
a, = 0,568 

Diese Werte stimmen mit den in Abschnitt 10.1 gefundenen überein; die 
Gütekriterien sind also 

R 2 = 0,959; dy = 0,270; dyrel = 18,9 % 

Wählt man als Ausgleichsfunktion ein Polynom höheren Grades, steigt 
die Anpassungsgenauigkeit; mit n = 5 erreicht man so z. B. ein Bestimmt¬ 
heitsmaß von 1,0. 


Literaturhinweise 

G, W: [14], S. 80 ff; G,W: [18]; G,W: [22], S.340 bis 343 
Program mbesch reibung 

Die Zeilennummern des Programms Polanpass beginnen erst bei 10000; 
auf diese Art und Weise soll erreicht werden, daß dieses Programm 
leichter mit dem zuvor beschriebenen Programm Ausgleich zusammenge¬ 
faßt werden kann, so daß alle Ausgleichsprogramme dann in einem 
Programmpaket zusammengefaßt sind. 
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In 

Zeile 10030 werden über das Unterprogramm in den 
Zeilen 10200 

bis 10310 die Ausgangswerte eingegeben. Entsprechend den Grund¬ 
lagen wird dann in den 

Zeilen 10500 

bis 10800 das Gleichungssystem aufgestellt, das der Ausgleichsforde¬ 

rung genügt. Als Ergebnis hat man die Lösungsmatrix 
Imat. Die 

Zeilen 10590 
bis 10610 

und 

Zeilen 10630 

bis 10650 berechnen lediglich xf bzw. x \; die etwas komplizierte 
Ausführung wurde gewählt, weil der CPC 464 wie schon 
erwähnt Schwierigkeiten beim Potenzieren negativer Zah¬ 
len hat; das Vorzeichen ist hier nicht immer richtig. In den 

Zeilen 11000 

bis 11200 werden die einzelnen Polynomglieder berechnet. Da es sich 
hierbei um die Lösung eines linearen Gleichungssystems 
handelt, wurde das entsprechende Programm aus Kapitel 8 
übernommen. Anschließend werden in den 

Zeilen 11400 

bis 11590 wie in Abschnitt 10.2 die Gütekriterien berechnet. Die 
Ergebnisse werden über die 

Zeilen 11800 

bis 11950 ausgegeben. Zwischen den 
Zeilen 12000 

und 12090 sind die gleichen Weiterführungsmöglichkeiten wie im 
Programm Ausgleich gegeben. 

Programmlänge: etwa 4,4 K. 

Die Güte des verwendeten Lösungsalgorithmus kann man einfach 
dadurch überprüfen, daß man Werte einer beliebigen Funktion dem 
Programm als Meßwerte eingibt. Wählt man dann den Polynomgrad 
mindestens so hoch wie den tatsächlichen Polynomgrad, so findet das 
Programm normalerweise die richtigen Faktoren. 
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10000 REM Programm Polanpass 1.0 

10010 DIM lmat <10, 10), x<10),y<10),al(10),a<10) 

10020 MODE 2 

10030 GOSUB 10200 : REM Eingabe der Wertetabelle 

10040 GOSUB 10500 : REM Loesungsmatr1xberechnung 

10050 GOSUB 11000 : REM Berechnung der Polynomgllede 
r 

10060 GOSUB 11400 : REM Guetekrlterlen 
10070 GOSUB 11800 : REM Ausgabe 
10080 GOSUB 12000 : REM We1terfuehrungsprogramm 
10090 CLS : LOCATE 20,10 : PRINT“Na endlich. Das wir 
d auch Zelt":PRINT:PRINT TAB<23>;“Ich bin schon ganz 
muede . “ 

10100 LOCATE 1,23 
10110 END 

10200 REM Eingabe der Wertetabelle 

10210 CLS : IF weiter =1 THEN RETURN 

10220 PRINT “Wieviele Wertepaare sollen eingegeben w 

erden“J:INPUT anz 

10230 ERASE x : ERASE y : DIM x(anz),y<anz) 

10240 FOR k = 0 TO anz-1 

10250 PRINT “x“»k+l;: INPUT x<k) 

10260 PRINT"y“Jk+l;"=“;:INPUT y(k) 

10270 NEXT k 

10280 PRINT:PRINT“Eingabe ok < /n)“;:INPUT a$ 

10290 IF a*=“n“THEN 10240 
10300 IF a$<>““THEN 10280 
10310 RETURN 

10500 REM Loesungsmatrixberechnung 
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10510 CLS 

10520 PRINT'Welchen Grad soll das Ausgleichspolynom 
haben“;:INPUT grad 

10530 PRINT:PRINT:PRINT"Bltte warten Sie. Ich bin be 
schaeftIgt.“ 

10540 ERASE lmat :DIM lmat(grad,grad+l> 

10550 FOR j =0 TO grad 

10560 FOR k=0 TO grad 

10570 FOR 1=0 TO anz-1 

10580 IF x< 1 ><0 THEN 10590 ELSE 10610 

10590 IF k/2=FIX<k/2) THEN 10610 

10600 xl=-ABS<x< lKk) ; GOTO 10620 

10610 x1= ABS(x(1>*k> 

10620 IF 1(1X0 THEN 10630 ELSE 10650 
10630 IF j/2=FIX(j/2> THEN 10650 
10640 x2=-ABS<x<1) Ä J) : GOTO 10660 
10650 x2= ABS<x<lKJ) : 

10660 lmat<J,k) = lmat<j,k) + xl * x2 
10670 NEXT 1 
10680 NEXT k 
10690 NEXT J 

10700 REM b-Wert-Berechnung 

10710 FOR J =0 TO grad 

10720 FOR 1=0 TO anz-1 

10730 IF 1(1X0 THEN 10740 ELSE 10760 

10740 IF J/2=FIX< J/2) THEN 10760 

10750 x2=-ABS<x<1)~J) : GOTO 10770 

10760 x2= ABS(x(1 ) 

10770 lmat<J,grad+l) = 1mat<J,grad+1) + y(l)*x2 
10780 NEXT 1 
10790 NEXT i 
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10800 RETURN 

11000 REM Determinantenberechnung nach Gauss 

11010 ERASE a : ERASE al : DIM a<grad,grad+1>,al<gra 

d) 

11020 FOR 1=0 TO grad 
11030 FOR k = 0 TO grad+1 
11040 a<1,k)=lmat<1,k) 

11050 NEXT k 

11060 NEXT 1 

11070 FOR k = 0 TO grad 

11080 FOR 1 = 0 TO grad 

11090 IF 1=k THEN GOTO 11150 

11100 f = -a(l,k)/a(k,k) 

11110 FOR 1 = 0 TO grad+1 

11120 b= a(k,1)*f 

11130 a(1 , 1>=a<1,1) + b 

11140 NEXT 1 

11150 NEXT 1 

11160 NEXT k 

11170 FOR 1=0 TO grad 

11180 a 1 (1) = a(l,grad+1) / ad,i> 

11190 NEXT 1 
11200 RETURN 

11400 REM Guetekrlterlen 

11410 mwy = 0 : sqmy = 0 : sqay = 0 

1 1420 FOR k = 0 TO anz-1 

11430 mwy = mwy + y<k) 

11440 NEXT k 
1 1450 mwy = mwy/anz 
11460 FOR k=0 TO anz-1 

11470 sqay = sqay + <y(k)-mwy)*(y<k)-mwy> 
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11480 x=x<k> : GOSUB 11600 : REM Funktlonswertberech 
nung 

11490 sqmy = sqmy + <yb-mwy>*<yb-mwy) 

11500 NEXT k 

11510 r2u = sqmy/sqay 

11520 REM Quadratische Abweichungen 

11530 FOR k=0 TO anz-1 

11540 x = x(k> : GOSUB 11600 : REM Funkt 1onswertberech 
nung 

11550 abw = abw + (y(k)-yb)*(y(k)-yb> 

11560 NEXT k 

11570 abw = SQRCabw/anz) 

11580 abwrel = abw/mwy*100 
11590 RETURN 

11600 REM Funktionswertberechnung 

11610 yb = 0 

11620 FOR 1=0 TO grad 

11630 yb = yb + al(l)*x Ä l 

11640 NEXT 1 

11650 RETURN 

11800 REM Ausgabe 

11810 CLS 

11820 PRINT"Anpassung durch ein Polynom“ 

11830 PRINT:PRINT"y = aO + al*x + a2*x Ä 2 +...+an*x A n 

M 

11840 PRINT*_ 


11850 PRINT 
11860 FOR k=0 TO grad 
11870 PRINT'a";k;": ";al(k) 
11880 NEXT k 
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11890 PRINT:PRINT:PRINT"r Ä 2 : "ir2u 
11900 PRINT:PRINT"Mlttlere Abweichung : ";abw 
11910 PRINT" = "5USING"###.##"; 

abwrel;:PRINT*%“ 

11920 LOCATE 1,23 

11930 PRINT”Welter mit beliebiger Taste" 

11940 IF INKEY*="* THEN 11940 
11950 RETURN 

12000 REM We1terfuehrungsprogramm 
12010 CLS 

12020 PRINT'Sollen die gleichen Messwerte fuer eine 
weitere Berechnung verwendet werden 

( /n) 

";:INPUT a$ 

12030 IF a*="" THEN weiter =1 : GOTO 10020 
12040 IF a$< >"n" THEN 12000 
12050 CLS 

12060 PRINT'Soll eine weitere Rechnung mit anderen M 
esswerten durchgefuehrt werden < /n> INPUT a$ 
12070 IF a*="" THEN RUN 10000 
12080 IF a*<>"n" THEN 12050 
12090 RETURN 
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folg anhand der Knobeltabellen 
festzustellen! 


Ein praktischer Kurs auf zwei Ebe¬ 
nen mit Beispielen und Lösungs¬ 
wegen für Schulen/Hochschulen, 
Aus-/Weiterbildung und für Hobby¬ 
programmierer. Mit jedem der ins¬ 
gesamt 20 Programme werden 
neue BASIC-Anweisungen einge¬ 
führt. An jedes Programm schlie¬ 
ßen sich zehn Übungen an, die das 
Verständnis für Programmstruktu¬ 
ren und Anweisungen vertiefen. 
Anfängerkenntnisse sind vorteilhaft. 


Merkel, Erich 
BAS IC-Intensiv¬ 
kurs I 

Sprachelemente, 
Strukturen, 
Programmaufbau 
Reihe CHIP WISSEN 
256 Seiten, 

25,- DM, 1985 
ISBN 3-8023-0775-5 


Merkel, Erich 
BAS IC-Intensiv¬ 
kurs II 

Massenspeicher, 
Drucker, Grafik, 
komplexe Strukturen 
Reihe CHIP WISSEN 
ca. 260 Seiten, 
ca. 28,- DM, 1985 
ISBN 3-8023-0869-7 
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Pol, Bernd 
Wie man in 
BASIC 

programmiert 

Reihe CHIP WISSEN 
368 Seiten, 

16 Abbildungen, 

3. Auflage 1984 
30,- DM 

ISBN 3-8023 0637 6 


Ein Buch für Praktiker, und mehr 
als nur eine Einführung! An zwei 
bis ins Detail ausgearbeiteten Fall¬ 
studien werden die Grundlagen des 
Programmierens verdeutlicht und 
die wichtigsten BASIC-Bestandteile 
eingehend besprochen. Vor allem: 
Wie ist ein Problem zu lösen? 
Warum ist das so formuliert? Wie 
wendet man Programmiertechniken 
mit BASIC an? Diese und ähnliche 
Fragen werden beantwortet. 


Schneller erfolgreich durch Computer-Bücher 







Mein Home-Computer 
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Und 


das bringt 


jeden Monat: 

* Programme für alle gängigen Home-Computer 

* Anwendungsbeispiele aus der Praxis 

$ Marktübersicht, Tests und Kaufberatung für Zusatzgeräte 
und Home-Computer 

* Schnellkurse für Einsteiger zum Sammeln 

* Tips und Tricks 

* Interessantes, Aktuelles und Unterhaltsames aus der 
Home-Computer-Szene 

* News, Clubnachrichten 

Holen Sie sich die neueste Ausgabe bei Ihrem Zeitschriften¬ 
händler oder fordern Sie ein Kennenlernheft direkt beim Vogel- 
Verlag, Leserservice HC, Postfach 67 40, 8700 Würzburg, an. 












CHIP Das Mikrocomputer-Magazin 
Jeden Monat umfassende Informationen — 
für alle, die Bescheid wissen müssen. 

CHIP ist das Mikrocomputer-Magazin, das über 
die ganze Welt der Mikrocomputer berichtet. 
Mit CHIP erfahren Sie mehr über die Arbeits¬ 
weise eines Mikrocomputers, über die 
Anwendungsmöglichkeiten und Leistungs¬ 
merkmale der Geräte. 



Mehr noch, Monat für Monat bringt CHIP 
aktuell: 

— Erfahrungsberichte aus der Praxis 

— Problemlösungen 

— Trends 

— Hard- und Soft¬ 
ware-Tests 

— Marktübersichten 

— Kaufberatung 

— und die CHIP-Börse 


CHIP gibt es für 6,50 DM monatlich bei Ihrem Zeit¬ 
schriftenhändler. Im Abonnement für 69,— DM pro Jahr 
(statt 78,— DM). Oder direkt beim CHIP- Leser-Service, 
Vogel-Verlag, Postfach 67 40. 8700 Würzburg 1. 





Mit diesem Buch entdecken Sie immer wieder 
neue, gute Seiten Ihres CPC 464 und überwinden 
den Frust vor mathematischen Problemen. 

Ob Schüler, Student, «fertiger« Techniker oder 
Ingenieur: Hier finden Sie die richtige 
Problemauswahl, verständliche und eindeutige 
Erläuterungen der Lösungswege, durchgerechnete 
Beispiele (ohne Benutzung der Programme) und — 
als Schwerpunkt - die Programme selbst. Damit 
können Sie natürlich auch unabhängig von den 
Grundlagenabschnitten arbeiten. Der Umgang 
mit den Programmen ist recht einfach, auch 
eventuelle Änderungen machen keine 
Schwierigkeiten dank der ausführlichen 
Beschreibungen. 



VOGEL- BUCHVERLAG 
WÜRZBURG 


ISBN 3-8023-0856-5 


Ein Buch von 


CHIP 


dem 


Mikrocomputer-Magazin 
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